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FOREWORD* 


The  Workshop  on  Directional  Acoustic  Sensors  provided  a  forum  for  the  technical 
exchange  of  ideas  and  experience  in  the  use  of  directional  acoustic  sensors  for  undersea 
applications.  Researchers  from  academia,  industry,  and  government  research  laboratories 
presented  their  current  work  on  the  theory,  design,  and  application  of  various  types  of  directional 
acoustic  sensors.  This  workshop,  which  was  held  in  Newport,  Rhode  Island,  builds  upon  the 
1995  Workshop  on  Acoustic  Particle  Velocity  Sensorst  by  addressing  progress  and  new  theories 
in  the  area  of  multi-component  directional  acoustic  sensors. 

Twelve  invited  papers  were  presented  during  the  two-day  workshop.  A  poster  session  was 
also  conducted  in  which  various  organizations  presented  overviews  of  specific  directional  sensor 
applications.  Spirited  discussions  explored  issues  important  in  the  development  and  application 
of  directional  sensors.  Several  of  the  invited  papers  pointed  out  the  Taylor  expansion 
representation  of  the  acoustic  field,  the  first  three  terms  of  which  represent  the  scalar,  vector,  and 
tensor  (or  dyadic)  components  of  the  field,  referred  to  by  some  as  zero-,  first-,  and  second-order 
tensor  components.  Topics  ranged  from  system  noise  characteristics  of  sensors  measuring  these 
components  to  the  signal  processing  of  their  outputs.  Studies  on  the  mechanical  design  and 
mounting  of  various  forms  of  these  sensors  were  described.  Some  designs  comprised  acoustic 
motion-sensors,  whereas  others  are  essentially  fixed-sensors  that  measure  gradients  of  various 
order  acoustic  components.  Other  papers  provided  results  of  current  and  past  research  on  the  use 
of  these  sensors  in  various  applications:  ocean  acoustics;  hull-mounted  and  towed-sonar  arrays; 
and  deployed  acoustic  sensors  (i.e.,  DIFAR  systems).  Development  plans  and  application  results 
of  forming  certain  combinations  of  these  acoustic  sensor  orders  were  also  presented,  including 
the  use  of  acoustic  intensity  probes  in  sound  scattering  measurements.  New  transduction 
materials  (single  crystals)  were  also  described,  showing  promise  for  advancing  the  performance 
and  design  of  vector  and  higher-order  acoustic  sensors. 

The  Office  of  Naval  Research  (ONR)  and  the  Naval  Undersea  Warfare  Center  (NUWC) 
Division  sponsored  this  workshop.  The  organizers  would  like  to  thank  Dr.  James  McEachem, 
ONR  Team  Leader  for  Sensors,  Sources,  and  Arrays,  and  Dr.  Richard  Nadolink,  NUWC  Chief 
Technology  Officer,  for  their  encouragement  and  financial  support.  We  also  thank  Ms.  Kathy 
Stark  of  America  House  Communications,  Newport,  RI,  for  the  planning  and  administration  of 
this  meeting. 


*The  following  proceedings  papers  were  published  “as  is”  and  were  not  reviewed  or  edited  by  the  Workshop 
Committee. 


t“Acoustic  Particle  Velocity  Sensors:  Design,  Performance,  and  Applications,”  A  IP  Conference  Proceedings  368, 
Mystic,  Connecticut,  1995,  M.  J.  Berliner  and  J.  F.  Lindberg  (eds.),  American  Institute  of  Physics  Press,  New  York, 
1996. 


It  is  hoped  that  these  proceedings*  will  stimulate  further  research  in  this  engaging  area  of 
underwater  acoustics.  Directional  acoustic  sensors  show  considerable  potential  for  enhancing 
undersea  acoustic  measurements  in  a  variety  of  commercial  and  military  applications.  Much 
work  remains  to  realize  the  full  promise  of  these  versatile  sensors. 


Dr.  Benjamin  A.  Cray,  NUWC 
Workshop  Co-Chairman 


'  !<?„  C,£/cU 

Dr.  Roy  C.  Els  wick,  ONR 
Workshop  Co-Chairman 


A  workshop  speaker,  Dr.  Boris  Aronov,  Chief  Scientist  at 
BTech  Acoustics,  Inc.,  and  Research  Professor  of  Acoustics  at 
the  University  of  Massachusetts,  Dartmouth,  describes  the 
features  of  a  directional  acoustic  sensor  component  currently 
being  evaluated  for  underwater  applications. 


Scientists  and  engineers  with  a  broad  range  of  experience  and  interests  discussed  the  design  and  use  of 
directional  acoustic  sensors  for  ocean  observations  and  organic  sensing  applications. 


*Bookmark  and  thumbnail  views  are  available  under  the  Window  menu  to  aid  the  electronic  reading  of  this  Adobe 
Acrobat  document.  (Please  note  that  this  document  was  prepared  for  screen  optimization.) 
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Interesting  Questions...  (we  hope ) 
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ABSTRACT 

Several  sensor  systems  that  simultaneously  measure  the  underwater  acoustic  pressure  and  particle 
velocity  fields  ("vector  sensors")  have  been  designed  and  deployed  at  sea  by  the  Marine  Physical  Labora¬ 
tory  over  the  past  few  decades.  These  systems  include  the  set  of  12  freely  drifting,  neutrally  buoyant  Swal¬ 
low  floats,  the  16-element  vertical  DIFAR  (actually  "TRIFAR")  array,  and  a  set  of  moored,  flow-shielded 
DIFAR  sonobuoy  systems  for  measuring  nearshore  acoustics.  In  addition,  standard  Navy  DIFAR  53D 
sonobuoys  have  been  deployed  in  many  experiments.  Examples  of  the  measurements  made  by  these  sensor 
systems  are  presented  to  illustrate  their  scientific  usefulness.  The  framework  for  the  analysis  of  the  data  is 
provided  both  by  the  energetics  of  acoustic  fields  and  by  array  processing  with  single- spatial-point  meas¬ 
urements.  Special  emphasis  here  is  placed  on  the  reactive  acoustic  intensity,  i.e.,  the  energy  flux  required 
to  support  the  spatial  structure  in  the  acoustic  field.  Ideas  for  extending  the  vector  sensor  concept  to  a  new 
"tensor  sensor"  also  are  presented. 

INTRODUCTION 

The  purpose  of  this  paper  is  to  summarize  work  done  at  the  Marine  Physical  Lab  (MPL)  in  direc¬ 
tional  acoustic  sensors.  Sect.  I  briefly  describes  some  of  the  ocean  acoustic  vector  sensor  systems 
designed,  built,  and/or  deployed  at  sea  by  MPL.  The  basic  equations  and  approach  used  in  the  analysis  of 
the  directional  sensor  data  are  presented  in  Sect.  II.  Sect.  Ill  contains  a  brief  overview  of  some  of  the 
scientific  results  obtained  with  the  use  of  vector  sensor  data.  At-sea  data  from  four  different  directional 
sensor  systems  deployed  in  four  different  ocean  acoustics  experiments  are  presented.  The  ideas  behind, 
and  the  benefits  of,  a  new  design  "tensor  sensor"  are  discussed  in  Sect.  IV.  We  hope  to  be  able  to  build 
such  a  system  in  the  near  future.  Finally,  a  few  concluding  remarks  are  given  in  Sect.  V. 

I.  OVERVIEW  OF  MPL  VECTOR  SENSOR  SYSTEMS 

A.  Freely  Drifting  Swallow  Floats 

The  Swallow  floats  are  autonomous,  neutrally  buoyant  sensor  systems  that  simultaneously  measure 
the  three  components  of  acoustic  particle  velocity  and  acoustic  pressure  in  the  0.3  to  20  Hz  frequency  band. 
They  drift  freely  along  with  the  prevailing  ocean  current  to  reduce  contamination  from  relative  flow  past 
the  sensor,  a  well-known  problem  with  measurements  in  the  infrasonic  band  from  fixed  sensor  systems. 
Fig.  1  shows  a  schematic  drawing  of  a  Swallow  float.  In  addition  to  the  infrasonic  acoustic  sensors  (i.e., 
the  infrasonic  hydrophone  below  the  glass  sphere,  the  3-component  geophone  assembly  on  the  underside  of 
the  north  pole  of  the  glass  sphere,  and  the  compass  inside  the  sphere  at  the  south  pole  used  to  orient  the 
horizontal  geophone  components),  each  float  was  equipped  with  an  8  kHz  hydrophone  system  to  provide 
acoustic  localization  capability.  Further  details  on  the  system  hardware  can  be  found  in  Ref.  1 . 

B.  The  Vertical  ’’DIFAR”  Array 

MPL’s  vertical  DIFAR  array,  re  Fig.  2,  is  composed  of  16  elements  with  an  inter-element  spacing  of 
15  m.  Each  of  the  elements  contains  three  orthogonally-oriented  geophones  to  measure  acoustic  particle 
velocity  and  a  hydrophone  to  measure  pressure.  Each  element  also  contains  a  flux-gate  compass  for 
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measuring  the  orientation  of  the  horizontal  geophones,  a  programmable  high-frequency  data  acquisition  for 
acoustic  element  localization,  a  preamp  with  programmable  gain  from  0  to  120  dB,  and  a  16-bit  A/D  con¬ 
verter  that  provides  additional  dynamic  range.  The  data  sampling  rate  is  programmable  to  either  150,  300, 
600,  1200,  or  2400  samples/sec.  Data  collected  to  date  were  sampled  at  600  samples/sec,  providing  an 
effective  data  bandwith  of  10-270  Hz.  Additional  details  on  the  vertical  DIFAR  array  hardware  are  con¬ 
tained  in  Ref.  2. 

C.  Nearshore  Moored  DIFAR  Sonobuoys 

As  part  of  MPL’s  Adaptive  Beach  Monitoring  program,  a  set  of  8  standard  Navy  SSQ  53D  DIFAR 
sonobuoys,  provided  by  the  Naval  Air  Warfare  Center,  were  modified  for  long-term  moored  deployments 
in  very  shallow  water  (less  than  50  m).  Their  purpose  was  to  measure  the  directional  properties  of  the 
nearshore  underwater  ambient  sound  field  as  a  function  of  range  offshore,  starting  just  outside  the  surf 
zone.  The  DIFAR  sensor  in  each  unit  was  suspended  by  an  elastic  rubber  tube  within  a  bottom-mounted, 
hemispherical  flow  shield.  The  multiplexed  sensor  signal  was  transmitted  up  an  electrical  cable  to  the  RF 
transmitter  mounted  on  a  surface  float  where  it  was  RF  telemetered  to  the  shore-based  recording  facility. 
The  surface  float  also  contained  sufficient  battery  power  for  continuous  3-5  day  operation.  System  design 
was  based  heavily  on  one  developed  by  researchers  at  the  Applied  Research  Lab,  University  of  Texas  at 
Austin  [3]. 

Another  sensor  system  whose  data  are  presented  in  Sect.  IIIC  is  a  variant  of  the  moored  DIFAR 
sonobuoy  system.  Referred  to  as  the  "OSZM"  (Outside- Surf-Zone  Mooring),  it  transmitted  its  data  by  elec¬ 
tromechanical  cable  deployed  through  the  surf  zone  to  a  land-based  recording  facility  rather  than  by  RF 
telemetry.  The  cable  also  allowed  electrical  power  to  be  transmitted  to  the  mooring,  thus  eliminating  the 
need  for  battery  replacement.  More  importantly,  instead  of  a  water-column  DIFAR  sensor,  its  sensor  pack¬ 
age  was  composed  of  three  orthogonally-oriented  geophones  and  a  hydrophone  mounted  on  a  flat  metal 
plate.  The  plate  was  diver-deployed  directly  on  the  ocean  bottom  just  outside  the  surf  zone.  Since  the  geo¬ 
phones  measured  the  three  components  of  sediment  motion  rather  than  those  of  the  water  column,  the  sys¬ 
tem  does  not  classify  as  a  directional  underwater  acoustic  sensor  per  se.  However,  the  unconsolidated  sand 
layer  covering  the  ocean  bottom  acts  nearly  like  a  fluid.  Therefore,  underwater  acoustic  signal  processing 
techniques  can  be  applied  to  these  data,  as  demonstrated  in  Sect.  IIIC. 

D.  Standard  U.S.  Navy  DIFAR  Sonobuoys 

A  final  sensor  system  whose  data  are  presented  in  Sect.  IIIB  is  the  standard  U.S.  Navy  DIFAR 
sonobuoy.  Not  only  are  these  systems  a  critical  part  of  the  Navy’s  air  ASW  effort  [4],  but  they  can  be  a 
valuable  tool  in  ocean  acoustics  scientific  studies  (e.g.,  [5]).  These  ingenious  sensor  systems  provide  a 
benchmark  against  which  newer  underwater  acoustic  directional  sensors  are  judged. 


II.  THEORETICAL  FRAMEWORK 


The  theoretical  basis  for  all  results  presented  in  this  paper  is  a  set  of  four  fundamental  equations. 
The  first  is  the  conservation  of  linear  momentum  under  the  linear  acoustic  approximation, 


dv  (x ,  t ) 

poW — ~ — 
3 1 


+  Vp(x,t)=0 


(1) 


The  ambient  density  is  p0(y)  and  the  acoustic  field  variables,  v(x,t)  and p(x,t ),  are  the  particle  velocity 
and  pressure,  respectively.  The  linearized  conservation  of  mass  with  the  equation  of  state  embedded  in  it  is 

dp  (x,t)  „ 

- “ - +Ks(x)  V-v(x,O  =  0  (2) 

3 1 


The  fluid  properties,  i.e.,  the  adiabatic  incompressibility  K^(y)  and  the  density,  are  assumed  to  be  fixed  in 
time  over  the  time  scales  of  acoustic  signal  propagation.  Deviations  from  an  ideal  fliud  (e.g.,  viscous  dissi¬ 
pation)  are  assumed  negligible.  These  two  equations  contain  all  the  physics  needed  for  this  paper.  They 
hold  in  regions  where  no  acoustic  sources  exist. 
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An  important  implication  of  Eq.  (1)  is  that  if  the  acoustic  pressure  field  is  measured  (or  otherwise 
known)  everywhere  within  a  given  region  of  space,  the  corresponding  acoustic  particle  velocity  field  can 
be  derived  everywhere  within  this  region.  That  is,  measurement  of  acoustic  particle  velocity  adds  no  new 
information  beyond  that  available  from  spatially  distributed  measurements  of  pressure. 


The  third  equation,  which  provides  the  basis  for  array  processing  with  measurements  at  a  single  point 
in  space,  is  an  expression  of  Taylor’s  theorem.  Applying  this  theorem  to  the  acoustic  pressure  field  in  a 
region  about  the  measurement  point,  x0 ,  where  no  acoustic  sources  exist,  gives 


p  (x,  t)  =  p  (x0 ,  t)  +  Vp  (x0 ,  t)  •  Ax  +  — Ax 


Matrix  of 
2nd  DerivativesJ 


Ax  + 


(3) 


where  Ax  =  x  -  xQ . 

Finally,  the  acoustic  strain  rate,  a  secondrank  tensor  denoted  by  [8  ]  and  which  is  used  in  the  concept 
of  a  "tensor  sensor"  presented  in  Sect.  IV,  is  equal  to 

[e](x,0  =  Vy(x,0  (4) 


This  expression  assumes  that  the  spatial  gradients  of  the  fluid  ambient  density  are  negligible  so  that  the 
acoustic  particle  velocity  field  is  curl-free  (i.e.,  V  x  y  =  0). 

By  combining  Eqs.  (1)  and  (2)  to  eliminate  particle  velocity,  the  acoustic  wave  equation  in  pressure 
can  be  derived.  Alternatively,  the  equations  can  be  transformed  from  first  order  to  second  order  in  the 
acoustic  variables  to  obtain  expressions  that  characterize  the  energetics  of  acoustic  fields  [6,7].  Just  as  two 
types  of  energy  density  exist  in  a  sound  field  (potential  and  kinetic),  two  types  of  energy  flux  occur.  One, 
the  "active"  intensity,  describes  the  net  transport  of  acoustic  energy  through  the  medium.  In  a  monotone 
sound  field,  it  is  proportional  to  the  spatial  gradient  of  the  wave  field  phase.  The  second  type  of  energy 
flow,  the  "reactive"  intensity,  is  proportional  to  the  spatial  gradient  of  the  pressure  autospectrum  (or  poten¬ 
tial  energy  density).  It  is  required  to  support  the  spatial  structure  of  the  field.  Its  existence  and  importance 
is  the  reason  that  ocean  acoustic  fields  cannot,  in  general,  be  modeled  as  the  superposition  of  uncorrelated 
plane  wave  arrivals. 

The  two  energy  densities  and  the  individual  components  of  the  active  and  reactive  intensities  appear 
in  the  4x4  data  cross  spectral  matrix  formed  from  vector  sensor  ("combined  receiver")  data.  The  properties 
of  the  3x3  particle  velocity  sub-cross  spectral  matrix  can  be  analyzed  in  terms  of  the  polarization  of  acous¬ 
tic  particle  motion  [7,8].  Therefore,  all  cross  spectral  matrix  terms  have  basic  physical  significance. 


III.  SOME  RESULTS  OF  AT-SEA  MEASUREMENTS 
A.  Vector  Array  Processing 

The  Taylor  series  expansion  of  the  acoustic  pressure  field  given  in  Eq.  (3)  says  that  the  measurement 
of  acoustic  pressure  and  its  higher-order  spatial  derivatives  at  a  single  point  in  space  is  equivalent  to  the 
measurement  of  acoustic  pressure  in  a  volume  about  the  measurement  point.  Therefore,  the  techniques 
used  in  beamforming  with  spatially  distributed  pressure  measurements  can  be  applied  to  single  point  meas¬ 
urements  of  pressure  and  its  spatial  derivatives.  In  particular,  since  Eq.  (1)  shows  that  acoustic  particle 
velocity  at  a  given  frequency  is  proportional  to  the  first  order  spatial  gradient  of  pressure,  then  high  resolu¬ 
tion  beamforming  techniques  can  be  used  with  simultaneous  pressure  and  particle  velocity  measurements. 
Application  of  minimum  variance  (Capon  [9])  as  well  as  conventional  beamforming  techniques  to  these 
types  of  single  point  measurements  were  discussed  in  Refs.  10  and  11  (see  also  [12,13]).  As  an  example,  a 
comparison  of  the  results  of  using  conventional  beamforming  (CBF)  and  adaptive  beamforming  (ABF)  at  a 
source  tow  frequency  of  225  Hz  with  data  from  all  properly-operating  channels  of  the  Vertical  DIFAR 
Array  is  presented  in  Fig.  3  (from  a  combination  of  Figs.  5  and  6  of  [10]).  The  elevation  angle  is  plotted 
along  the  vertical  axis,  with  negative  angles  indicating  downward-pointing  beams,  and  the  azimuth,  in 
degrees  clockwise  from  true  north,  is  given  along  the  horizontal  axis.  These  data  were  collected  when  the 
towed  source  was  located  23°  east  of  north. 
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Most  of  the  dark-colored  areas  in  the  figure  are  due  to  the  bottom-bounce  arrival,  predicted  from 
numerical  modeling  to  be  the  most  predominant  arrival  at  this  range.  Vertically-aliased  replicas  of  this 
arrival  also  occur  because  the  array  element  spacing  of  15  m  represents  more  than  two  acoustic 
wavelengths  at  225  Hz.  The  horizontal  resolution  of  the  ABF  method  is  superior  to  that  using  conventional 
DIFAR  processing  in  this  case  since  the  signal-to-noise  ratio  (SNR)  was  large. 

To  combine  the  data  from  DIFAR  elements  spaced  over  an  interval  of  depth  in  the  water  column,  a 
number  of  approaches  can  be  taken.  For  the  results  presented  in  Fig.  3,  a  conventional  delay-and-sum 
beamforming  in  the  vertical  direction  first  was  performed.  The  delay-and-summing  was  done  with  the  time 
series  from  sensors  of  the  same  type;  that  is,  all  the  pressure  time  series  were  delay-and-summed,  then  the 
vertical  geophone  data  were  delay-and-summed,  and  so  on.  This  initial  step  created,  in  effect,  the  four- 
component  time  series  for  an  equivalent  single  DIFAR  element  with  extended  vertical  aperture  correspond¬ 
ing  to  each  vertical  look  direction.  Then,  conventional  Bartlett  processing  (upper  panel)  and  minimum 
variance  Capon  processing  (lower  panel)  were  applied  to  the  four-component,  delayed-and-summed  time 
series. 

An  attempt  to  perform  global  Capon  processing,  i.e.,  estimating  and  inverting  the  full  64  x  64  cross 
spectral  matrix,  failed  due  to  the  problem  of  correlated  arrivals  from  different  vertical  directions.  Corre¬ 
lated  arrivals  result  in  a  non-negligible  reactive  energy  flux  which  is  the  cause  of  the  problem.  Real  ocean 
acoustic  wave  fields  in  most  cases  cannot  be  modeled  as  the  superposition  of  uncorrelated  plane  wave 
arrivals. 

As  discussed  in  Ref.  11,  the  distance,  Ay,  that  the  pressure  field  can  be  extrapolated  from  the  meas¬ 
urement  point  with  a  given  error  and  with  an  expansion  to  a  given  order  is  dependent  upon  the  degree  of 
spatial  variability  of  the  pressure  field.  The  maximum  spatial  variability  in  any  direction  at  a  given  fre¬ 
quency  cannot  exceed  that  determined  by  the  acoustic  wavelength,  at  least  for  non-evanescent  acoustic 
fields.  Otherwise,  the  velocity  of  the  energy  flow  required  to  support  the  spatial  structure  would  exceed  the 
speed  of  sound  in  the  fluid,  which  is  physically  impossible.  Therefore  the  effective  spatial  aperture  of  a 
single  point  array  is  related  to  the  acoustic  wavelength,  rather  than  being  determined  by  a  fixed  inter¬ 
element  spacing  as  with  conventional  arrays.  Thus,  single  point  arrays  are  frequency-adaptive;  their  effec¬ 
tive  aperture  decreases  with  increasing  frequency.  The  result  of  this  property  is  that  the  plane  wave 
response  (or  beampattem)  for  single  point  arrays  is  independent  of  frequency.  This  result  also  implies  that 
grating  lobes  do  not  exist;  spatial  aliasing  cannot  occur  since  a  sampling  in  space  is  not  being  performed. 

The  phenomenon  of  superdirectivity  is  directly  related  to  the  Taylor  series  expansion  of  the  pressure 
field.  That  is,  superdirectivity  arises  for  a  spatially  distributed  hydrophone  array  when  the  directivity  index 
is  maximized  as  a  function  of  the  element  weights,  and  the  interelement  spacing  becomes  smaller  than  half 
the  acoustic  wavelength  [14].  One  can  show  that,  as  the  ratio  of  the  interelement  spacing  to  the  acoustic 
wavelength  approaches  zero,  the  weights  for  a  "linear  point  array"  approach  the  finite  difference  approxi¬ 
mations  to  the  spatial  derivatives  of  pressure  given  in  Eq.  (3).  The  instability  that  results  when  the  weights 
become  large  and  of  opposite  sign  can  be  avoided  by  the  use  of  alternative  transduction  methods  suggested 
by  the  physical  interpretation  of  the  spatial  derivatives  of  pressure,  e.g.,  when  the  measurement  of  a  com¬ 
ponent  of  acoustic  particle  velocity  replaces  the  measurement  of  pressure  at  two  closely-spaced  points. 

Note  that  high  resolution  beamforming  techniques  have  been  applied  to  single  point  measurements  of 
other  wave  fields,  e.g.,  the  estimation  of  ocean  surface  gravity  wave  directional  spectra  from  "pitch-and- 
roll"  buoys  [15]. 

B.  Bioacoustics 

In  1996,  the  Marine  Physical  Lab  conducted  an  experiment  just  south  of  San  Miguel  Island  bordering 
the  Santa  Barbara  Channel  as  part  of  its  Marine  Mammal  Vocalization  program.  The  purpose  of  this 
ONR-sponsored  program  was  to  study  the  natural  calling  behavior  of  large  baleen  whales  (blue,  fin,  and 
humpback  whales)  and  to  investigate  how  their  calls  could  be  used  as  loud,  low-frequency  sources  of 
opportunity  to  determine  properties  of  the  ocean  environment.  During  one  period,  recordings  of  vocalizing 
blue  whales  were  made  simultaneously  by  a  vertical  line  array  of  48  hydrophones  and  a  freely  drifting, 
standard  Navy  SSQ  53B  sonobuoy.  Fig.  4  shows  a  5-min  spectrogram  from  the  omnidirectional  com¬ 
ponent  of  the  DIFAR  sonobuoy.  Two  characteristic  Type  A  calls  (the  second  call  after  the  start  of  the  5- 
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min  period  and  the  one  at  the  4-min  mark)  and  6  Type  B  calls  with  high  SNR  (17-Hz  FM  sweeps  and  asso¬ 
ciated  higher  harmonics)  are  present.  The  issue  is  to  determine  the  number  and  locations  of  animals  contri¬ 
buting  to  this  sequence  of  calls.  If  more  than  one  individual  is  vocalizing,  an  interesting  question  is  how,  if 
at  all,  calling  behavior  is  altered  by  the  presence  of  other  calling  animal(s). 

The  5  min  of  data  from  the  sonobuoy  were  divided  into  consecutive  5 -sec  blocks  and  the  average 
active  intensity  per  frequency  within  each  block  was  calculated.  (Sonobuoy  data  permit  only  the  horizontal 
component  of  active  intensity  to  be  determined).  The  upper  panel  in  Fig.  5  shows  a  plot  of  the  horizontal 
active  intensity  in  a  2.4-Hz-wide  bin  centered  at  51  Hz  as  it  evolved  over  the  5-min  period.  The  51-Hz  bin 
contains  the  energy  in  the  3rd  harmonic  of  the  Type  B  calls  and  little  of  the  energy  from  the  Type  A  calls. 
The  figure  was  created  by  plotting  the  sequential  average  active  intensity  vector  estimates  consecutively 
head  to  tail,  starting  at  the  origin  of  the  compass.  The  beginning  and  end  of  each  vector  is  indicated  by  a 
small  "x"  (as  are  the  10°  intervals  on  the  outer  compass  circle).  The  plot  clearly  shows  that  the  first  two  B 
calls  created  net  energy  flux  in  the  WNW  direction,  1 80°  from  the  source  to  the  ESE  (approximate  bearing 
of  120°).  The  third  B  call  generates  energy  flow  to  the  SE,  indicating  a  second  vocalizing  animal  to  the 
NW  (bearing  340°).  The  final  three  Type  B  calls  created  progressively  weaker  flux  to  the  ENE,  in  the 
same  direction  as  the  first  two  calls.  Therefore,  two  animals  appear  to  be  the  sources  of  the  Type  B  calls, 
one  to  the  NW  and  a  second,  more  "talkative"  one  to  the  ESE. 

The  corresponding  plot  for  the  91  Hz  bin,  which  contains  the  energy  from  the  Type  A  calls  but  no 
Type  B  energy,  is  shown  in  the  lower  panel  in  Fig.  5.  The  91  Hz  bin  also  contains  energy  from  a  spectral 
line  (re  Fig.  4),  radiated  by  the  R/P  FLIP,  the  platform  from  which  the  experiment  was  conducted.  Upon 
deployment  from  FLIP,  the  sonobuoy  drifted  almost  due  north,  towards  San  Miguel  Island  away  from 
FLIP.  The  91 -Hz  energy  flow  at  the  beginning  of  the  5-min  period,  and  at  other  times  when  neither  of  the 
Type  A  signals  were  present,  is  to  the  north  in  the  direction  away  from  FLIP.  The  first  Type  A  signal  inter¬ 
rupts  this  northward  flow  by  creating  flux  to  the  SE,  away  from  a  source  at  bearing  335°.  It  most  likely 
was  created  by  the  same  animal  that  generated  the  third  Type  B  call  in  the  previous  figure.  The  second  A 
call  clearly  is  coming  from  the  same  direction  as  the  more  "talkative"  animal  from  the  previous  figure. 
Therefore,  two  animals  most  likely  are  contributing  to  the  calling  sequence  in  Fig.  4. 

Range  and  depth  estimates  of  these  two  calling  animals,  determined  by  matched-field  processing 
with  the  vertical  hydrophone  array  data,  are  discussed  on  pp.  1291-1292  and  presented  in  Fig.  8  of  Ref.  16. 

These  DIFAR  sonobuoy  data  demonstrate  one  great  benefit  of  such  systems;  they  allow  measure¬ 
ment  of  the  directionality  of  very  low  frequency  signals.  DIFAR  sonobuoy s  along  with  DIFAR  demulti¬ 
plexing  systems  and  display  units,  are  beginning  to  be  routinely  used  in  marine  mammal  acoustics  studies 
[17]. 

C.  Seismology 

Earthquake-generated  T  phase  arrivals  have  been  recorded  in  almost  every  Swallow  float  deploy¬ 
ment.  Association  of  these  arrivals  with  specific  events  was  accomplished  by  calculating  the  travel  times 
to  epicentral  locations  of  major  events  posted  in  earthquake  bulletins  and  comparing  the  listed  origin  times 
with  those  predicted  from  the  Swallow  float  data.  Agreement  of  the  two  origin  times  typically  was  within  a 
few  secs.  An  example  of  a  Swallow  float  recording  of  an  earthquake  is  presented  in  Ref.  6.  During  the 
IONEX  92  experiment,  over  50  aftershocks  from  a  subduction  zone  event  off  Crete  were  recorded  [18].  In 
addition,  several  aftershocks  of  the  1992  Lander’s  earthquake  in  southern  California  were  recorded  during 
a  vertical  DIFAR  array  deployment  in  the  southern  California  Bight  region. 

Another  source  of  land-based  seismic  signals  that  can  couple  into  the  water  columm  is  shown  in  Fig. 
6.  Offshore  underwater  sounds  from  these  tracked  vehicles  operating  on  land  were  recorded  during  MPL’s 
Adaptive  Beach  Monitoring  program.  As  an  example,  the  adaptive  plane  wave  beamforming  results  over 
the  30-70  Hz  frequency  band  as  a  function  of  time  during  the  transit  of  four  such  tracked  vehicles  down- 
coast  along  the  beach  is  shown  in  Fig.  7  [19].  These  results  were  obtained  from  data  collected  by  a  2D 
hydrophone  array  on  the  ocean  bottom  1.5  km  offshore.  The  beach,  to  the  NE  from  the  array  deployment 
site,  runs  NW-SE  at  this  site.  The  directionality  indicated  by  the  acoustic  results  agrees  well  with  visual 
observations  during  the  experiment. 
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Data  collected  over  this  same  time  period  by  the  OSZM’s  bottom-mounted  geophone/hydrophone 
package  (Sect.  IC)  were  processed  using  acoustic  vector  intensity  techniques,  similar  to  the  approach  used 
for  the  blue  whale  calls  in  Fig.  5.  The  horizontal  projection  of  the  active  intensity  component  integrated 
across  the  bands  from  34  to  56  Hz  ("x'j  and  from  64  to  70  Hz  (triangles)  as  a  function  of  time  over  the 
same  time  period  as  Fig.  7  is  shown  in  Fig.  8.  For  the  34-56  Hz  band,  the  background  energy  flow  before 
the  start  of  the  vehicle  run  is  from  the  southwest  to  the  northeast.  The  start  of  the  run  resulted  in  flow  to  the 
south,  which  evolves  into  a  flow  to  the  southwest  as  the  vehicles  traveled  down  the  beach.  The  majority  of 
the  flow  in  this  lower  band  occurs  in  a  more  southerly  direction  than  that  in  the  higher  frequency  band. 
The  latter  is  dominated  by  flow  1 80°  away  from  the  bearing  to  the  coastline,  corresponding  to  the  location 
where  the  vehicles  passed  through  CPA.  Whereas  the  64-to-70-Hz  flow  continues  to  evolve  towards  a  wes¬ 
terly  direction  as  the  vehicles  proceeded  downcoast,  it  does  not  do  so  in  the  lower  frequency  band. 
Although  changes  in  land-to-water  coupling  characteristics  may  play  a  role,  these  frequency  band  differ¬ 
ences  are  probably  due  to  the  fact  that  the  frequency  content  of  the  vehicle  signals  shifted  to  higher  fre¬ 
quencies  as  they  increased  speed  upon  reaching  the  open  part  of  the  beach.  Results  of  adaptive  beamform¬ 
ing  with  the  hydrophone  array  data  over  these  two  frequency  bands  separately  are  consistent  with  those  in 
Fig.  8. 

As  the  previous  results  illustrate,  the  geophone/hydrophone  package  clearly  provides  useful  informa¬ 
tion  on  the  time-evolving  directionality  of  the  land  vehicle  sounds.  It  does  so  without  the  need  to  know  the 
phase  velocity  of  the  arriving  energy  since  the  individual  sensor  channels  are  located  at  a  single  point  in 
space.  To  obtain  the  broadband  adaptive  plane  wave  beamforming  results  in  Fig.  7,  the  phase  velocity  as  a 
function  of  frequency  first  had  to  be  measured.  The  phase  velocity  in  this  case  is  a  strong  function  of  fre¬ 
quency  since  the  received  energy  is  carried  by  the  lowest  order  normal  mode  which  is  near  cutoff  at  these 
low  frequencies  in  these  very  shallow  waters  [19]. 

D.  Physical  Acoustics 

Dislocations  are  places  in  the  acoustic  field  where  the  amplitude  goes  to  zero  and  the  total  field  phase 
becomes  undetermined  [20].  Although  they  are  curious  features  of  the  field  in  their  own  right,  they  may 
also  be  sensitive  indicators  of  fluctuations  in  the  environment  [21].  Three  candidate  dislocations  in  the 
pressure  field  are  indicated  by  vertical  dashed  lines  in  Fig.  9.  The  solid  curve  is  the  acoustic  pressure  spec¬ 
trum  and  the  dotted  curve  is  the  vertical  component  particle  velocity  spectrum.  These  data  were  collected 
by  a  Swallow  float  deployed  at  nearly  1600  m  depth  during  the  1990  NATIVE  1  experiment  near  the  Blake 
Plateau  in  the  northwest  Atlantic  Ocean.  Those  in  Fig.  9  were  obtained  during  a  2.7-hour  period  (record 
number  along  the  abscissa  is  a  measure  of  time)  as  the  source,  transmitting  a  7.0  Hz  tone,  was  towed  out  to 
distance  of  nearly  19  km. 

At  a  dislocation  in  the  pressure  field,  not  all  the  correponding  particle  velocity  components  will  be 
zero.  That  is,  although  the  pressure  field  amplitude  may  be  zero,  its  spatial  gradient,  upon  which  the  parti¬ 
cle  velocity  also  depends,  will  not  be  unless  the  field  amplitude  is  negligible  throughout  a  finite  volume 
(i.e.,  a  shadow  zone  rather  than  dislocation).  Fig.  9  shows  that  at  many  of  the  pressure  field  minima,  the 
vertical  geophone  component  spectrum  actually  attains  a  local  maximum. 

The  vector  properties  of  acoustic  intensity  can  be  used  to  determine  what  happens  in  the  neighbor¬ 
hood  of  a  dislocation.  The  directions  of  the  active  and  reactive  intensities  change  quickly  in  such  regions. 
In  particular,  the  radial  component  of  reactive  intensity  changes  sign  in  passing  by  a  dislocation.  Also,  the 
magnitude  of  the  reactive  intensity  becomes  large  since  greater  energy  flow  is  required  to  support  the 
existence  of  the  acoustic  hole.  In  addition,  acoustic  particle  motion  near  a  dislocation  is  strongly  non- 
rectilinear  since  the  active  and  reactive  flow  directions  become  misaligned  (discussed  below). 

Fig.  10  shows  plots  of  the  projections  of  the  intensity  vectors  in  the  radial/vertical  plane  (radial  in  the 
forward  direction  is  the  direction  180°  away  from  the  source)  over  a  15-min  period  (20  records)  about  the 
first  pressure  field  minimum,  occurring  at  record  996  in  Fig.  9.  The  active  intensity  (upper  panel)  is  dom¬ 
inated  mostly  by  the  horizontal  component,  representing  the  direction  of  net  flux  of  energy  down  the 
waveguide.  However,  near  the  candidate  dislocation  at  record  996,  the  radial  active  intensity  component 
decreases  significantly  so  that  net  flux  occurs  predominantly  in  the  vertical  direction.  Also,  in  passing  by 
this  location,  the  vertical  active  intensity  component  flips  from  the  downward  to  the  upward  direction.  The 
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lower  panel  of  Fig.  10  shows  the  corresponding  reactive  intensity.  It  is  mostly  vertical  since  the  spatial 
structure  in  ocean  waveguides  typically  is  significantly  greater  in  that  direction  than  in  the  radial.  The 
reactive  intensity  can  vary  wildly  in  places  where  the  spatial  variations  in  the  pressure  field  are  very  small, 
as  near  the  beginning  of  the  20-record  period.  Of  interest  here,  however,  is  how  both  the  radial  and  vertical 
reactive  intensity  components  flip  in  direction  at  record  996.  The  behavior  of  the  two  types  of  energy  flow 
seen  in  Fig.  10  is  consistent  with  that  expected  near  a  dislocation. 

The  polarization  of  the  particle  motion  in  a  way  combines  the  information  contained  in  the  active  and 
reactive  intensity  components.  As  discussed  in  Ref.  8  (see  pp.  218-219  of  [7]  for  derivation),  the  following 
relationship  exists  between  the  two  types  of  vector  acoustic  intensity  and  the  imaginary  components  of  the 
3x3  particle  velocity  data  cross  spectral  matrix: 


Im 


j[S  Y](f) 


Sp(f) 


(5) 


This  relationship  can  be  used  to  define  the  following  two  nondimensional  measures  of  the  "degree  of  cir¬ 
cularity"  of  the  particle  motion: 

£n>(Ox&v</) 

0!  =  2  — - — -  (6a) 

Sp(f)USv(f ) 


and 


Imj|SvK/)j 

02  =  2 -  (6b) 

Purely  circular  motion  results  in  a  value  for  0  of  ±1,  with  the  sign  distinguishing  between  prograde  and 
retrograde  directions,  and  purely  rectilinear  motion  has  a  value  of  zero.  These  two  measures  of  the  degree 
of  circularity  are  plotted  in  Fig.  11  over  the  same  time  period  as  in  Fig.  9;  Eq.  (6a)  is  plotted  with  a  dotted 
curve  and  Eq.  (6b)  with  a  solid  curve.  Both  measures  agree  extremely  well  over  the  full  duration  of  the 
source  tow.  As  indicated  by  the  vertical  dashed  line  at  record  996  in  Fig.  11,  the  particle  motion  reverses 
from  being  nearly  perfectly  circular  in  one  direction  to  circular  motion  in  the  opposite  direction.  Again, 
this  result  is  consistent  with  that  expected  in  the  neighborhood  of  a  dislocation. 

These  and  additional  results  on  possible  dislocations  in  Swallow  float  data  were  presented  at  the  most 
recent  Acoustical  Society  of  America  meeting  [22] . 


IV.  ’TENSOR  ”  SENSOR 

To  increase  the  directivity  index  (DI;  defined  as  the  array  gain  of  a  single  plane  wave  in  an  isotropic 
noise  field)  for  a  conventional,  spatially  distributed  pressure  sensor  array,  more  sensors  at  the  appropriate 
spacing  simply  are  added  (assuming  spatial  coherence  of  the  signal  across  the  increased  array  aperture). 
For  an  array  that  theoretically  exists  at  a  single  point  in  space,  an  increase  in  DI  is  obtained  by  making 
measurements  of  increasingly  higher  order  spatial  derivatives  of  the  acoustic  pressure  field.  Pressure  is  a 
scalar  quantity,  equivalent  to  a  tensor  of  rank  zero.  The  gradient  of  pressure  is  a  vector  quantity,  a  first 
rank  tensor,  and  the  gradient  of  a  vector  is  a  2nd  rank  tensor. 

A  physical  interpretation  of  the  term  at  second  order  in  the  Taylor  series  expansion  in  Eq.  (3)  is 
given  by  Eq.  (4).  That  is,  for  an  acoustic  field  in  an  otherwise  stationary  fluid,  the  spatial  gradient  of  the 
acoustic  particle  velocity  is  equal  to  the  acoustic  rate  of  strain.  The  concept  of  strain  rate  usually  appears 
in  fluid  mechanics  with  reference  to  viscous  dissipation.  However,  a  sound  field  causes  all  kinds  of  spatial 
distortions  of  the  fluid  medium  from  its  ambient  state  so  that  strain  rate  also  is  a  useful  concept  in  acoustics. 

Results  of  investigating  the  benefits  of  measuring  the  quantities  out  to  second  order  in  the  Taylor 
series  expansion  were  presented  in  [13]  and  summarized  in  [11].  The  investigation  involved  both  the 
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beamforming  properties  of  a  "tensor  sensor"  and  the  vector  and  tensor  properties  of  the  terms  in  the  data 
cross  spectral  matrix  in  the  context  of  the  energetics  of  acoustic  fields.  Those  results  are  repeated  here. 

A.  Tensor  Sensor  Beamforming  Properties 

The  beamformer  output  time  series  for  a  tensor  sensor  (originally  called  the  "double  DIFAR"  in  [13]) 
can  be  written  as 

3  i  3  3 

d(2\x ,  t)  =  a0p(x,  t)  +  '£jajZjvj(x,  t )cos(py- )  +  —'Z'Lbij  Wij £ij (* > f  )cos(p, )cos(Py )  (7) 

/-!  2  >=1/=1 

where  8,-y  is  the  ij- th  term  in  the  2nd  rank  strain  rate  tensor,  [e],  the  Zy-  ’s  and  s  are  conversion  factors 
(not  necessarily  purely  real),  the  cos((37  )’s  are  the  direction  cosines,  and  the  cij  ’ s  are  s  are  arbitrary 
weights.  Different  approaches  to  beamforming  can  be  defined  by  determining  the  arbitrary  weights  in  vari¬ 
ous  ways.  The  conventional  beamforming  approach  sets  the  weights  to  fixed  and  "equal"  (properly  nor¬ 
malized)  values.  The  standard  DIFAR  beamformer  with  a  cardioid  beampattern  (re  the  dashed  curve  in 
Fig.  12b)  is  an  example  of  this  approach.  It  provides  a  DI  of  4.8  dB,  has  a  main  lobe  width  (defined  by  the 
3-dB  down  points)  of  131°,  has  no  side  lobes  and  just  one  null. 

Alternatively,  the  weights  can  be  determined  under  various  optimization  criteria.  For  example,  they 
can  be  determined  so  that  the  DI  is  maximized.  A  DIFAR  beamformer  at  first  order  has  a  maximum  DI  of 
6  dB,  a  main  lobe  width  of  105°,  a  single  side  lobe  of  level  -6  dB  (180°  from  direction  of  maximum  sensi¬ 
tivity),  and  2  nulls.  The  plane  wave  response  of  this  beamformer  is  plotted  as  a  dotted  curve  in  Fig.  12a. 
In  comparison,  a  tensor  sensor  has  a  maximum  DI  of  9.5  dB,  a  main  lobe  width  of  65  deg,  side  lobe  level 
of  -9.5  dB,  and  4  nulls.  The  resulting  plane  wave  response  having  maximum  DI  is  shown  as  the  solid  curve 
in  Fig.  12a.  The  solid  curve  in  Fig.  12b  is  the  tensor  sensor  beampattern  that  results  from  applying  the  con¬ 
straint  that  a  null  exist  180°  from  the  direction  of  maximum  sensitivity.  It  has  a  maximum  DI  of  9.0  dB,  a 
main  lobe  width  of  70°,  side  lobes  that  are  13.0  dB  down  in  level,  and  3  nulls. 

A  second  optimization  approach  is  to  determine  the  weights  using  a  data-adaptive  approach  of 
minimizing  the  beamformer  output  variance  under  various  constraints  (e.g.,  Capon  processing).  Fig.  13 
shows  the  results  of  a  tensor  sensor  simulation  in  which  the  acoustic  field  is  created  by  four  sources  having 
the  following  combination  of  directions  of  arrival  (DOA)  and  signal-to-noise  ratios  (SNR);  (DOA:  0°, 
SNR:  12  dB),  (DOA:  45°,  SNR:  9  dB),  (DOA:  180°,  SNR:  6  dB),  and  (DOA:  270°,  SNR:  3  dB).  The  solid 
curve  is  the  beamformer  output  using  the  conventional  beamforming  approach  and  the  dotted  curve  shows 
the  Capon  processor  output.  The  presence  of  4  sources  is  clearly  indicated  in  the  Capon  results  since  this 
processor  has  some  control  over  the  direction  in  which  it  steers  its  nulls. 


B.  Energetics  of  the  Cross  Spectral  Matrix 


By  using  tensor  and  vector  algebra  and  the  fundamental  physical  laws  expressed  in  Eqs.  (1)  and  (2), 
relationships  between  the  various  quantities  in  the  tensor  sensor  10x10  data  cross  spectral  matrix  can  be 
derived.  This  approach  is  exactly  the  same  used  to  derive  the  vector  properties  of  the  two  components  of 
acoustic  intensity.  One  property  of  the  acoustic  pressure/strain  rate  tensor  is 


P 


[e]  =  V(pv)  + 


p0(v)(v) 


(8) 


Also,  it  can  be  shown  using  the  strain  energy  function  that  the  trace  of  p  [8]  (i.e.,  the  sum  of  the  diagonal 
elements)  yields  the  conservation  of  acoustic  energy  equation.  For  stationary,  stochastic  fields  in  the  fre¬ 
quency  domain: 


Re 


5 

P  [e] 


(f) 


(9) 


Therefore,  the  cross  spectrum  between  the  acoustic  pressure  and  the  acoustic  rate  of  strain  provides  infor¬ 
mation  on  the  local  spatial  heterogeneity  of  the  active  acoustic  intensity.  Likewise,  a  property  of  the 
velocity/strain  rate  second  rank  tensor  is 
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p0v*  [e]  =  V{K.E.D.|  (10) 

where  "K.E.D."  signifies  the  kinetic  energy  density. 

The  existence  of  relationships  between  some  of  the  terms  in  the  tensor  sensor  data  cross  spectral 
matrix  implies  that  they  do  not  provide  entirely  independent  information  on  the  acoustic  field.  However, 
these  relationships  provide  a  way  of  quality  checking  the  data  from  such  systems.  We  hope  to  have  the 
opportunity  to  build  a  prototype  tensor  sensor  in  the  near  future. 

An  interesting  question  is  -  do  high  quality  ocean  acoustic  strain  rate  sensor  systems  already  exist  in 
nature  [23]? 

V.  CONCLUSIONS 

A  sensor,  or  array  of  sensors,  that  simultaneously  measures  acoustic  pressure  and  the  three  com¬ 
ponents  of  acoustic  particle  velocity  adds  no  new  information  to  what  can  be  obtained  by  a  set  of  spatially 
distributed  pressure  measurements.  However,  vector  and  "tensor"  sensors  have  some  distinct  advantages; 
they  are  more  compact  in  size,  can  be  easier  to  deploy,  and  do  not  suffer  from  those  errors  introduced  by 
element  location  uncertainty  or  uncertainty  in  the  phase  speed  of  propagation  across  the  array  aperture. 
They  may  be  the  only  practical  way  of  obtaining  directional  properties  of  the  underwater  acoustic  field  at 
infrasonic  frequencies.  In  addition,  a  single  point  vector  or  tensor  sensor  has  a  beam  pattern  that  is 
independent  of  frequency,  making  it  convenient  for  broadband  studies.  Finally,  the  terms  in  these  sensors’ 
data  cross  spectral  matrices  are  of  basic  physical  significance  so  that  results  often  are  more  easily  inter¬ 
preted  and  insights  more  readily  obtained. 

The  greatest  challenge  in  this  field  is  making  high  quality  ocean  acoustic  measurements  with  these 
types  of  sensor  systems.  Motion-induced  self  noise  is  a  particularly  vexing  problem. 

Finally,  the  beautiful  and  complicated  interference  patterns  displayed  by  real  ocean  acoustic  fields 
rely  upon  the  reactive  energy  flow  for  their  existence.  Such  patterns,  which  contain  very  useful  informa¬ 
tion  on  source  location  and  environmental  properties,  are  completely  absent  in  fields  that  can  be  modeled 
as  the  superposition  of  uncorrelated  plane  wave  arrivals. 
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FIGURE  CAPTIONS 

Figure  1.  Schematic  of  an  MPL  Swallow  float. 

Figure  2.  Schematic  of  2  elements  of  MPL’s  16-element,  225-m-aperture,  Vertical  DIFAR  Array. 

Figure  3.  Comparison  of  conventional  and  adaptive  (Capon)  processing  at  225  Hz  from  at-sea  data  col¬ 
lected  by  the  Vertical  DIFAR  Array  in  1991. 

Figure  4.  5-min  spectrogram  from  the  omnidirectional  component  of  a  53B  DIFAR  sonobuoy  deployed 
during  MPL’s  1996  Marine  Mammal  Vocalization  experiment. 

Figure  5.  Upper  panel  is  the  horizontal  active  intensity  at  51  Hz  over  the  same  5-min  time  period  and  from 
the  same  DIFAR  sonobuoy  as  in  Fig.  4.  Lower  panel  shows  the  corresponding  result  at  91  Hz. 

Figure  6.  Amphibious  tracked  vehicle  ("amtrack")  traveling  down  the  beach  during  MPL’s  1996  Adaptive 
Beach  Monitoring  experiment  (ABM  96). 

Figure  7.  Adaptive  beamforming  over  the  30-70  Hz  frequency  band  during  a  10-min  period  in  ABM  96 
using  data  from  the  64-element  hydrophone  line  array  deployed  quasi-parallel  to  the  coast  1.5  km  offshore 
in  12-m-deep  water. 

Figure  8.  Horizontal  active  intensity  from  the  OSZM  data  integrated  over  the  bands  from  34  to  56  Hz 
("x")  and  from  64  to  70  Hz  (triangles)  as  a  function  of  time  during  the  10-min  period  shown  in  Fig.  7.  Also 
plotted  as  a  smooth  curve  is  the  north-south  geophone  velocity  spectral  density  against  that  from  the  east- 
west  component  as  a  function  of  time,  both  integrated  across  the  34  to  56  Hz  band. 

Figure  9.  Hydrophone  (solid  curve)  and  vertical  geophone  (dotted  curve)  spectral  levels  at  7  Hz  as  a  func¬ 
tion  of  time  (20  records  equals  15  min)  recorded  by  a  Swallow  float  at  nearly  1600  m  depth  during  source 
tow  in  the  NATIVE  1  experiment. 

Figure  10.  Time  series  of  the  active  (upper  panel)  and  reactive  (lower  panel)  intensities  at  7  Hz  over  a 
1 5-min  time  period  about  the  first  pressure  spectral  null  (marked  by  the  lefthand-most  vertical  dashed  line) 
in  Fig.  9.  Plotted  is  the  projection  of  the  intensity  vectors  in  the  radial/vertical  plane.  The  compasses  in  the 
upper  right  corners  indicate  Up,  Down,  Forward,  and  Back. 

Figure  11.  Degree  of  circularity  of  particle  motion  in  the  radial/vertical  plane  at  7  Hz  over  the  same  time 
period  as  Fig.  9. 

Figure  12.  2D  plane  wave  responses  for  various  component  weightings  for  first  order  DIFAR  sensors  and 
second  order  tensor  sensors.  In  Fig.  12a,  the  dotted  curve  is  for  the  DIFAR  sensor  with  maximum  DI  of 
6.0  dB  and  the  solid  curve  is  for  the  tensor  sensor  with  maximum  DI  of  9.5  dB.  In  Fig.  12b,  the  2D  plane 
wave  response  for  the  standard  DIFAR  with  equal  component  weighting  (DI  =  4.8  dB)  is  plotted  with  a 
dashed  curve,  and  the  tensor  sensor  with  maximum  DI  under  the  constraint  that  a  null  exist  in  the  180° 
direction  (DI  =  9.0  dB)  is  plotted  with  a  solid  curve. 

Figure  13.  Simulation  of  conventional  Bartlett  (solid  curve)  and  Capon  (dotted  curve)  beamformer  outputs 
of  a  tensor  sensor  for  a  field  composed  of  four  sources  having  directions  of  arrival  (DOA)  and  signal-to- 
noise  ratios  of  (DOA:  0°,  SNR:  12  dB),  (DOA:  45°,  SNR:  9  dB),  (DOA:  180°,  SNR:  6  dB),  and  (DOA: 
270°,  SNR:  3  dB). 
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Abstract:  Vector-based  directional  sensors  present  an  effective  solution  to  an  important  subset  of 
measurement  problems  in  underwater  acoustics.  There  is  an  inherent  directionality  in  the  vector 
field  that  cannot  be  obtained  from  a  pressure  measurement  at  a  single  point.  Inertial  sensors  in 
neutrally  buoyant  packages  respond  directly  to  acoustic  particle  motion  without  the  dynamic- 
range  penalty  of  gradient  sensors;  however,  the  conflicting  requirements  of  supporting  the  sensor 
without  interfering  with  its  response  to  the  acoustic  field  present  interesting  design  challenges. 
Furthermore,  their  susceptibility  to  non-acoustic  excitation  makes  high  performance  difficult  to 
achieve  in  some  circumstances.  The  gradient  sensor  relies  on  a  finite-difference  approximation  to 
the  pressure  gradient.  The  elements  of  the  gradient  sensor  do  not  need  to  move  with  the  velocity 
field  and  the  elements  can  be  made  much  smaller  allowing  higher  frequency  operation.  However, 
the  output  of  a  gradient  sensor  depends  on  subtraction  of  two  nearly  equal  quantities,  which 
severely  limits  the  dynamic  range.  Success  in  the  application  of  vector  sensors  depends  on  an 
understanding  of  their  physics  and  their  limitations. 


I.  INTRODUCTION 

There  are  at  least  two  fundamentally  different  approaches  to  achieving  directional 
discrimination  in  acoustic  sensors.  One  approach  is  to  use  arrays  of  single  elements  with 
their  outputs  combined  to  form  beams.  Such  arrays  are  typically  several  wavelengths  in 
dimension.  Another  approach  is  to  measure  vector  properties  of  the  acoustic  field.  A  sensor 
that  responds  to  either  the  local  pressure  gradient  or  the  local  acoustic  particle  velocity  has  an 
intrinsic  directionality.  Large-aperture  arrays  easily  outperform  individual  vector  sensors  in 
degree  of  directionality  but  incur  an  obvious  size  penalty.  While  a  great  body  of  literature 
exists  on  the  design  and  performance  of  large-aperture  arrays,  vector  sensors  are  less  well 
documented.  Military  and  civilian  applications  consume  millions  of  gradient  or  velocity 
sensors  of  various  types,  so  there  is  a  substantial  body  of  corporate  knowledge,  though. 

Terminology  is  confusing  and  inconsistent  in  the  literature.  In  this  summary  paper,  a 
sensor  that  responds  directly  to  incident  pressure  will  be  called  a  scalar  or  pressure  sensor.  A 
sensor  that  responds  directly  to  some  vector  measure  of  the  acoustic  field  will  be  called  a 
vector  sensor.  If  a  vector  sensor  measures  by  finite-difference  approximation  of  the  field 
gradient  it  will  be  called  a  gradient  sensor;  if  a  vector  sensor  responds  directly  to  motion  of 
the  sensor  body  with  an  output  directly  related  to  acoustic  displacement,  velocity,  or 
acceleration,  it  will  be  called  an  inertial  sensor.  Sometimes  the  terms  velocity  sensor  and 
pressure-gradient  sensor  are  used  interchangeably  in  the  literature,  but  this  obscures  the 
physics.  The  manner  in  which  gradient  and  inertial  sensors  respond  to  the  field  is 
substantially  different. 
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A  short,  superficial  history  of  the  development  of  vector  sensors  illustrates  some  of  the 
important  design  issues.  The  ribbon  microphone  was  one  of  the  first  widely  known  vector 
sensors.  In  1931,  Olson1  describes  the  relationship  between  pressure  gradient  and  velocity 
and  outlines  the  near-field  correction  for  the  amplitude  of  the  inferred  pressure  field.  Two 
years  later,  Wolff  and  Massa2  broach  the  possibility  of  measuring  acoustic  intensity  from  the 
product  of  velocity  and  pressure.  A  neutrally  buoyant  inertial  sensor  for  underwater 
acoustics  was  built  and  tested  in  1941  by  Kendall.3  The  concepts  embodied  by  Kendall’s 
design  are  common  to  most  inertial  vector  sensors  used  underwater  since  that  time. 

Much  of  the  research  directed  toward  vector  sensors  was  funded  by  the  US  Navy  toward 
development  of  air-deployed  underwater  sensors  (as  sonobuoys).  For  expendable  sensors,  the 
small  size  of  the  vector  sensor  outweighs  the  lower  directionality  in  comparison  to  large- 
aperture  arrays.  Bauer  and  DiMattia,4  supported  by  the  Naval  Air  Development  Center, 
described  the  design  of  a  two-component  inertial  vector  sensor  using  moving-coil  transducers 
in  1966.  The  directional  sonobuoy  (DIFAR5)  was  approved  for  operational  service  for  US 
Navy  Maritime  Patrol  Aircraft  by  1969.  Since  that  time,  US  manufacturers  alone  have 
produced  more  than  four  million  directional  sonobuoys.  The  DIFAR  sonobuoy  has  become 
the  principal  tool  in  maritime  patrol  activity  with  a  high  performance-to-price  ratio. 

One  of  the  critical  limitations  of  the  vector  sensor  reached  the  open  literature  in  1977 
with  Keller’s  summary6  of  measured  flow-induced  noise  on  gradient  and  inertial  sensors. 

The  sensitivity  to  non-acoustic  disturbance  is  much  greater  than  for  a  simple  pressure  sensor. 
These  measurements  showed  20-dB  higher  flow  noise  on  a  vector  sensor  compared  to  a 
pressure  sensor  at  100  Hz  for  a  flow  speed  of  0.5  knots  (increasing  to  70  dB  higher  at  10  Hz). 
Much  of  the  subsequent  work  on  vector  sensors  has  been  aimed  at  controlling  the  response  to 
flow-induced  noise. 

In  1980,  the  vertical-line-array  DIFAR  (VLAD)  received  approval  for  service.  This 
sonobuoy  combined  a  large-aperture  vertical  array  with  a  two-axis  vector  sensor.  This 
sonobuoy  is  particularly  noteworthy  because  the  design  philosophy  was  not  to  steer  a  narrow 
beam  toward  the  intended  target.  Instead,  it  was  designed  to  reject  prevailing  ambient  noise. 
As  a  result,  in  proper  application,  it  achieves  an  array  gain  far  higher  than  its  directivity  index. 

Perhaps  the  ultimate  in  low-frequency  performance  of  an  underwater  vector  sensor  was 
achieved  by  the  Marine  Physical  Laboratory’s  acoustic  Swallow  Float  system.  Also  notable 
is  the  in-depth  analysis  of  limits  to  performance  published  in  1991  by  D’Spain,  et  al.7  The 
quasi-free  sensor  body  was  incorporated  in  the  neutrally  buoyant  infrasonic  DIFAR 
sonobuoy  (NBIDS)  developed  in  the  1990’s. 

Several  themes  are  evident  throughout  this  history.  The  vector  sensor  was  and  is  highly 
successful  as  an  air-expendable  instrument.  The  sonobuoy  industry  markets  a  two-axis 
vector  sensor  complete  with  RF  telemetry,  a  sophisticated  sea-surface  motion  isolation 
system,  electronic  magnetic  compass,  and  military  qualified  air-launch-hardened  package  for 
less  than  $500.  Perhaps  the  most  serious  shortcoming  of  the  vector  sensor  is  its  sensitivity  to 
motion  (in  the  case  of  the  inertial  sensor)  or  to  non-acoustic  pressure  fluctuations  (in  the  case 
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of  the  gradient  sensor).  Certainly,  the  success  of  a  technology  in  one  arena  does  not 
guarantee  success  in  every  application.  Success  depends  on  careful  matching  of  the 
technology  to  the  operational  requirements.  Such  matching  can  only  be  done  through 
understanding  both  the  capabilities  and  limitations. 

What  follows  is  a  summary  of  some  of  the  important  lessons  learned  in  development  of 
vector  sensors  for  underwater  systems.  The  list  of  issues  is  not  exhaustive  nor  are  all  the 
issues  of  equal  importance  but  these  issues  should  be  considered  whenever  new  applications 
are  explored. 


II.  FUNDAMENTALS 


Much  of  the  potential  and  many  of  the  limitations  of  vector  sensors  are  clear  from  the 
fundamental  equations  for  linear  acoustic  fields.  Euler’s  equation  or  Newton’s  Law  per  unit 
volume  when  there  is  no  mean  flow  is 
vt  dv 

-V P  =  p—  =  jcopv  .  (1) 

at 

This  vector  equation  is  actually  three  scalar  equations.  Any  one  of  those  equations 
shows  that  a  measurement  of  one  component  of  the  pressure  gradient  is  equivalent  to  a 
measurement  of  the  acoustic  acceleration  in  that  direction.  The  term  on  the  left  represents  a 
gradient  measurement  while  the  center  term  or  the  term  on  the  right  represents  an  inertial 
measurement  (hence  the  distinction  between  gradient  and  inertial  sensors).  The  two  are 
equivalent  as  long  as  the  field  is  entirely  acoustic  in  origin ;  the  equation  is  not  true  in 
general. 


The  equation  of  continuity  (merged  with  the  equation  of  state)  is 


1  dp  _  j  (0  p 
pc2  dt  pc2 


(2) 


This  is  a  single  scalar  equation.  (Multiply  Eq.  2  by  p/j  (p  to  convert  it  to  an  expression  of 
the  potential-energy  density.  In  this  form,  it  is  clear  that  only  simple  fluid-A,lK£  compression 
is  accounted.)  If  we  know  all  three  components  of  the  divergence  of  the  velocity  field,  we 
can  determine  the  pressure  as  long  as  the  field  is  entirely  acoustic  in  origin. 


The  curl  of  Euler’s  equation, 

-VxVp  =  0  -  p— (Vxv)  -  jcopV  xv  (3) 

dt 


shows  that  the  rate  of  change  of  the  curl  of  acoustic  velocity  is  zero  as  long  as  the  field  is 
entirely  acoustic  in  origin.  Since  any  constant  component  cannot  contain  signal  information, 
this  is  equivalent  to  the  statement  that  the  curl  of  acoustic  velocity  is  zero.  This  is  indicative 
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more  of  a  limitation  in  the  linear  acoustic  approximation  than  it  is  a  statement  of  an 
experimentally  relevant  relationship.  Here  is  the  first  direct  evidence  that  non-acoustic  flows 
are  excluded.  Turbulence  -  a  source  of  flow  noise  -  is  excluded  by  the  underlying 
assumptions,  but  turbulence  is  critical  to  understanding  the  limits  to  inertial  sensing  in  a 
fluid.  In  the  derivation  of  the  linear  acoustics  equations,  it  is  assumed  that  there  is  no  bulk 
flow  or  viscosity,  that  there  is  a  specific  equation  of  state  entirely  dependent  on  the  acoustic 
compressibility  of  the  medium,  and  that  all  quantities  are  small. 

What  we  learn  from  these  equations  is  that  if  we  can  measure  pressure  gradient  we  can 
make  inferences  about  the  acoustic  particle  velocity.  Conversely,  if  we  can  measure  acoustic 
particle  velocity  then  we  can  make  inferences  about  pressure.  But  these  inferences  follow 
Eqs.  1  and  2  only  if  the  field  is  acoustic.  Structure-borne  vibration  or  turbulence  (or  even 
mean  flow)  can  invalidate  the  process.  Even  if  the  field  is  entirely  acoustic,  care  must  be 
exercised  regarding  assumptions.  We  often  assume,  with  little  error,  that  the  pressure  from 
an  acoustically  small  source  varies  as  the  reciprocal  of  distance  from  that  source.  This  does 
not  mean  that  it  is  safe  to  assume  the  same  dependence  for  the  radial  component  of  particle 
velocity.  While  the  fundamental  governing  equations  are  reliable  for  linear  acoustic  fields, 
other  relations  like p  =  pcv  are  not  valid  in  general  and  can  introduce  significant  error  (for 
example)  in  calibration  or  application  of  vector  sensors.  It  is  also  important  to  recognize  that 
these  relations  are  specific  to  fluid  media.  If  we  introduce  a  medium  that  has  any  shear 
stiffness,  we  must  often  account  for  motion  coupled  into  shear. 


III.  SENSOR  TYPES 

A.  Pressure  sensor 

The  pressure  sensor  produces  an  output  that  is  directly  related  to  the  incident  pressure. 
Ideally,  there  is  no  sensitivity  to  acceleration,  so  motion  of  the  sensor  body  does  not  produce 
signal.  This  is  fundamentally  a  point  measurement.  If  the  responsivity  of  the  sensor  could  be 
held  constant,  reducing  the  size  of  the  sensor  would  not  affect  the  output.  There  is  no 
information  about  direction  in  the  pressure-sensor  output.  This  is  a  scalar  or  “zero¬ 
dimensional”  sensor. 

B.  Pressure-gradient  sensor 

Subtraction  of  the  outputs  from  two  elements  produces  a  result  that  is  directly  related  to 
the  pressure  gradient.  The  air-acoustic  intensity  probe  functions  this  way  as  do  most  noise¬ 
canceling  microphones.  Some  early  directional  sonobuoys  used  gradient  measurement.  The 
“sensor”  consists  of  at  least  two  pressure-sensing  elements  (or  both  sides  of  a  single 
element).  The  elements  can  be  rigidly  fixed  in  space:  motion  is  not  central  to  sensing.  The 
gradient  sensor  implements  the  finite-difference  equation 

/>(x  +  <5)  -  p(x-8 )  ~  2  <5—^.  (4) 
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The  pressure  difference  is  directly  proportional  to  acoustic  acceleration,  and  the  magnitude  of 
the  resultant  is  intimately  related  to  the  separation,  2  5,  between  elements.  This  is  an 
“aperture”  sensor.  Its  performance  is  fundamentally  associated  with  its  spatial  extent. 
Reducing  the  size  reduces  the  response  (whereas  increasing  the  size  increases  the  error)  even 
if  the  element  response  is  unchanged.  This  results  in  an  unavoidable  loss  in  dynamic  range. 
The  gradient  sensor  is  fundamentally  directional. 

C.  Inertial  sensors 

The  output  of  an  inertial  sensor  is  directly  related  to  the  motion  of  the  sensor  body. 
Depending  on  the  specific  type  of  transduction,  the  output  may  be  more  nearly  proportional  to 
acceleration,  velocity,  or  displacement,  but  this  aspect  is  not  fundamental.  More  important: 
motion  of  the  sensor  body  is  essential.  Ideally,  there  is  no  sensitivity  to  pressure.  The  inertial 
measurement  is  fundamentally  a  point  measurement.  The  size  of  the  sensor  body  does  not 
affect  the  response  magnitude  in  a  fundamental  way  (but,  like  the  pressure  sensor,  does  affect 
the  output  in  practical  ways).  The  inertial  sensor  responds  to  motion  regardless  of  origin:  it 
cannot  discriminate  between  acoustic  waves,  turbulence,  or  structure-borne  vibration. 

D.  Inertial  gradient  sensors 

In  analogy  with  pressure-gradient  sensors,  it  is  also  possible  to  measure  the  spatial  rate 
of  change  of  the  velocity  field.  Two  inertial  elements  spaced  some  distance  apart  can 
provide  one  component.  In  contrast  to  the  pressure-gradient,  there  are  nine  unique  elements 
of  the  spatial  derivative  of  velocity: 


dvy 

dx 

dx 

dx 

dv* 

dvy 

dv^ 

dy 

dy 

dy 

dvx 

dvy 

dvz 

dz 

dz 

dz 

The  inertial  pair  implements  the  following  finite-difference  equation: 


vx(x  +  <5)  -  vx(x-S) 


2S  dv*  _  d2  P 
dx  cop  dx2 


(6) 


As  with  the  pressure-gradient  sensor,  the  output  depends  fundamentally  on  the 
separation  between  elements,  and  there  is  an  unavoidable  loss  in  dynamic  range.  The  sum  of 
the  main  diagonal  elements  of  the  derivative  matrix  is  the  divergence  of  the  velocity  field, 
which  is  directly  proportional  to  the  second  time  derivative  of  the  pressure.  This  provides  a 
useful  test  of  the  fidelity  of  the  sensor  system  or  of  the  acoustic  nature  of  the  field. 
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In  a  purely  acoustic  field,  the  matrix  would  be  symmetric  since  symmetry  is  equivalent 
to  zero  curl.  In  many  important  cases  of  simple  waves  (and  appropriate  alignment  of  the 
reference  axes  with  the  field)  all  of  the  off-diagonal  terms  are  zero  (again  providing  a  test  of 
either  the  sensor  system  or  the  field).  This  would  be  true  for  a  pure  plane  wave  in  any  set  of 
Cartesian  coordinates  or  for  a  spherical  wave  in  coaxial  spherical  coordinates.  (An  inertial 
rotation  sensor  would  measure  quantities  related  to  the  off-diagonal  terms.) 


IV.  REQUIREMENTS  FOR  INERTIAL  SENSORS 

A.  Fundamental  requirements 

By  far  the  most  important  aspect  of  the  inertial  sensor  is  that  it  must  move  freely  with 
the  acoustic  fluid  motion.  Furthermore, 

•  it  must  respond  faithfully  in  magnitude  and  phase  to  the  local  acoustic  particle 
velocity; 

•  its  presence  must  not  perturb  the  acoustic  field  or  that  perturbation  must  be 
correctable;  and, 

•  it  must  be  insensitive  to  or  shielded  from  non-acoustic  excitation  such  as 
turbulence  in  the  local  flow  field,  turbulence  produced  by  the  sensor  body, 
vibrations  conducted  through  the  support  structure,  or  local  radiation  of  nearby 
structural  vibration. 

Intuition  developed  with  pressure  sensors  is  dangerous  when  applied  to  inertial  sensors. 

Inertial  sensors  are  typically  far  more  sensitive  to  motion  than  pressure  sensors. 

B.  Sensitivity  to  motion 

Sensitivity  to  motion  of  a  pressure  hydrophone  is  normally  given  by  its  acceleration 
sensitivity,  Sa.  Typical  acceleration  sensitivity  for  a  hydrophone8  is  of  the  order  of  one  pascal 
per  meter  per  second-squared.  Consider  a  pressure  hydrophone  and  an  inertial  sensor 
immersed  in  an  acoustic  plane  wave.  Assume  that  both  are  free  to  move  with  the  acoustic  field 
and  that  the  pressure-associated  output  of  the  hydrophone  and  the  velocity  output  of  the 
inertial  sensor  are  accurate.  The  pressure-equivalent  signal  produced  by  the  hydrophone  in 
response  to  the  acoustic  acceleration  would  be  Sa  times  the  acoustic  velocity  multiplied  by  the 
radian  frequency,  co.  The  pressure-equivalent  signal  produced  by  the  inertial  sensor  would  be 
pc  times  the  acoustic  velocity.  The  ratio  of  the  pressure-sensor  response  to  the  velocity-sensor 
response  is  then  Saco/  pc.  This  ratio  is  4  •  10'6  at  1  Hz  and  4  •  10'3  at  1000  Hz.  In  this  situa¬ 
tion,  the  pressure  sensor  is  substantially  less  sensitive  to  motion  at  all  but  very  high  frequency. 

An  important  consequence  is  that  isolation  mounts  developed  for  pressure  hydrophones 
are  entirely  inadequate  for  inertial  sensors  from  the  standpoint  of  isolation  from  structural 
vibration.  (And,  for  an  inertial  sensor,  the  mount  must  allow  freedom  of  motion  with  the 
acoustic  field.) 
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Another  instructive  comparison  concerns  the  response  to  vertical  motion.  Vertical 
oscillatory  motion  is  fundamentally  different  from  horizontal  motion  for  underwater  sensors. 
In  vertical  motion,  there  is  a  potentially  large  pressure  signal  produced  simply  from  motion 
in  the  hydrostatic  pressure  gradient.  For  an  amplitude,  h,  of  vertical  motion,  there  is  a 
corresponding  amplitude,  pgh,  of  pressure  change.  One  millimeter  of  vertical  motion 
corresponds  to  10  pascals  of  pressure  amplitude  (140  dB  with  respect  to  one  micropascal). 
The  amplitude  of  such  an  oscillation  with  velocity,  v,  is  v/  CO. 

For  a  pressure  hydrophone  and  an  inertial  sensor  undergoing  the  same  vertical 
oscillation,  their  outputs  converted  to  equivalent  pressure  can  be  compared.  The  effective 
pressure  amplitude  is  pg  v/co.  The  pressure-equivalent  output  of  the  velocity  sensor  is  pc  v. 

o  z: 

The  ratio  of  pressure-sensor  output  to  velocity-sensor  output  is  10'  at  1  hertz  or  10'  at  1000 
Hz.  Again,  the  amplitude  of  the  spurious  signal  is  far  larger  on  the  inertial  sensor  than  on  the 
pressure  sensor.  For  example,  even  with  an  excellent  suspension  system,  it  is  difficult  to 
make  low-frequency  vertical-axis  inertial  measurements  with  a  sonobuoy. 

A  practical  consequence  of  the  motion  sensitivity  of  the  inertial  sensor  is  that  a  complex 
surface-motion  isolation  system  is  a  necessity  for  directional  sonobuoys  even  though  low  cost 
is  critical.  Furthermore,  any  application  that  specifies  acceleration-canceling  hydrophones 
would  be  challenging  to  adapt  to  inertial  sensors. 

C.  Multi-axis  systems 

If  more  than  one  component  of  the  velocity  or  acceleration  vector  is  measured  at  the 
same  location,  some  care  is  required  in  the  sensor  design.  The  design  of  mounts  and 
suspensions  is  particularly  difficult  when  the  sensor  body  must  be  free  to  move  in  all  three 
axes  as  discussed  in  a  separate  section  below.  Another  issue  associated  with  multi-axis 
sensors  is  inter-element  phase  error.  Three-axis  accelerometers  are  commercially  available 
and  consist  of  a  triad  of  orthogonal  single-axis  accelerometers  in  a  single  mounting  block.  A 
similar  design  for  an  underwater  sensor  leads  to  a  frequency-dependent  phase  error  in  the 
relationship  between  components  because  the  response  centers  of  the  elements  are  not 
collocated.  (This  phase  error  also  exists  in  the  commercial  accelerometer,  but  most 
structural-vibration-measurement  applications  do  not  involve  component-to-component 
phase  reconstruction,  so  the  small  package  size  is  more  desirable  than  phase  accuracy.) 

If,  instead  of  three  elements,  six  elements  are  arranged  so  that  each  pair  straddles  a 
common  center  point,  then  the  phase  error  associated  with  the  response  centers  is  zero.  The 
common  phase-center  design  is  also  advantageous  in  that  the  center  point  can  be  collocated 
with  the  center  of  mass  of  the  sensor  body.  If  the  phase  center,  the  center  of  mass,  and  the 
center  of  acoustic  pressure  are  all  collocated,  several  potential  errors  in  fidelity  can  be 
avoided.  The  elements  in  the  directional-sonobuoy  sensor  body  have  the  common-centroid 
design  (but  the  separate  pressure  hydrophone  is  vertically  offset  from  the  phase  center  of  the 
directional  elements).  Another  illustration  of  the  common-centroid  design  is  the  miniature 
directional  element  shown  in  Fig.  1.  Here,  four  differential  capacitance  elements  are 
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arranged  with  a  common  center  and  connected  electrically  as  a  full  bridge.  There  are  a 
number  of  practical  advantages  in  this  design  besides  the  common  phase  center.9 


FIG.  1 .  Full-bridge  differential-capacitance  inertial  sensor  for  two-axis  operation.  Each  axis  comprises  a  full 
bridge  made  from  two  opposing  structures.  The  two  orthogonal  pairs  of  structures  provide  two-axis  response 
with  a  common  phase  center.  The  structure  is  7  x  7  x  3  millimeters  in  size.  Each  bridge  is  driven  with  a  high- 
frequency  AC  signal  and  the  output,  modulated  by  the  sensed  acceleration,  is  detected  synchronously. 


V.  REQUIREMENTS  FOR  PRESSURE-GRADIENT  SENSORS 

A.  Fundamentals 

The  pressure-gradient  sensor  relies  (either  explicitly  or  implicitly)  on  subtraction  of  two 
nearly  equal  quantities.  This  concept  is  essential  to  its  operation  and  is  the  root  of  its  most 
serious  limitations.  The  gradient  sensor 

■  must  have  excellent  element-to-element  matching  in  both  phase  and  amplitude  (or 
have  a  reliable  calibration  and  correction); 

■  must  have  low  self-noise  in  all  elements  since  the  dynamic  range  is  compromised 
fundamentally  by  the  subtraction;  and 

■  must  not  perturb  the  acoustic  field  significantly. 

B.  Finite  differences 

The  magnitude  of  the  gradient  calculated  by  finite  difference  is  reduced  in  amplitude  by 
a  factor  of  the  order  of  the  element  separation,  d,  in  wavelengths.  Small  kd  (where  k  is  the 
wave  number)  is  essential  for  accuracy  in  determining  the  derivative  but  small  kd  also  means 
a  large  reduction  in  dynamic  range.  This  is  completely  analogous  to  the  loss  of  precision  in 
numerical  finite-difference  operations. 

There  are  also  phase  errors  associated  with  the  orientation  of  the  elements.  A  sensor 
designed  to  measure  the  gradient  in  two  dimensions  could  be  made  from  three  elements  as 
shown  in  Fig.  2.  However,  the  phase  centers  of  both  pairs  and  of  the  equivalent  monopole 
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are  not  coincident.  If,  instead,  four  elements  are  used,  then  all  of  the  phase  centers  can  be 
collocated.  This  construction  is  equivalent  to  using  centered  differences  instead  of  one-sided 
differences  in  finite-difference  computations. 


FIG.  2.  Two  configurations  for  measuring  the  gradient  in  two  dimensions.  For  the  arrangement  on  the  left,  the 
phase  centers  of  the  two  dipoles  and  of  the  equivalent  pressure  formed  by  summing  all  three  elements  are  not 
coincident.  In  contrast,  the  phase  centers  for  both  dipoles  and  the  equivalent  pressure  for  the  arrangement  on 
the  right  are  coincident.  The  four-sensor  arrangement  has  zero  in-plane  phase  error. 


C.  Basic  dipole 


The  basic  “dipole”  element  is  shown  in  Fig.  3.  This  is  only  a  true  dipole  for 
vanishingly  small  separation,  but,  in  practice,  its  output  is  not  essentially  different  than  a 
dipole  because  the  spacing  is  kept  small.  If  the  outputs  from  the  two  elements  are  subtracted, 
the  resultant  pressure  is 


Pd 


^■Pa-^Pb 


4y>osin 


kd 


cos  6 


\ 

J 


(7) 


~  2  j  p0  k d cos 6  ( kd  «  l) 


d 

FIG.  3.  Basic  two-element  gradient  sensor.  The  outputs  of  the  two  elements  are  subtracted. 

For  separation,  d,  much  less  than  an  acoustic  wavelength,  the  system  behaves  as  a  dipole. 

In  order  to  facilitate  comparison  of  monopole,  dipole,  and  quadrupole  sensors,  a 
consistent  normalization  of  four  units  (two  for  each  element  in  this  case)  will  be  used 
throughout.  The  comparison  monopole  would  have  a  scale  factor  of  four.  This  normalization 
maintains  the  same  effective  self-noise  level  for  each  type  of  sensor  so  that  the  progressive 
loss  in  dynamic  range  is  immediately  evident. 

Of  critical  importance  in  understanding  the  performance  of  the  gradient  sensor  is  the 
influence  of  noise.  Internal  self-noise  in  the  elements  is  normally  not  coherent  from  element 
to  element  (in  contrast  to  ambient  noise).  Consequently,  even  though  the  outputs  are 
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subtracted,  the  self-noise  adds  incoherently.  For  the  dipole,  the  signal  amplitude  is  2  kd 
while  that  of  the  monopole  is  4  but  the  self-noise  floor  (with  our  normalization)  is  the  same. 
As  a  result,  the  signal-to-noise  can  be  poor  even  if  each  element  taken  independently  has 
high  signal-to-noise. 

Within  these  constraints,  three  orthogonal  dipoles  (using  six  elements,  for  example) 
give  all  of  the  first-order  terms  in  the  Taylor  series  expansion  for  pressure  about  the 
measurement  point  and  completely  define  the  vector  gradient  in  pressure. 

The  outputs  of  two  orthogonal  dipoles  (and  a  simple  pressure  measurement)  can  be 
combined  to  form  a  cardioid,  which  can  be  steered  in  any  direction  by  proper  phasing  of  the 
two  dipoles: 


Pd  +  ikdPu  ~  4.JPo  kd[l  +  cos#]  (8) 

where  pm  is  four  times  po  to  maintain  the  “four-unit”  normalization.  Notice  that  we  retain  the 
signal-to-noise  disadvantage  (kd  factor).  In  addition,  a  frequency-independent  phase  shift  of 
90  degrees  must  be  provided. 

D.  Basic  quadrupole 

If  a  dipole  is  useful,  then  perhaps  a  quadrupole  is  even  more  useful.  In  finite-difference 
form,  the  longitudinal  quadrupole  is  equivalent  to  the  second  spatial  derivative: 

p(x  +  8)  -  2 p(x)  +  p(x-8)  ~  4 82  .  (9) 

ox 

The  arrangement  for  a  longitudinal  quadrupole  is  shown  in  Fig.  4.  If  the  outputs  are 
combined  as  follows 


Pq 


=  2Pb 
“  4  p0 


Pa  ~  Pc 
f  kd 


V  2  J 


cos  0 


2 p0  -  2p0  cos(kdcosO) 
(kd  «  l) 


(10) 


then  we  obtain  an  approximation  to  the  second  derivative  of  the  field.  The  directivity  (see 
Fig.  5)  is  sharper  than  the  dipole,  but  the  signal-to-noise  penalty  (in  amplitude)  goes  as  (kd)2, 
so  the  dynamic  range  is  seriously  limited.  From  Fig.  5,  if  the  single  element  noise  is  60  dB 
below  the  signal,  the  dipole  peak  is  less  than  30  dB  above  the  noise  and  the  quadrupole 
would  produce  no  directivity  at  all.  In  effect,  the  quadrupole  sacrifices  almost  70  dB  of 
dynamic  range  (11  bits  in  a  digital  system)  and  would  require  an  extraordinarily  high  element 
signal-to-noise  ratio  to  estimate  the  second  derivative  accurately. 
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FIG.  4.  Basic  configuration  of  a  three-element  sensor  for  estimating  the  second  space 
derivative.  Two  times  the  output  of  the  center  element  is  subtracted  from  the  sum  of  the 
outputs  of  the  outer  elements.  If  d  is  much  less  than  a  wavelength,  the  configuration  is  a 
longitudinal  quadrupole. 
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FIG.  5.  Comparison  of  the  directional  response  of  a  monopole,  a  dipole,  and  a  quadrupole 
(/ ka  =  0.05).  The  sensors  are  compared  in  such  a  way  as  to  have  equal  self-noise  power,  so 
the  reduction  in  level  of  the  main  lobe  corresponds  to  the  sacrifice  in  dynamic  range. 


One  hindrance  to  the  development  of  good  engineering  judgment  in  directional  systems 
is  the  tendency  to  normalize  directivity  plots  so  that  the  maximum  response  is  unity.  Fig.  5  is 
constructed  so  that  the  self-noise  floor  is  identical  in  each  case  and  the  performance  penalty 
incurred  in  the  finite-difference  process  is  clearly  evident.  For  transmitting  arrays, 
normalization  to  input  power  produces  more  instructive  comparisons  than  peak 
normalization.  For  example,  comparisons  between  the  beam  patterns  of  the  primary  and 
difference  beams  in  a  parametric  sonar  done  in  this  manner  would  still  reveal  the  beam  shape 
but  also  highlight  the  inefficiency  of  the  process. 

The  signals  from  orthogonal  quadrupoles  can  be  combined  to  form  a  cardioid-squared 
pattern  as  follows 

4 pQ  —  ljkdpD  +  (kd)2  pM  ~  Ap0  (kd)2  [l  +  cos#]2  (11) 

which  can  be  steered  to  any  orientation  by  proper  phasing. 
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Because  of  the  dynamic-range  loss,  the  quadrupole  may  be  completely  impractical  with 
pressure  elements  and  second-order  finite-difference  construction;  however,  the  equivalent 
can  be  formed  through  finite-difference  combination  of  inertial  sensor  outputs  (see  Eq.  6). 
This  approach  avoids  the  harsh  (kdy  dynamic-range  penalty  but  introduces  the  motional 
sensitivity  of  the  inertial  sensor  that  is  exacerbated  by  the  finite-difference  operation. 


VI.  COMPARISON  OF  GRADIENT  AND  INERTIAL  SENSORS 

A.  Inertial  sensors 

The  sensor  output  is  directly  related  to  acceleration  (or  velocity  or  displacement),  and 
good  dynamic  range  is  readily  achievable.  The  effective  noise  floor  of  an  inertial  sensor  can 
usually  be  designed  to  be  far  below  that  of  an  equivalent  gradient  sensor.  However, 

■  the  sensor  body  must  be  free  to  respond  accurately  to  the  acoustic  wave; 

■  it  must  be  protected  from  non-acoustic  disturbances; 

■  it  is  difficult  to  miniaturize  sufficiently  to  make  good  measurements  at  high 
frequency;  and, 

■  it  is  difficult  to  isolate  sufficiently  to  make  good  measurements  at  low  frequency. 

B.  Pressure-gradient  sensors 

The  elements  of  the  gradient  sensor  can  be  compact  and  rigidly  fixed  in  space.  The 
elements 


■  must  be  well  matched  in  both  phase  and  amplitude,  or  some  provision  for 
calibration  and  correction  must  be  made; 

■  must  have  low  self  noise;  and, 

■  must  have  a  large  basic  dynamic  range  to  compensate  for  the  subsequent  loss  in 
range  associated  with  the  finite-difference  operation. 


VII.  SUSPENSIONS 
A.  Basic  requirements 

Suspensions  are  critical  to  proper  operation  of  inertial  sensors  much  more  so  than  for 
pressure  or  pressure-gradient  sensors  so  the  focus  here  is  on  inertial  sensors.  The  fundamentals 
of  inertial- sensor  suspensions  have  been  covered  in  detail  elsewhere10  so  only  a  brief  summary 
of  the  basics  is  presented.  Fundamentally,  the  suspension 

■  must  have  a  natural  frequency  well  below  the  intended  range  of  operation; 

■  must  fix  the  time-average  position  and  orientation  of  the  sensor  body; 

■  must  permit  movement  of  the  sensor  body  with  the  acoustic  field; 
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■  must  isolate  the  sensor  from  structure-borne  vibration; 

■  must  not  distort  the  response  of  the  sensor  either  in  magnitude,  phase,  or  apparent 
angle;  and 

■  must  withstand  operational  shock. 

In  order  to  categorize  various  systems,  we  will  consider  both  the  dimensionality  of 
the  sensor  and  the  number  of  position  and  orientation  constraints.  Since  pressure  sensing 
produces  a  scalar  measure  of  the  field,  that  will  be  called  a  “zero-dimensional”  measurement. 
A  single-axis  vector  probe  constitutes  a  one-component  measurement  and  so  on  to  three 
dimensions  of  measurement.  More  dimensions  can  be  defined  for  inertial-gradient  sensors, 
but  we  can  decompose  these  higher  order  systems  into  sets  of  one-,  two-,  or  three-  component 
sensors  insofar  as  suspension  design  is  concerned.  We  will  not  discuss  mounting  zero- 
dimension  sensors  further. 

Since  the  suspension  must  support  the  sensor  in  some  repeatable  position,  it  imposes 
constraints  at  least  on  the  time-averaged  position.  There  are  six  degrees  of  freedom  for  such 
support:  three  axes  of  translation  and  three  axes  of  rotation.  Depending  on  the  application,  it 
may  not  be  necessary  to  constrain  all  six  degrees  of  freedom. 

B.  One-component  sensors 

If  the  vector  sensor  is  intended  to  measure  only  one  component  of  the  field,  the 
suspension  design  is  straightforward  with  a  number  of  practical  variants.  The  ribbon 
microphone  and  the  modern  air-acoustic  intensity  probes  are  examples  of  one-component 
sensors.  Presently  under  commercial  development11  are  inertial  velocity  sensors  with 
compact  suspensions  of  spiral-cut  leaf  springs.  By  making  use  of  the  fact  that  low  stiffness  is 
not  required  in  directions  orthogonal  to  the  sense  axis,  the  suspension  is  much  more  space 
effective  and  durable.  Such  a  suspension  introduces  anisotropy,  which  can  have 
consequences  for  reconstruction  of  the  acoustic  field. 

Single-axis  sensors  embedded  in  compliant  coatings  rely  on  the  characteristics  of  the 
coating  for  suspension.  Freedom  of  motion  along  the  sense  axis  depends  on  the  relatively  low 
shear  stiffness  of  the  coating.  Perpendicular  to  the  sense  axis  and  in  the  plane  of  the  coating, 
the  compressional  elastic  properties  of  the  coating  material  are  important.  If  the  compressional 
stiffness  is  close  to  that  of  water  and  the  boundaries  are  far  removed,  then  reasonable  isotropy 
and,  therefore,  fidelity  in  one  dimension  is  possible.  Two  varieties  of  single-component 
suspension  (and  one  type  of  two-component  suspension)  are  shown  in  Fig.  6. 

C.  Two-component  sensors 

In  underwater  acoustics,  two-component  sensing  is  far  more  common  than  any  other 
type.  For  two-component  sensors,  the  third  axis  can  be  used  as  a  “stiff’  support  axis.  This 
results  in  a  robust  suspension  with  a  very  low  resonance  frequency  in  either  sense  direction. 
The  sonobuoy  (see  Fig.  6)  uses  a  uniaxial  support:  a  weight  hangs  below  the  sensor  body  and 
pulls  the  connecting  cable  taut  between  sensor  body  and  whatever  is  above  the  sensor  body. 
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This  taut-line  suspension  is  very  weak  for  small  transverse  displacements,  so  the  sensor  body 
is  free  to  follow  the  acoustic  motion  with  good  fidelity  above  the  transverse  resonance. 

There  are  three  translational  constraints  but  only  two  rotational  constraints.  Rotation  about 
the  support  axis  is  not  constrained.  Instead,  a  magnetic  compass  is  used  to  stabilize  the 
azimuthal  orientation  electronically. 

Good  low-frequency  performance  is  difficult  to  obtain  with  any  inertial  sensor  because 
the  amplitudes  of  the  corrupting  excitations  tend  to  increase  with  decreasing  frequency.  The 
sonobuoy  suspension  is  designed  to  isolate  from  vertical  and  horizontal  surface  motion  but  it 
cannot  remove  steady  horizontal  flow  that  results  from  current  shear.  Consequently,  there  is 
some  irreducible  level  of  flow  past  the  sensor  body  that  results  in  flow-induced  noise. 
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FIG.  6.  Primary  varieties  of  velocity-sensor  suspension.  The  suspended  (sonobuoy)  two-component  sensor  is 
shown  on  the  left.  A  sensor  supported  by  mechanical  springs  at  each  end  is  shown  in  the  center,  and  a  sensor 
embedded  in  a  compliant  suspension  material  is  shown  on  the  right. 


D.  Three-component  sensors 

Three-component  sensors  with  six-degree-of-freedom  constraint  are  difficult  to  build 
with  reasonable  size.  Probably  the  best  known  three-component  sensor  system  is  the 
Swallow  Float12  system.  This  system  has  only  one  translational  constraint  provided  by 
buoyancy  adjustment  for  stability  in  the  ocean  density  profile.  Two  rotational  constraints  are 
imposed  by  shifting  the  center  of  mass  below  the  center  of  buoyancy;  rotation  about  the 
vertical  axis  is  unconstrained.  By  relaxing  the  translational  constraint,  exceptional  low- 
frequency  performance  is  obtained. 

In  order  for  a  system  like  the  Swallow  Float  to  reach  vertical  equilibrium  in  a  density 
gradient,  the  body  must  be  significantly  less  compressible  than  the  fluid.  As  a  result,  the 
body  will  scatter  the  local  acoustic  field  substantially  if  the  frequency  of  operation 
approaches  ka=  1 . 

Another  approach  for  three-component  sensing  is  to  stage  single-component  sensors  in 
nested  suspensions.  This  is  difficult  to  accomplish  without  a  serious  volume  penalty. 
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A  third  technique  is  to  embed  the  sense-element  cluster  in  a  compliant,  elastic  material. 
Special  care  must  be  exercised  in  the  design  of  the  shear  stiffness  and  the  external  mounts, 
however.  If  the  design  relies  on  shear  stiffness  to  control  the  suspension  resonance,  the 
mount  must  be  designed  to  avoid  any  axis  being  constrained  by  compression  of  the  material 
since  the  compressional  modulus  is  normally  much  higher.  With  such  a  system,  it  is 
extremely  difficult  to  avoid  considerable  anisotropy  in  the  suspension. 

E.  Symmetry  And  Anisotropy 

The  fundamental  tenet  of  suspension  design  is  that  the  primary  resonance  of  the 
suspension  should  be  well  below  the  lowest  frequency  of  interest.  While  this  is 
straightforward  to  achieve  in  one-  and  two-component  vector  sensors,  it  can  be  very  difficult 
to  achieve  in  practical  three-component  sensors.  If  all  of  the  fundamental  resonances  are  not 
well  below  the  frequency  of  operation,  then  we  must  consider  the  effects  of  three- 
dimensional  anisotropy  in  the  suspension  on  the  signals  produced  by  the  sensor. 

For  example,  the  suspension  used  in  the  sonobuoy  directional  package  has  resonances 
with  respect  to  horizontal  displacements  that  are  very  low,  but  the  resonance  with  respect  to 
vertical  displacement  can  be  in  the  band  of  interest.  This  anisotropy  affects  the  fidelity  of  the 
sensor  response  with  respect  to  the  vertical  component.  In  the  sonobuoy  application, 
accurate  resolution  of  the  vertical  component  is  immaterial.  The  anisotropy,  in  this  case, 
only  serves  to  increase  the  apparent  level  of  signals  at  low  frequency  over  that  that  would  be 
produced  by  a  perfect  horizontal-component  sensor. 

However,  this  anisotropy  would  be  problematic  if  it  were  necessary  to  sense  the  vertical 
component  with  the  same  suspension.  Not  only  would  the  isolation  be  poor,  but  also  the 
vertical  sensor  would  be  overly  constrained  and  the  vector  resultant  would  have  the  wrong 
direction.  The  Swallow  Float  and  the  neutrally  buoyant  infrasonic  DIFAR  sonobuoy  with 
their  quasi-free  “suspensions”  circumvent  this  problem  but  only  by  incurring  substantial 
limitation  in  employment. 

A  soft  polymeric  material  can  be  used  as  a  suspension  if  the  restoring  force  is 
arranged  to  be  through  shear  in  the  polymer.  If  the  polymer  body  is  fixed  at  two 
diametrically  opposite  points,  then  shear  is  the  operative  restoring  force  for  displacement 
perpendicular  to  the  support  axis.  However,  along  the  support  axis,  the  restoring  force  is 
compression/extension.  The  modulus  may  be  much  higher  for  compression/extension,  so  it 
is  more  difficult  to  keep  this  resonance  low.  Consequently,  the  isolation  may  be  much  worse 
and  the  motion  more  restricted  for  the  sensor  aligned  along  the  support  axis.  Clearly,  this 
also  results  in  anisotropy  in  the  overall  suspension  scheme  and  distortion  of  angle  of  arrival. 
In  evaluating  such  systems,  it  is  vital  to  test  the  sensor  for  arrivals  from  all  directions.  If  it  is 
only  evaluated  for  arrivals  perpendicular  to  the  support  axis  (for  example),  problems  can  go 
undetected.  The  sensor  should  be  tested  in  its  eventual  mount  and  should  be  evaluated  for  all 
directions  of  arrival  not  only  with  respect  to  the  sensor  elements  but  also  with  respect  to  the 
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mount.  (It  is  also  necessary  to  measure  the  excitation  of  the  mount  itself  to  avoid  confusing 
sensor  suspension  problems  with  problems  associated  with  laboratory  fixtures.) 

Interesting  problems  can  also  arise  when  periodic  arrays  of  single-component  sensors 
are  embedded  in  a  compliant  layer.  The  suspension  problems  for  a  single  sensor  are  not 
essentially  different  from  those  discussed  above,  but  the  periodic  array  adds  a  complication. 
The  interaction  of  the  sensor  bodies  through  either  shear  or  compressional  coupling  in  the 
compliant  layer  introduces  effective  bulk  properties  to  the  layer.  The  periodic  structure 
introduces  “conduction  bands”  and  “stop  bands”  for  propagation  through  the  layer. 
Coincidence  between  wave  numbers  in  the  material  and  incident  tangential  wave  numbers  for 
acoustic  waves  produces  strong  coupling.  The  conduction/stop  band  structure  can  make  this 
coupling  complicated  and  introduce  spurious  structure  to  the  vector  sensor  output. 


VIII.  PERFORMANCE  LIMITS 

There  are  many  limits  to  performance  in  any  sensor  technology.  We  have  already 
discussed  problems  associated  with  response  to  non-acoustic  disturbance  and  limitations 
imposed  by  suspensions.  Here,  a  few  additional  problems  are  summarized.  As  is  the  spirit 
throughout  this  paper,  this  is  not  an  exhaustive  list.  The  intent  is  to  ensure  that  several 
important  considerations  are  not  ignored  in  the  design,  development,  and  application  of 
vector  sensors. 

A.  Noise 

Sensor  self-noise  can  limit  the  performance  of  a  system  dramatically.  Those  noise 
components  associated  with  the  basic  sensor  structure,  materials,  or  first-stage  electronics 
should  never  be  overlooked  even  if  they  have  been  ignored  successfully  in  the  past.  These 
components  are  generally  uncorrelated  from  sense  element  to  sense  element  and,  as  such,  are 
straightforward  to  account  in  a  performance  budget.  The  impact  of  such  noise  on  dipole  and 
quadrupole  sensors  was  discussed  above.  In  short,  these  noise  components  should  not  be 
compared  to  the  omnidirectional  ambient-noise  level  to  determine  their  importance.  The 
more  array  gain  desired  from  the  system  or  the  more  dynamic  range  sacrificed  by  finite 
difference,  the  further  below  ambient  these  noise  sources  must  be. 

Vector  sensors  are  particularly  susceptible  to  flow  noise.  Flow  noise  is  noise  associated 
with  relative  motion  between  the  sensor  body  and  the  fluid  in  which  it  is  immersed.  Flow 
noise  may  consist  of  pre-existing  turbulence  advected  past  the  sensor  body  by  the  relative 
flow,  or  it  may  consist  of  turbulence  produced  by  the  presence  of  the  body  in  the  flow.  In 
either  case,  there  is  a  non-acoustic  field  to  which  the  sensor  responds.  The  effects  of  flow 
noise  can  be  reduced  by  increasing  the  effective  sensor  size  (with  a  wind-screen,  a  flow  shield, 
or  many  sensors  connected  together).  In  any  case,  though,  the  overall  sense  dimension  must 
be  small  enough  to  respond  properly  to  the  local  acoustic  field.  If  the  correlation  distance  of 
the  turbulence  is  much  smaller  than  an  acoustic  wavelength,  extra  sensor  dimension  provides 
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some  averaging  over  the  flow  disturbances.  Combining  the  outputs  of  several  closely  spaced 
elements  is  common  practice  in  towed  arrays  of  pressure  hydrophones. 

Any  motion  transmitted  from  a  structure  to  the  sensor  will  also  produce  a  signal  that  is 
not  associated  with  the  incident  acoustic  field.  One  of  the  critical  aspects  of  suspension 
design  is  transmissibility  or  vibration  isolation. 

For  a  “good”  sensor  design,  the  ocean  ambient  should  control  the  achievable 
performance  rather  than  the  sensor  noise,  flow  noise,  or  structure-borne  noise.  Consideration 
of  the  effects  of  ambient  noise  is  difficult,  though,  because  the  ambient-noise  field  is 
complicated.  Even  in  the  deep  ocean,  the  ambient  is  strongly  anisotropic;  in  important 
shallow-water  areas,  the  anisotropy  can  be  extreme.  In  order  to  facilitate  calculations,  it  is 
often  assumed  that  the  ambient  noise  is  isotropic.  However,  unless  the  frequency  is  above 
50  kHz  or  so,  ocean  ambient  (particularly  in  shallow  water)  can  be  strongly  anisotropic  and 
this  will  affect  performance  predictions.  Performance  numbers  obtained  using  the 
assumption  of  isotropic  noise  are  neither  upper  nor  lower  limits  to  ultimate  performance.  In 
some  cases,  considerable  performance  improvement  can  be  obtained  by  exploiting  the 
anisotropy  of  the  noise  field.13 

B.  Fidelity  of  response 

If  the  suspension  has  been  designed  well,  then  its  influence  on  the  acoustic  performance 
of  the  sensor  should  be  negligible.  However,  the  sensor  itself  can  corrupt  the  received  signal. 
Although  many  lab  measurements  and  calibration  focus  on  the  magnitude  response  of 
sensors,  the  phase  response  is  often  more  important  in  determining  the  eventual  system 
performance. 

Consider,  for  example,  a  single-degree-of-freedom  accelerometer  that  might  be  used  in 
an  inertial  sensor.  Conventional  wisdom  suggests  designing  the  Q  to  be  rather  small  so  that 
the  response  is  “flat”  over  as  wide  a  band  of  frequency  as  possible.  Such  a  design 
philosophy  results  in  substantial  phase  error  well  below  the  resonance.  While,  in  principle, 
these  phase  errors  can  be  calibrated  and  removed,  the  success  of  such  a  correction  process 
relies  on  the  stability  of  the  resonance  frequency  as  well  as  on  the  stability  of  the  damping. 
Small  shifts  in  the  resonance  frequency  (which  would  be  difficult  to  detect  from  an 
examination  of  the  magnitude  response)  can  produce  substantial  phase  changes  even  a 
decade  below  the  resonance. 

A  higher  Q  system  results  in  much  smaller  in-band  phase  error.  Shifts  in  the  resonance 
frequency  or  damping  (temperature  dependence,  for  example)  have  reduced  effect  away  from 
the  resonance.  Since  phase  errors  usually  create  more  problems  in  array  beamforming  or  in 
vector  component  resolution  (particularly  for  gradient  sensors),  the  higher  Q  design  is 
important.  The  price  of  the  higher  Q  is  that  the  system  must  handle  a  larger  dynamic  range 
of  signals. 
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An  often-overlooked  aspect  of  the  phase  response  results  from  the  low-frequency  roll¬ 
off.  The  lower  edge  of  the  pass  band  for  a  piezoelectric  sensor  (for  example)  is  normally  set 
by  the  relationship  between  the  sensor  capacitance  and  the  input  resistance  of  the  first  stage 
preamplifier.  This  is  a  single-pole  filter  so  the  effects  of  phase  extend  well  into  the  pass  band 
and  are  only  as  stable  as  the  element  values  in  the  filter. 

C.  Calibration  and  operation  in  the  “near”  field 

Often,  it  is  assumed  that  acoustic  pressure  is  equal  to  the  fluid  density  times  the  fluid 
sound  speed  times  the  acoustic  particle  velocity  but  this  is  only  true  in  certain  circumstances 
(most  commonly  a  plane-wave  field).  In  “free-field”  calibration,  the  range  between  source 
and  receiver  is  often  such  that  kr  is  not  very  large.  For  comparison  calibration  (for  example) 
of  two  pressure  sensors,  the  near-field  details  are  unimportant.  However,  for  comparison  of 
pressure  sensors  and  inertial  sensors  a  correction  is  required. 

If  the  source  is  small  with  respect  to  a  wavelength,  then  it  can  be  treated  as  a  simple 
source.  If  we  use  Euler’s  equation  to  find  the  relationship  between  pressure  and  particle 
velocity  using  the  pressure  field  produced  by  a  simple  source,  we  find  that  the  normalized 
ratio  is 


P 


(12) 


For  kr  =  12  (almost  three  meters  at  1000  Hz),  the  phase  is  still  5  degrees  different  from  the 
simplistic  plane-wave  assumption. 

While  it  is  simple  to  remove  this  correction  from  measured  data,  it  is  often  useful  to 
measure  the  phase  relationship  as  a  function  of  source-receiver  separation.  The  tangent  of 
the  measured  phase  (with  the  phase  bias  error  subtracted)  as  a  function  of  kr  should  be  a 
straight  line.  By  iterative  fitting,  the  phase  bias  can  be  determined  and  the  quality  of  the 
calibration  can  be  assessed.  Departures  from  the  expected  straight-line  behavior  should  be 
suspect. 
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Abstract:  Calculations  are  presented  that  demonstrate  that  a  single  highly 
directive  acoustic  receiver  can  have  a  directivity  index  of  up  to  9.5  dB.  This 
compares  to  a  directivity  index  (NDI)  of  6  dB  for  an  acoustic  vector  sensor;  a 
single  pressure  sensor  is  omnidirectional  and  has  no  directivity.  In  addition  to 
measuring  pressure  (p)  and  particle  velocity  (u,v,w),  as  in  a  vector  sensor,  the 
highly-directive  acoustic  receiver  measures  also  the  three  gradients  of  acoustic 
particle  velocity  (du/dx,  dv/dy,  dw/dz).  The  sum  of  these  gradients  is  propor¬ 
tional  to  the  instantaneous  density  of  the  acoustic  field.  There  is  a  price  to  pay 
with  highly  directive  sensors;  these  sensors  are  more  sensitive  to  nonacoustic 
(subsonic)  noise  sources. 

I.  INTRODUCTION 

In  1996,  the  Office  of  Naval  Research  funded  a  study  at  the  Naval  Undersea 
Warfare  Center  to  compare  the  benefits  of  a  vector-sensing  passive  sonar  array  to  a 
conventional  pressure-sensing  array.1  This  initial  investigation  showed  that  directivity 
and  self-noise  gains  could  be  realized  when  vector  sensors  are  used  in  place  of 
conventional  pressure-sensing  hydrophones.  These  sensors  can  be  used  to  create 
cardioid-directive  array  elements,  which  will  provide  6.0  dB  of  additional  directivity  gain 
(Ndi)  and  as  well  as  provide  rejection  of  noise  from  specific  directions. 

Utilizing  acoustic  intensity  measurements  in  submarine  sonars,  such  as  the  Low 
Frequency  Hull  Array,  will  be  challenging.  There  are  additional  sonar  self-noise 
mechanisms  that  are  not  present,  or  are  ineffective,  with  conventional  pressure-sensing 
hydrophones  —  though  there  are  many  potential  benefits.  For  example,  Hickling  and 
Morgan  and  Wei  present  a  formulation  to  determine  the  positions  of  sound  sources 
using  underwater  acoustic  intensity  measurements.  The  derivation  considers  two  steady- 
state  acoustic  sources,  located  at  arbitrary  positions  Q\  and  Qi  with  Cartesian  coordinates 
(xi,  yi,  z\)  and  (x2,  V2,  zj),  respectively.  Two  vector  sensors  are  used  to  determine  the 
sound  intensity  field  at  known  positions  A  and  B.  From  these  measurements,  it  is  simple 
to  derive  a  set  of  eight  nonlinear  equations  with  eight  unknowns  (the  coordinate  positions 
and  rms  pressures  of  the  two  sources).  The  task  then  is  to  solve  this  system  of  equations; 
the  authors  investigated  a  possible  solution  via  polynomial  continuation. 

II.  SIGNAL  GAIN  IN  ISOTROPIC  NOISE 

The  directivity  index  (Ndi)  of  an  array  (or  directional  sensor)  is  defined  as  a  decibel 
measure  of  the  improvement  in  the  signal-to-noise  ratio  (SNR)  that  a  beamformed  array 
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(or  sensor)  provides  in  an  ideal  isotropic  noise  field  with  a  perfectly  correlated  signal, 
relative  to  an  omnidirectional  array  element  in  the  free  field.4  That  is 

{SNR) 

DI  =  lOLogDF  =  lOLog  "‘y - ,  (1) 

yoly  K  )  omni_ element 

where  DF  is  the  directivity  factor.  This  ratio  reduces  to 
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where  2?(0s,(|)s)  is  the  power  sum  of  the  array  response  steered  to  0s,(|)s.  Hence, 
determining  the  signal  gain  of  a  directive  sensor  in  isotropic  noise  reduces  to  evaluating 
the  above  double  integral  over  0,(|). 

The  directional  dependence  of  the  acoustic  particle  velocity  components,  as  well  as 
the  velocity  gradient  components,  may  be  obtained  from  the  linearized  momentum  and 
continuity  equations.5  For  brevity,  this  derivation  will  be  omitted  here. 

Suppressing  the  harmonic  time  dependence  and  assuming  a  planewave  amplitude  of 
unity,  the  power  sum  of  the  weighted  quantities  can  be  written  as5 

5{1O)(0,(]))=|  w  +wxa  +  w  b  +  wzc+wxxa2+w  b2+wzzc2 
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where  a(0,(|))  =  cos(0)sin((|)),  /i(0,(j))  =  sin(0)sin((|)),  and  c(4>  )=  cos((|)),  and  0,4>  was  defined 
in  Fig.  1. 


The  objective  now  is  to  determine  the  weights  that  optimize  the  directivity  expression 
shown  in  Eq.  2.  It  is  straightforward  to  derive  the  directivity  of  a  single  omnidirectional 
pressure  sensor.  Namely,  from  Eqs.  2  and  3,  with  all  weights  except  wp  zeroed 
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As  expected,  the  directivity  index  is  zero  and  the  single-pressure  sensor  is 
omnidirectional.  A  similar  calculation  for  a  sensor  that  measures  only  the  u  component  of 
particle  velocity  gives 
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Upon  substitution  in  Eq.  (2), 

E)F(x)=3cos2(0s)sin2(^).  (6) 

Equation  (6)  indicates  that  a  single  velocity  sensor  may  have  a  maximum  directivity 
factor  three  times  greater  than  an  omnidirectional  pressure  sensor,  that  is,  a  gain  in 
directivity  index  of  4.8  dB. 

The  maximum  response  of  a  vector  sensor  is  given  as 

^(4)(0J»<|)5)=|wp  +  fl(05,<|)J)wJC  +b(<ds,§s)wy+c(§s)wz  |2,  (7) 

and  the  two-fold  integral  becomes 
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Hence,  the  directivity  factor  of  a  single  element  that  measures  all  three  components  of 
particle  velocity  and  acoustic  pressure  is 
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Maximum  directivity  is  obtained  with  weights 
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wz  =  3c((])J. 

(12) 

Substituting  the  optimal  real  weights  into  Eq.  (9)  yields 

DFpv  =  l  +  3(a(0„^)2  +  b(Qs,$s)2  +  c(( M2)=4. 

(13) 

Therefore,  the  maximum  directivity  for  a  single- vector  sensor  is  DI  =  10  log(4)  =  6  dB, 
and  this  holds  for  any  steering  angle  and  frequency. 

The  optimal  weights  given  in  Eqs.  (10)  through  (12)  were  obtained  by  solving  for 
the  set  of  weights  that  produced  an  extrema  value  for  DFpv.  In  practice,  the  procedure  can 
be  extended  to  determine  the  optimal  weights  for  higher-order  sensors,  that  is,  the 
weights  wxx,  Wyy,  w-z,  Wxy,  wxz,  wyz.  The  algebra  becomes  cumbersome.  The  optimal 
weights  for  a  directional  sensor  that  measures  pressure,  the  three  components  of  particle 
velocity,  and  a  single  gradient  of  velocity  (i.e.,  the  weights  wp,  wx,  wy,wz,  Wxx)  were 
determined.  These  optimal  weights  generate  the  directivity  shown  in  Fig.  2,  which  is 
noted  to  be  as  much  as  9.5  dB  (or  3.5  dB  greater  than  that  of  a  vector  sensor)  along  the 
boresight  of  the  x-axis.  Also  shown  are  the  steering  directions  at  which  the  velocity 
gradient,  u',  did  not  provide  additional  gain.  Not  surprisingly,  if  the  gradient  v'  is 
measured  as  well,  the  directivity  increases  along  the  boresight  of  the  y-axis,  at  0S  =  ±90°. 

As  with  a  vector  sensor,  it  can  be  shown  that  Bl0(Qs,§s)  simplifies  to  a  constant 
value  (independent  of  sensor  steering  direction  (0s,(j)s))  when  the  weights  are  chosen  to 
match  each  component’s  angular  response.  However,  it  has  not  been  determined  whether 
such  a  weight  set  is  optimal.  It  is  believed  that  a  highly  directional  sensor  will  achieve  a 
constant  9.5  dB  of  directivity,  independent  of  angular  steering,  with  optimal  weights. 

III.  NON-ACOUSTIC  (SUBSONIC)  NOISE  SOURCES 

Maidanik  and  Becker6  compare  the  wavevector  filtering  characteristics  of  a 
pressure  sensor  to  that  of  a  single-axis  velocity  sensor.  The  sensors  are  assumed  to  be 
embedded  in,  or  on,  the  surface  of,  a  compliant  boundary  exposed  to  noise  generated  by 
turbulent  flow.  The  results  presented  clearly  demonstrate  that  a  velocity  sensor  can  be 
more  sensitivity  to  subsonic  noise  sources  than  an  omnidirectional  pressure  sensor. 
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AZIMUTH  STEERING  ANGLE  (6° ) 

FIG.  2.  Directivity  of  a  directional  sensor  that  measures  pressure,  the 
acoustic  velocity  vector,  and  a  single  gradient  of  velocity. 


To  clarify  these  differences,  the  response  of  a  directional  sensor  to  the  flexural 
vibration  of  a  fluid-loaded  plate  (Fig.  3)  was  determined.  Acoustic  and  non-acoustic 
flexural  waves  are  generated  from  the  plate’s  vibration.  In  Fig.  4,  the  theoretical 
wavenumber  response  of  an  infinite  plate,  at  a  frequency  of  1000  Hz  is  given.  Clearly, 
near  the  surface  of  the  plate,  the  predominate  source  of  noise  is  centered  about  the  plate’s 
(flexural)  wavenumber,  denoted  kb.  The  magnitude  of  the  flexural  wavenumber,  for  the 
frequencies  and  plate  geometries  of  interest  here,  is  greater  than  that  of  the  acoustic 
wavenumbers,  and  hence  subsonic.  These  are  evanescent  waves  that  do  not  propagate, 
but  decay  exponentially  away  from  the  surface  of  the  plate.  In  submarine  sonar 
applications,  subsonic  noise  sources  are  common. 
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FIG.  3.  Fluid-loaded  plate  and  directional  sensor  configuration. 
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FIG.  4.  Normalized  magnitude  response  at  fluid-loaded  plate  at  1000  Hz. 


Three  SNRs  have  been  calculated,  each  corresponding  to  measurements  of  pressure, 
velocity,  and  the  gradient  of  velocity.  Signal  is  defined  as  a  planewave,  of  amplitude  P0, 
arriving  along  the  boresight  of  the  z-axis.  Noise  is  defined  to  be  that  due  solely  to 
flexural  vibrations  of  the  plate,  w(kb).  These  ratios  are  written  as 
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The  SNRs,  SNRV  and  SNRV',  at  standoff  distances  of  z0  =  0. 1  m  and  z0  =  0.3  m,  are 
shown  in  Figs.  5  and  6,  respectively  (where,  for  comparison,  Pa  =  F0).  Notice  that  each 
SNR  has  equivalent  exponential  noise  decay;  the  SNRs  improve  quickly  with  standoff. 
Near  the  plate,  however,  the  velocity  gradient  sensor  is  much  more  sensitive  to  flexural 
noise,  particularly  at  low  frequency.  Hence,  for  submarine  sonar  applications,  additional 
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efforts  will  be  required  to  mitigate  against  sources  of  subsonic  self-noise  with  directional 
sensors.  Viscoelastic  sensor  coatings  (passive  wavevector  filters),  standoff,  and  active 
filtering  techniques  may  be  necessary. 


FIG.  5.  SNR  for  a  velocity  and  velocity-gradient 
sensor  4  in.  above  a  steel  plate. 


FIG.  6.  SNR  for  a  velocity  and  velocity-gradient  sensor 
12  in.  above  a  steel  plate. 
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IV.  CONCLUSION 


Underwater  acoustic  intensity  measurements,  along  with  the  requisite  development 
of  directional  sensors,  has  been  an  area  of  active  research  in  Russia  since  the  early  1970s. 
Many  literature  references,  including  textbooks,  are  available.  Shchurov  ’  has  published 
over  50  academic  papers,  describing  many  aspects  of  the  use  of  intensity  in  underwater 
acoustics. 

In  the  West,  less  has  been  published  in  the  area  of  underwater  acoustic  intensity 
measurements  and  the  design  of  directional  sensors.  It  is  a  given  that  directional  sensors 
can  be  more  susceptible  to  certain  noise  mechanisms.  For  example,  at  600  Hz,  the  SNR 
of  the  velocity  gradient  component  (Fig.  5)  of  a  directional  sensor  is  -5  dB  below  that  of 
a  conventional  pressure  sensor.  Below  600  Hz,  the  deficit  is  even  greater.  Reducing  a 
directional  sensor’s  response  to  subsonic  self-noise  will  remain  a  challenging  task. 

The  gains  possible,  though,  with  this  technology  are  significant.  The  point 
directional  sensors  described  here  have  directivity  gains  of  as  much  as  9.5  dB.  In 
addition,  multiple  nulls  may  be  generated  within  a  directional  sensor’s  angular  response. 
These  nulls  can  then  be  used  to  reduce  signal  interference,  acoustic  and  nonacoustic 
sensor  self-noise,  and  scattering  and  diffraction  effects. 

One  specific  goal  for  developing  directional  sensor  technology  is  to  improve  the 
detection  and  localization  capabilities  of  the  AN/BSY-1  low-frequency  hull  array 
(LFHA)  in  SSN  6881  class  submarines  by  replacing  the  array’s  pressure  sensors  with 
directional  receivers. 
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Abstract.  We  consider  a  single  acoustic  sensor  located  at  some  measurement  point 
r0  =  (x0,y0,z0).  A  Taylor  series  for  the  scalar  acoustic  pressure  field  p(t,r )  about  this 

point  would  include  the  scalar  (tensor  of  order  zero)  pressure  p{t,r0  )  as  the  zero-order 
term,  the  pressure  gradient/vector  (tensor  of  order  one)  Vp(t,r0 )  at  the  point  as  part  of  the 
first-order  term,  the  dyadic  (tensor  of  order  two)  Wp(t,r0 )  at  the  point  as  part  of  the 

second  order  term,  and  so  on.  Using  this  Taylor  series,  we  define  a  general  class  of 
directional  acoustic  sensors  as  follows.  A  scalar  acoustic  pressure  sensor  (e.g.,  a 
hydrophone)  will  be  referred  to  as  a  directional  acoustic  sensor  of  order  zero.  This  sensor 
only  measures  the  scalar  acoustic  pressure  p(t,r0 )  at  the  point  r0 ;  its  Taylor  series  about 

r0  assumes  that  the  acoustic  pressure  field  p(t,  r)  about  that  point  is  independent  of  the 
field  point  r  —  (x,  y,  z).  A  vector  acoustic  pressure  sensor  will  be  referred  to  as  a 
directional  acoustic  sensor  of  order  one.  This  sensor  measures  both  p{t,r0 )  and  the 
pressure  gradient/vector  Vp(t,r0 )  at  the  point  r0 ;  its  Taylor  series  about  r0  assumes  that 
the  acoustic  pressure  field  p{t,r )  about  that  point  is  a  linear  function  of  the  field  point  r  . 
Similarly,  a  dyadic  acoustic  pressure  sensor  will  be  referred  to  as  a  directional  acoustic 
sensor  of  order  two.  This  sensor  measures  p(t,r0),  Vp(t,r0 )  and  the  dyadic  VVp(t,r0)  at 

the  point  r0 ;  its  Taylor  series  about  r0  assumes  that  the  acoustic  pressure  field  p(t,  r) 

about  that  point  is  a  quadratic  function  of  the  field  point  r  .  It  will  be  shown  that  the 
higher  the  order  the  directional  acoustic  sensor,  the  more  directional  the  sensor.  We  also 
consider  several  signal  processing  methodologies  for  processing  the  outputs  of  this 
general  class  of  directional  acoustic  sensors.  One  approach  is  a  physics-based  approach, 

which  makes  use  of  the  acoustic  intensity  vector  I  ( t,r )  =  -pit,  r)v(t,r )  ,  where  v(t,  r)  is 
the  acoustic  particle  velocity  vector.  In  addition  to  the  acoustic  intensity  vector,  we 

dn(t  r) 

introduce  the  useful  physics-based  field  intensity  vector  J{t,r )  = - —a(t,r),  where 

dt 

a(t,r )  is  the  acoustic  particle  acceleration  vector.  Another  signal  processing  strategy 
that  will  be  discussed  is  a  multi-channel  filter  approach,  where  the  outputs  of  the 
directional  acoustic  sensor  are  considered  as  the  inputs  to  a  multi-channel  filter,  designed 
for  some  desired  sensor  response.  It  will  be  shown  that  the  element  beam  pattern  of  a 
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directional  acoustic  sensor  can  be  shaped  or  optimized  for  array  gain  by  the  use  of  multi¬ 
channel  filtering.  We  close  the  paper  with  a  discussion  of  how  to  perform  acoustic  signal 
processing  on  an  array  of  directional  acoustic  sensors  and  provide  some  experimental 
results  using  real  data. 

L  INTRODUCTION  AND  HISTORICAL  BACKGROUND  OF 
DIRECTIONAL  ACOUSTIC  SENSORS 

A  single  acoustic  sensor,  located  at  some  measurement  point  r0 ,  that  can  accurately 
detect  and  estimate  the  direction  of  arrival  (DOA)  from  an  arbitrary  (narrowband  or 
broadband)  acoustic  source,  has  always  been  of  interest  to  the  underwater  [1],  [2]  and  in¬ 
air  [3]  acoustic  communities.  For  example,  about  thirty  years  ago,  the  United  States 
Navy  developed  a  single  acoustic  sensor  for  use  in  their  SSQ-53  series  of  Directional 
Low-Frequency  Analysis  and  Recording  (DIFAR)  sonobuoys  [4].  This  DIFAR  sensor  is 
composed  of  a  single  hydrophone  and  two  orthogonal  x-y  accelerometers,  all  co-located 
in  a  single  sensor  package.  The  acoustic  data  out  of  the  DIFAR  sensor  is  basically  a 
multi-channel  time  series,  derived  from  the  simultaneous  measurement  of  the  scalar 
acoustic  pressure  field  at  the  measurement  point  r0 ,  namely  p(t,r0),  and  the  x-y 
components  of  the  acoustic  particle  acceleration  vector  a(t,r0).  By  performing  a 
temporal  low-frequency  Fourier  analysis  on  the  output  of  the  scalar  acoustic  pressure 
sensor  (hydrophone)  or  the  output  of  a  cardioid  beamformer  formed  by  a  linear 
combination  of  the  pressure  channel  and  the  x-y  accelerometer  channels,  one  could  detect 
the  presence  of  an  acoustic  source.  By  finding  the  unit  vector  associated  with  a(t,r0)  and 
using  the  hydrophone  as  a  phase  reference,  one  could  also  unambiguously  find  the  DOA 
of  the  source.  In  the  sonobuoy  community,  this  is  called  DIFAR  signal  processing  and  is 
generally  referred  to  as  SINE  ( sin  <9  beam  pattern  for  the  y-axis  accelerometer),  COSINE 
( cos@  beam  pattern  for  the  x-axis  accelerometer)  and  OMNI  (omnidirectional  beam 
pattern  for  the  hydrophone)  processing.  Based  on  our  brief  definition  in  the  abstract,  the 
DIFAR  sonobuoy  sensor  is  basically  a  two-dimensional  (2-D)  vector  sensor  or  2-D 
directional  acoustic  sensor  of  order  one.  It  is  important  to  note  that  the  above  DIFAR 
sensor  measures  the  exact  scalar  acoustic  pressure  p(t,  r0)  and  the  exact  values  of  the  x 
and  y  components  of  the  pressure  gradient  Vp(t,  r0 ) .  The  latter  follows  from  the  fact  that 

Vp(t,r  )=-p(r)a(t,r)  (1) 

for  linear  acoustic  waves  propagating  in  an  inhomogeneous  acoustic  medium.  Here, 
p(r)  is  the  mass  density  of  the  medium  at  the  field  point  r  .  Thus,  in  order  to  know  the 
exact  value  of  the  x  and  y  components  of  the  pressure  gradient  at  the  measurement  point 
r0 ,  we  assume  that  the  DIFAR  sensor  exactly  measures  the  x  and  y  components  of 
a(t,r0)  and  knows  p(r)  at  the  point  r0 .  The  exact  value  of  Vp(t,r0)  follows  from  (1). 

About  the  same  time  the  U.S.  Navy  developed  the  DIFAR  sonobuoy  sensor,  they 
also  developed  another  2-D  directional  acoustic  sensor  of  order  one  called  the  multimode 
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hydrophone,  which  is  used  as  a  directional  sensor  in  the  AN/WLR-9  series  of  acoustic 
intercept  receivers  [5].  This  2-D  vector  sensor  can  be  modeled  by  two  orthogonal  pressure 
dipoles  in  the  x-y  plane.  Figure  1  gives  a  pictorial  representation  of  the  x  and  y-axis 
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FIG.  1.  The  multimode  hydrophone  can  be  modeled  by  two 
orthogonal  pressure  dipoles  in  the  x-y  plane. 


pressure  dipoles.  The  x-axis  pressure  dipole  consists  of  two  scalar  pressure  sensors 
separated  by  a  distance  Ax  along  the  x-axis.  Similarly,  the  y-axis  dipole  consists  of  two 
scalar  pressure  sensors  separated  by  a  distance  Ay  along  the  y-axis.  Unlike  the  DIFAR 
sensor,  the  multimode  sensor  cannot  measure  the  exact  values  of  the  x  and  y  components 
of  the  pressure  gradient  Vp(t,r0),  since  it  uses  four  scalar  pressure  sensors  instead  of  two 
orthogonal  accelerometers.  However,  we  can  approximate  these  pressure  gradient 
components  at  the  measurement  point  r0  =  (0,0,0)  by  using  the  finite  difference 
approximations: 
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Similarly,  we  can  approximate  the  scalar  acoustic  pressure  at  the  point  r0  -  (0,0,0)  by 
averaging  the  four  pressure  sensors;  that  is, 
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Using  Eqs.  (1),  (2a),  and  (2b),  we  can  also  approximate  the  x-y  components  of  the 
acoustic  particle  acceleration  vector  a(r,0,0,0),  provided  we  know  the  mass  density  p(r) 
at  the  point  r0  =  (0,0,0) .  Referring  to  the  previous  DIFAR  discussion,  we  see  that  the 
multimode  acoustic  sensor  data  is  approximately  in  the  same  multi-channel  time  series 
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form  as  the  DIFAR  acoustic  sensor  data.  Thus,  in  principle,  we  can  use  the  SINE, 
COSINE  and  OMNI  DIFAR  signal  processing  techniques  to  process  the  outputs  of  the 
multimode  hydrophone.  The  acoustic  intercept  community  refers  to  (2b)  as  the  SINE 
channel  (approximate  sin<9  beam  pattern  for  the  y-axis  pressure  dipole),  (2a)  as  the 
COSINE  channel  (approximate  cos 9  beam  pattern  for  the  x-axis  pressure  dipole)  and  (3) 
as  the  OMNI  channel  (omnidirectional  beam  pattern  for  (3))  and  uses  DIFAR-like 
processing  to  detect  and  estimate  the  DOA  of  acoustic  emissions. 

For  about  30  years,  researchers  in  the  field  of  in-air  acoustics  have  used 
directional  acoustic  sensors  of  order  one  (vector  sensors)  to  estimate  the  magnitude  and 
direction  of  the  acoustic  intensity  vector  in  order  to  obtain  the  acoustic  power  per  unit 
area  (Watts  /  to2)  radiated  by  a  source  as  a  function  of  time  t  and  position  r  . 

Specifically,  researchers  at  General  Motors  Corporation  have  used  the  instantaneous 
acoustic  intensity  vector 

l{t,r  =-p(t,r))v(t,r)  (4) 


to  measure  the  in-air  acoustic  power  per  unit  area  generated  by  various  vehicles  [6]. 

Here,  v(t,r )  is  the  acoustic  particle  velocity  vector,  which  can  be  obtained  by  performing 
a  time  integration  of  -Vp(t,r)/ p(r) .  It  is  important  to  note  that  (4)  is  valid  for  any  type 
of  linear  acoustic  wave  propagating  in  an  inhomogeneous  acoustic  medium;  it  does  not 
have  to  be  a  plane  wave.  The  magnitude  of  (4)  gives  the  instantaneous  acoustic  power 
per  unit  area  and  the  unit  vector 


ii(t,r)  = 


1 7 


(5) 


gives  the  instantaneous  DOA  of  the  source’s  radiated  acoustic  power. 

Nehorai  [7]  and  Shchurov  [8]  have  done  an  extensive  amount  of  research  on  the 
signal  processing  associated  with  acoustic  vector  sensors.  For  a  single  vector  sensor, 
both  of  these  researchers  use  the  acoustic  intensity  vector  (4)  as  the  basis  for  their 
processing  methodologies.  Their  basic  approach  is  to  perform  a  time  average  of  (4)  and 
use  this  time-averaged  intensity  vector  to  derive  the  corresponding  unit  vector,  which 
results  in  the  DOA  estimate  of  the  acoustic  source.  In  this  paper,  we  shall  refer  to  this 
approach  as  the  acoustic  intensity  algorithm.  Nehorai  [9]  has  also  developed  a 
localization  algorithm  for  estimating  the  position  vector  of  an  acoustic  source  (range, 
azimuth  angle  and  elevation  angle  of  the  source)  using  an  arbitrary  distributed  array  (e.g., 
the  array  does  not  have  to  be  linear)  of  acoustic  vector  sensors.  His  approach  is  to  first 
use  the  above  acoustic  intensity  algorithm  at  each  sensor  location  to  derive  a  unit  vector 
to  the  source  relative  to  a  specific  sensor  in  the  array.  If  the  conditions  were  ideal,  all  the 
various  unit  vector  estimates  from  each  sensor  would  intersect  at  one  point,  giving  the 
position  vector  of  the  source.  In  practice,  when  the  conditions  are  not  ideal,  Nehorai  uses 
a  least-squares  approach  to  combine  all  the  unit  vectors  from  each  sensor  to  give  the 
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‘best’  estimate  of  the  source’s  position  vector.  In  this  paper,  we  shall  refer  to  this 
algorithm  as  the  intensity-based  ranging  algorithm. 

Cray  and  Nuttall  [10]  have  investigated  linear  signal  processing  approaches  for 
beamforming  single  and  multiple  acoustic  vector  sensors  that  measure  the  acoustic 
pressure  p{t,r0 )  and  the  three  acoustic  particle  velocity  vector  components  vx(t,r0), 

vy  (t,  r0 )  and  vz  (t,  r0  ) .  (Note  that  this  is  equivalent  to  measuring  p(t,  r0 )  and  Vp(t,  r0 ) , 
using  (1)  to  get  the  acoustic  particle  acceleration  vector  a(t,r0),  then  performing  a  time 
integration  on  a(t,r0)  to  get  the  acoustic  particle  velocity  vector  v  (t,^)).  Specifically, 

they  show  that  by  forming  a  linear  combination  of  these  vector  sensor  measurements, 
namely, 

b(t,  r0 )  =  w0p(t,  r0)+  wlxvx{t,  F0  )+  wlvvv  ( t ,  F0  )+  wlzvz  (t,  F0  )  (6) 

=  w0p(t,P0)+w1-v(t,P0) 


one  can  choose  the  real  weights  vv0 ,  vv,  =  {wu ,  w1  v ,  w, . )  so  that  this  single  vector  sensor 

behaves  like  a  spatial  filter.  In  this  sense,  a  single  vector  sensor  can  perform  spatial 
filtering  or  beamforming,  with  b(t,r0)  representing  the  output  of  the  beamformer. 

In  this  paper,  we  make  use  of  the  Taylor  series  expansion  of  the  pressure  field 
p{t,r )  about  the  point  r0in  order  to  define  a  general  class  of  directional  acoustic  sensors. 

(Refer  to  the  abstract).  We  then  investigate  sensors  that  approximate  the  pressure 
gradient  by  using  orthogonal  pressure  dipoles  and  finite  difference  approximations  and 
show  that  extreme  caution  must  be  exercised  when  doing  so.  In  order  to  deal  with 
directional  sensors  that  measure  the  pressure  gradient  directly  (e.g.,  directional  sensors 
that  use  accelerometers),  we  introduce  the  useful  acoustic  field  intensity  vector 
dp(t  r) 

j(t<r)  = - —ci(t,r),  which  deals  directly  with  the  pressure  gradient  or  the 

dt 

acceleration  vector.  Also,  because  it  involves  the  time  derivative  of  the  acoustic  pressure, 
this  field  intensity  vector  naturally  provides  for  equalization  of  an  ambient  acoustic 
pressure  background,  like  the  underwater  environment,  that  has  an  acoustic  noise 
spectrum  that  falls  off  at  6  dB  per  octave.  We  also  relate  the  work  of  Cray  and  Nuttall  to 
the  aforementioned  Taylor  series  expansion  and  generalize  their  concept  of  beamforming 
to  the  general  problem  of  shaping  the  element  pattern  of  a  directional  acoustic  sensor.  A 
discussion  of  how  to  beamform  a  practical  array  of  directional  sensors  will  be  given.  In 
the  latter  part  of  the  paper,  we  provide  experimental  results  that  show  the  practical  value 
of  the  acoustic  intensity  vector,  the  usefulness  of  Nehorai’s  ranging  algorithm  and  how 
the  directional  sensor’s  spatial  filtering  ability  could  be  used  in  a  noisy  environment. 

II.  PRESSURE  GRADIENT  ESTIMATION 

In  both  the  underwater  and  in-air  communities,  the  pressure  gradient/vector,  needed 
for  a  directional  acoustic  sensor  of  order  one,  is  found  in  a  variety  of  different  ways.  As 
discussed  above,  one  can  measure  the  acoustic  particle  acceleration  vector  a(t,r0 ) 
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directly  and  use  (1)  and  knowledge  of  the  mass  density  at  the  measurement  point  r0  to 
obtain  a  fairly  accurate  value  of  the  pressure  gradient  Vp(t,r0).  Alternatively,  one  could 
measure  the  acoustic  particle  velocity  vector  v(t,  r0)  directly,  perform  a  time  derivative  to 
get  the  acceleration  vector  and  proceed  as  above  to  get  the  pressure  gradient.  (The 
accuracy  of  this  method  would  depend  on  how  well  one  implemented  the  time 
derivative).  If  one  decides  to  use  scalar  pressure  sensors  to  approximate  the  pressure 
gradient,  then  the  finite  difference  approximations  in  (2)  would  have  to  be  computed. 
However,  great  care  must  be  exercised  when  approximating  the  pressure  gradient  by  (2). 
For  example,  (2a)  approximates  the  x-component  of  the  pressure  gradient  at  the 
measurement  point  r0  —  (0,0,0) .  (Refer  to  Fig.  1).  If  the  pressure  field  was  the  arbitrary 
(narrowband  or  broadband)  plane  wave 


p(t*r)=f 


t  +  ‘ 


n  •  r 


(7) 


which  is  assumed  to  be  propagating  toward  the  origin  of  the  coordinate  system  in  Fig.  2, 
then  the  exact  value  of  the  temporal  Fourier  transform  of  the  pressure  gradient  would  be 


wavefront 


FIG.  2.  The  propagation  of  a  planar  wavefront  towards  the  origin 
of  a  rectangular  coordinate  system. 


VP(®,r)=  jF(oj)kefk  r  =  jP(co,r)k  .  (8) 

Here,  the  wavenumber  or  propagation  vector  k  =kn has  the  magnitude  k  =  0)/c  =  2n / X 
and  the  direction  n ,  where  h  is  the  unit  vector 
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n  =  sin0r  cos0rx  +  sin  (j)T  sin0ry  +  cos0r£. 


(9) 


The  position  vector 

r  =  xx  +  yy  +  zz  =  rsirufrcosdx  +  r  sirup  cosdy  +  rcos(f>z  (10) 

is  the  vector  from  the  origin  to  an  arbitrary  field  point  r  -  (x,  y,  z)  in  rectangular 
coordinates  or  r  -  (rsin0cos0,rsin0sin0,rcos0)  in  spherical  coordinates.  The 
temporal  Fourier  transform  of  the  plane  wave  pressure  field  (7)  is  found  by  using 

P(©,r)  =  J  p(t,r)e~jmdt .  (11) 

Now  if  we  use  (2a)  to  approximate  the  x-component  of  the  pressure  gradient  at  the  origin 
(the  point  r0  -  (0,0,0) ),  then  the  temporal  Fourier  transform  of  this  approximation  can  be 
expressed  as 


dx 


=  Px  .  ?0  )  «  Px  {(O,r0  )  =  jF  (to)- 


(  Ax 

2sin|  n — sin^rcos0r 

A 


Ax 


(12) 


Since  the  exact  value  of  PX((Q,  r0)  is  given  by 


2tt 

Px  (co,  r0 )  =  JF  ( (O ) —  sin  (j)T  cos  0T , 

A 


(13) 


then  the  fractional  error  in  estimating  the  x-component  of  Px(co,r0)  can  be  expressed  as 


1- p /p  =1- 


(  Ax 

sin  K — sin0r  cos0r 
_  1 _ _ _ 

it — sin^T  cos0T 
A  Yt  t 


( 


Ax  . 


\ 


=  l-sinc  K — ^sin^cos^ 

l  A  J 


(14) 


Notice  that  (14)  is  frequency  dependent  and  is  valid  only  when  the  two  hydrophones  in 
the  x-axis  dipole  (Refer  to  Fig.  1)  are  exactly  matched  in  their  amplitude  and  phase 
frequency  responses.  In  practice,  this  is  seldom  the  case,  so  we  must  modify  our  analysis 
to  include  the  effects  of  mismatch  in  the  frequency  responses  of  the  hydrophone  and 
corresponding  signal  conditioning  associated  with  each  hydrophone  in  the  dipole.  It  can 
be  shown  that  when  the  two  hydrophones  of  the  x-axis  dipole  are  mismatched  in  the 
magnitudes  of  their  combined  (sensor  plus  electronics)  frequency  responses,  the 
magnitude  of  the  fractional  error  in  estimating  the  x-component  of  Px((D,r0)  can  be 
expressed  as  [11] 
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1  -PJPr 


=  1- 


82  +  4(l  +  c>)sin2 

(Ax  Y 

K — sin<j)T  cos0r 

L 

Ax 

2  71  —  sin(j)T  cos0r 
X 

(15) 


which  is  valid  for  sin  </>,  cos  0T  =£  0 .  Here,  8  is  a  frequency-dependent  mismatch 
parameter  defined  by 

c>(®)=  [5,1(m)G1(tt))-52(tt))G2(ft))]/52(ft))G2(tt))  ,  (16) 

where  5j(®)  and  G,  (ft))  are  the  magnitudes  of  the  frequency  responses  of  the  hydrophone 
and  signal  conditioning,  respectively,  for  the  hydrophone  located  on  the  positive  x-axis 
(refer  to  Fig.  1)  and  S2(co)  and  G2(co)  are  the  magnitudes  of  the  frequency  responses  of 
the  hydrophone  and  signal  conditioning,  respectively,  for  the  hydrophone  located  on  the 
negative  x-axis.  The  contour  plot  in  Fig.  3  gives  the  magnitude  of  the  fractional  error  in 
estimating  the  x-component  of  the  pressure  gradient  or  the  x-component  of  the  acoustic 
particle  acceleration  vector  as  a  function  of  both  8  and  A x/X  under  the  condition 
dT  =  0° ,  (j)T  =  90° .  [Note  that  this  condition  is  for  a  plane  wave  traveling  in  the  negative 
x-direction  in  the  x-y  plane,  which  means  that  the  wave  hits  the  main  response  axis  of  the 
dipole  or  approximate  accelerometer]. 


Fi««t«w  EweMMApHucMei  a«cmio**i*i  n**oo*s* 
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FIG.  3.  Contour  plot  of  fractional  error  (magnitude)  as  a  function  of  8  and  Ax/A,. 
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The  mismatch  parameter  S  is  plotted  on  the  y-axis.  Ax/ X  is  plotted  on  the  x-axis 
and  each  contour  represents  1.2%  of  fractional  error  magnitude.  For  example,  suppose 
we  wanted  to  estimate  the  x-component  of  the  pressure  gradient  of  an  acoustic  field  at 
100  Hz  using  a  pressure  hydrophone  dipole,  with  spacing  Ax  =3  inches,  with  a  fractional 
error  magnitude  of  10%.  From  Fig.  3,  this  implies  that  the  amplitude  sensitivities  of 
these  hydrophones  must  be  matched  better  than  0.01  dB.  This  is  very  difficult  to  achieve 
in  practice.  If  the  same  dipole  was  used  to  estimate  the  pressure  gradient  at  1000  Hz, 
2000  Hz,  and  10,000  Hz,  then  the  amplitude  sensitivities  of  these  hydrophones  must  be 
matched  better  than  approximately  0.03  dB,  0.04  dB,  and  0.8  dB,  respectively.  Thus,  if 
one  had  a  mechanical  packaging  constraint  to  build  a  low-frequency,  directional  sensor  in 
a  small,  3-inch,  diameter  x  6-inch  length  right  circular  cylinder,  one  might  consider  using 
accelerometers  or  velocity  sensors  vice  orthogonal  pressure  dipoles  to  estimate  the 
pressure  gradient. 

III.  TAYLOR  SERIES  AND  DIRECTIONAL  ACOUSTIC  SENSORS 

D’ Spain  [12]  was  the  first  to  point  out  a  qualitative  relationship  between  the 
Taylor  series  expansion  of  the  pressure  field  and  a  vector  sensor.  Following  his 
observation,  we  consider  a  single  acoustic  sensor  located  at  some  measurement  point 
r0  -  (x0,y0,z0).  A  Taylor  series  for  the  scalar  acoustic  pressure  field  p(t,r )  about  this 
point  can  be  written  as 


p(*>  ?)  =  pit,  r0 ) ■ +  (r  -  r0 )  ■ ■  Vp(t,  r0 ) ■ + 1  (r  -  r0 ) ■ •  VVp(t,  r0 ) ■  (r  -  r0 ) ■ +  h.o.t.  (17) 

The  interpretation  of  (17)  is  as  follows.  A  directional  acoustic  sensor  is  capable  of 
extrapolating  the  acoustic  pressure  field  beyond  the  measurement  point  r0  so  that  it 
actually  knows  this  field  at  every  point  inside  a  sphere  of  radius  R  =  \r  -  r0\ ,  where  R  is 
defined  by  the  type  of  directional  sensor  and  the  error  associated  with  the  extrapolation. 
The  Taylor  series  (17)  is  the  formula  used  to  do  the  wavefield  extrapolation  with  some 
specified  error  e(t,r).  For  example,  a  directional  sensor  of  order  zero,  can  only  measure 
the  pressure  at  r0 ,  so  its  estimate  of  the  field  beyond  this  point  is  p{t,r)  =  p(t,  r0 )  and  the 
corresponding  estimation  or  extrapolation  error  is  e(t,  r)  =  p(t,r)—  p(t,r0).  If  the  error  is 
required  to  be  small  (less  than  10%),  then  R  will  be  small.  This  implies  that  the  aperture 
(2R)  of  a  single  pressure  sensor  is  small,  so  by  the  theory  of  spatial  Fourier  transforms 
[13],  the  zero-order  sensor  is  essentially  omnidirectional.  However,  a  sensor  of  order  one 
(vector  sensor)  measures  both  the  pressure  and  pressure  gradient  at  r0 ,  so  its  estimate  of 
the  field  beyond  this  point  is  p(t,r  )  =  p(t,  r0 )+  (r  -  r0 )•  Vp(t, r0 )  and  the  corresponding 
estimation  error  is 

e(t,r)=  p(t,r)-p(t,r )  (18) 

e(t,  r  )  =  pit,  r ) ■ -  [p{t, r0 )■ +  (r  -  r0 ) ■ ■  Vp(f ,  r0 )] 
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For  the  same  error,  the  vector  sensor  will  have  a  larger  aperture  than  the  scalar  sensor,  so 
it  should  be  more  directional.  Thus,  the  dyadic  sensor  should  be  better  than  the  scalar  and 
vector  sensors  at  estimating  the  field,  so  it  should  be  even  more  directional. 

We  define  the  mean-squared  error  in  estimating  the  pressure  field  beyond  r0  by 

1  K  T 

MSE  =  —  ffl e(t,rfdtdp  ,  (19) 

^00 

where  T  is  a  suitable  integration  time  (e.g.,  the  temporal  period  of  a  plane  wave)  and  j5  is 
the  angle  between  h  and  f .  Note  that  r  is  a  unit  vector  in  the  direction  of  r  .  (Refer  to 
Fig.  2).  Dividing  (19)  by  the  average  power  in  the  acoustic  pressure  field  we  obtain  the 
normalized  MSE.  Figure  4  shows  the  normalized  mean-squared  as  a  function  of  R/  A 
for  the  scalar,  vector  and  dyadic  sensors.  For  a  specified  normalized  MSE  of  10%,  notice 
that  the  acoustic  aperture  (2R)  of  the  scalar  sensor  is  about  A/ 10,  whereas  the  apertures 
for  the  vector  and  dyadic  sensors  are  A/3  and  A/2,  respectively. 


FIG.  4.  Normalized  mean-squared  estimation  error  (Eq.  (19))  vs.  R/A. 
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IV.  MULTICHANNEL  FILTERING  AND  DIRECTIONAL 
ACOUSTIC  SENSORS 

For  the  purposes  of  this  section,  it  is  convenient  to  deal  with  the  temporal  Fourier 
transform  of  the  Taylor  series  (17).  Rewriting  (17)  in  the  frequency  domain  we  obtain 

p((0,  r  )  =  W0 (co,  r  )P{(0,  r0 ) ■ +  W{ (co,  r  ) •  VP(co,  r0)+W2((O,r)- V VP(co,  r0)-W2(co,r),  (20) 

where  we  have  retained  only  the  first  three  terms  of  the  series  (dyadic  sensor)  and  the 
weights  W0 ,  Wt  and  W2  are  chosen  such  that  the  right  side  of  (20)  is  consistent  with  the 
Taylor  series  approximation  of  the  pressure  field.  Specifically,  JV0  =  1 ,  W,  =  r  -  r0  and 
W2  =  (r  -  r0 )/  V2  are  the  only  weights  that  allow  (20)  to  be  a  Taylor  series  extrapolation 
of  the  pressure  field  by  a  dyadic  sensor.  If  we  pick  another  set  of  weights,  we  cannot  be 
assured  that  (20)  will  be  a  valid  Taylor  series.  However,  this  does  raise  the  question  that 
maybe  we  could  select  a  different  set  of  weights  that  might  be  useful  for  another  purpose, 
like  shaping  the  beam  pattern  of  a  directional  sensor  or  finding  a  set  of  weights  that 
results  in  the  maximum  array  gain  in  isotropic  noise.  That  is,  the  linear  combination  on 
the  right  side  of  (20)  could  be  expressed  as 

B(co,r)  =  W0 (co,  r  )p{co, r^+W^CO,?)-  VP(co,  r0)+W2((O,r)- VVP(co,  r0)-W2 (ft),  r )  (21) 

where  the  function  B{co,r )  could  represent  the  output  of  a  frequency  domain 
multichannel  filter  with  P(oj,  r() ) ,  V Z3 (co ,  7; , )  and  WP(co,r0 )  as  the  inputs.  Refer  to  Fig.  5 
for  a  pictorial  representation  of  this  concept. 


p(uro) 

- ► 

Vp{t,rt ■  )  ^ 

Multi-channel  filter 

- ►  b(t,r) 

VV/?(/,r0)^ 

FIG.  5.  Multichannel  Altering  and  directional  acoustic  sensors. 


Thus,  (21)  is  a  generalization  of  (20)  and  can  be  viewed  as  Multichannel  filtering 
approach  to  directional  sensors. 

Following  Silvia  and  Schmidlin,  we  could  select  the  above  weights  in  (21)  to  steer 
the  beam  of  a  vector  sensor  so  that  we  have  maximum  sensitivity  in  the  look  direction 
and  place  a  null  in  another  direction.  Further,  we  could  select  the  weights  so  that  the 
directional  sensor  produces  a  beam  pattern  that  results  in  maximum  array  gain  [14]. 
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Figure  6  shows  the  two  weight- selection  criteria  for  the  vector  sensor,  and  Fig.  7 
does  the  same  for  the  case  of  the  dyadic  sensor. 


105° 


FIG.  6.  Selection  of  weights  for  a  vector  sensor  for:  (a)  maximum  array  gain, 

and  (b)  optimum  null  placement. 


For  maximum 
Array  gala 


For  optimum 
Null  placemeat 


270 


Array  Gala  "7dB 

(b) 


FIG.  7.  Selection  of  weights  for  a  dyadic  sensor  for:  (a)  maximum  array  gain, 

and  (b)  optimum  null  placement. 
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V.  EXPERIMENTAL  RESULTS 


The  experimental  results  for  the  acoustic  intensity  algorithm,  intensity-based 
ranging  algorithm,  and  beam  steering  algorithms  are  discussed  in  the  oral  presentation. 
Refer  to  the  viewgraph  material. 

VI.  CONCLUSIONS 

We  have  shown  that  the  Taylor  series  expansion  of  the  pressure  field  provides  an 
elegant  mathematical  framework  for  the  study  of  directional  acoustic  sensors  of  all 
orders.  The  acoustic  intensity  algorithm  generalizes  the  DIFAR  and  acoustic  intercept 
receiver  SINE,  COSINE  and  OMNI  signal  processing.  Pressure  gradient  estimation 
using  orthogonal  pressure  dipoles  must  be  done  with  caution.  Multichannel  filtering  of 
the  directional  sensor  outputs  provides  a  useful  tool  for  obtaining  a  desired  sensor 
response.  For  example,  the  weights  of  the  filter  can  be  chosen  to  give  a  vector  and  dyadic 
sensor  a  maximum  array  gain  relative  to  a  pressure  sensor  of  6  dB  and  9  dB,  respectively. 
The  higher  the  order  the  directional  sensor,  the  more  directive  the  sensor,  as  seen  by  the 
fact  that  a  vector  sensor  has  a  minimum  beamwidth  of  105  degrees  whereas  a  dyadic 
sensor  has  a  minimum  beamwidth  of  65  degrees.  Steering  a  vector  sensor  and  placing 
nulls  in  the  element  response  must  also  be  done  with  caution.  In  summary,  the  proposed 
theoretical  structure  allows  one  to  achieve  up  to  9.5  dB  array  gain  by  measuring  the 
acoustic  field  properties,  including  scalar,  vector  and  dyadic  at  a  single  point. 
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Abstract:  An  underwater  acoustic  intensity  probe  containing  a  hydrophone  and 

accelerometer  is  described.  The  design  process  starts  from  first  principles  in  which  the 
dynamics  of  an  accelerometer  that  is  imbedded  in  a  compliantly  suspended  sphere  are 
evaluated  using  theory  and  experiment.  The  results  are  extended  to  the  case  of  the 
intensity  probe,  which  has  the  geometry  of  a  right  circular  cylinder.  Lumped  parameter 
circuit  analysis  is  used  to  assess  the  impact  that  external  suspension  systems  and  internal 
sensor  dynamics  have  on  the  fidelity  of  the  measurement.  Experimental  data  is  presented 
which  shows  the  performance  of  the  sensor  in  a  standing  wave  field.  Limitations  of  the 
sensing  technology  are  discussed. 

INTRODUCTION 

The  motivation  for  this  research  concerns  the  need  to  develop  an  underwater  acoustic  intensity  probe 
that  exhibits  a  bandwidth  that  extends  well  into  the  kilohertz  region  of  the  audible  sound  spectrum. 
Neutrally  buoyant  pressure-velocity  (p-u )  intensity  probes1,2  are  very  novel  transducers,  but  are  somewhat 
limited  to  frequencies  below  2  kHz.  The  performance  of  the  geophone  is  a  chief  limitation  that  sets  the 
operational  bandwidth  of  such  probes.  Accordingly,  the  development  of  a  pressure-accelerometer  (p-a ) 
intensity  probe  has  been  pursued,  although  it  should  be  stated  that  the  concept  is  not  new.  Sykes3  and 
Schloss4  patented  such  devices  in  the  mid- 1960’s,  however,  each  design  is  relatively  primitive  with  respect 
to  today’s  technological  advances.  Nevertheless,  the  p-a  probe  described  in  this  paper  is  intended  for  use 
as  a  diagnostic  measurement  tool  to  aid  in  the  optimization  of  underwater  sound  projectors  that  are  utilized 
by  air-deployed  sonobuoys.5  Other  applications  involve  the  combination  of  three  orthogonal  inertial 
sensors  (e.g.,  either  velocity  or  acceleration  sensitive  devices)  in  conjunction  with  a  co-located  pressure 
sensor.  Such  sensors  are  termed  ‘acoustic  vector  sensors’  and  have  application  to  a  variety  of  tactical 
acoustic  sensing  platforms  that  are  operated  by  the  U.  S.  Navy. 

I.  SPHERICAL  VELOCITY/ACCELERATION  SENSORS 

The  fidelity  of  underwater  acoustic  intensity  measurements  is  largely  dependent  on  the  ability  to 
determine  the  particle  velocity  with  a  high  degree  of  precision.  This  is  true  for  p-u ,  p-a ,  p-p,6  and  m-m7  8 
intensity  probes.  To  illustrate  the  basic  physics  of  the  problem,  a  spherical  sensor  containing  an 
accelerometer  is  evaluated  using  theory  and  experiment.  The  object  of  this  analysis  is  to  quantify  the  effect 
of  the  sphere  density,  external  suspension  systems,  and  internal  sensor  dynamics  on  the  fidelity  of  the 
measurement.  The  results  are  used  to  aid  in  the  development  of  the  p-a  probe. 

A.  Theory 

Many  researchers9'13  have  derived  the  equation  of  motion  for  a  small,  unconstrained  solid  sphere  when 
subjected  to  an  acoustic  plane  wave  in  an  unbounded  fluid  medium.  The  analysis  shows  that  the  velocity 
amplitude  of  the  sphere  is  related  to  that  of  the  acoustic  wave  by: 


Vs  =  mo  +mi 
V0  ms+mi  ’ 


(1) 
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where  Vs  is  the  velocity  amplitude  of  the  sphere,  V0  is  the  velocity  amplitude  of  the  acoustic  wave,  ms  is  the 
mass  of  the  sphere,  is  the  induced  mass  that  is  created  as  a  result  of  the  sphere  translating  in  the  fluid, 
and  mQ  is  the  mass  of  fluid  having  the  same  volume  as  the  sphere.  For  the  case  of  the  sphere,  the  induced 
mass  can  be  computed  analytically  in  closed  form  and  is  equal  to  mJ2.  This  leads  to  well  known 
formulation: 


Vs  _  3p0 

K  2  ps+po 


where  p0  is  the  density  of  the  fluid  and  ps  is  the  density  of  the  sphere.  It  can  be  concluded  from  Eqs.  (1) 
and  (2)  that  a  neutrally  buoyant  sphere  has  the  same  velocity  amplitude  as  the  incoming  acoustic  wave,  a 
positively  buoyant  sphere  responds  with  an  amplitude  that  is  greater  than  that  of  the  acoustic  wave,  and  a 
negatively  buoyant  sphere  responds  with  an  amplitude  that  is  less  than  that  of  the  acoustic  wave.  These 
results  serve  as  the  basis  for  nearly  every  underwater  transducer  that  employs  an  inertial  sensor  to  measure 
the  acoustic  particle  velocity,  or  time  derivatives  thereof,  e.g.,  particle  acceleration. 

The  dynamics  of  the  problem  are  somewhat  different  when  the  sphere  is  constrained  by  a  compliant 
suspension  system.  Fig.  1  shows  a  schematic  of  a  suspended  sphere  along  with  a  lumped  parameter  circuit 
model  of  the  system.  In  Fig.  1,  the  spring  has  compliance  Cs  and  the  damper  has  resistance  Rs.  The 
equation  of  motion  can  be  inferred  from  the  circuit  and  ultimately  leads  to  the  following  formulation: 


Vo 


3  Po 

2 Ps  +  Po 


1 


®v 

CO 


(3) 


where  co  is  the  excitation  frequency,  cos  is  the  resonance  frequency,  and  £y  is  the  damping  factor.  The 
resonance  frequency  and  damping  factor  are  further  defined  as  cos=[(ms+  m ,•)  Cs]'m  and  =  RJ2(ms+  m^)(Os. 


FIG.  1.  Schematic  of  a  compliantly  suspended  sphere  undergoing  rigid  body  translation  when  subjected  to  an  acoustic  plane  wave. 
The  equivalent  mechanical  impedance  circuit  is  also  shown. 


Fig.  2  presents  the  result  of  plotting  Eq.  (3)  versus  dimensionless  frequency  cola\  over  damping  factors 
that  span  two  orders  of  magnitude.  This  analysis  is  specific  to  the  case  of  neutral  buoyancy.  It  can  be 
inferred  from  the  figure  that  about  a  decade  above  resonance,  the  magnitude  of  the  sphere’s  motion  is 
identical  to  that  of  the  acoustic  wave.  The  phase  response,  however,  depends  on  the  level  of  damping 
present  in  the  system,  and  in  particular,  only  lightly  damped  systems  result  in  good  phase  fidelity  above 
resonance  (e.g.,  frequencies  where  the  sphere  exhibits  little  or  no  phase  shift  with  respect  to  that  of  the 
acoustic  wave).  These  results  indicate  that  the  motion  of  a  lightly  damped  sphere  can  be  described  by  Eqs. 
(1)  and  (2)  at  about  a  decade  above  resonance. 

The  dynamics  of  the  system  are  further  exacerbated  when  an  inertial  sensor  such  as  a  piezoceramic 
accelerometer  or  moving  coil  geophone  is  imbedded  in  the  sphere.  In  either  case  the  system  can  be 
modeled  with  lumped  parameters  to  compute  the  open-circuit  receiving  response  of  the  transducer.  Figs.  3 
and  4  present  schematics  of  acceleration-based  and  velocity-based  spherical  transducers  along  with  their 
corresponding  lumped  parameter  circuit  models.  As  before,  the  spheres  contain  an  external  suspension 
system  and  are  subjected  to  an  acoustic  plane  wave  in  an  unbounded  medium. 
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FIG.  2.  Response  of  a  compliantly  suspended  sphere  undergoing  rigid  body  translation  when  subjected  to  an  acoustic  plane  wave. 

In  Fig.  3,  mt ,  Ct,  and  Rt  are  the  mechanical  mass,  compliance,  and  damping  of  the  accelerometer.  The 
transduction  constant  is  denoted  as  (/)  and  C0  is  the  clamped  capacitance.  The  parameters  A0,  A^and  At, 
correspond  to  the  acceleration  amplitudes  of  the  in-coming  wave,  sphere,  and  proof-mass,  respectively. 
Similarly,  in  Fig.  4,  mt ,  Ct,  and  Rt  are  the  mechanical  mass,  compliance,  and  damping  of  the  geophone.  The 
transduction  constant  is  denoted  as  Bl ,  Re  is  the  coil  resistance,  and  Le  is  the  coil  inductance.  The  parameter 
Vt  corresponds  to  the  velocity  amplitude  of  the  geophone  proof-mass.  The  open  circuit  output  voltage  for 
each  transducer  is  denoted  as  eoc. 


ms+mi  Cs  Rs  mt 


FIG.  3.  Schematic  of  an  acceleration-based  spherical  sensor  with  external  suspension  system.  Sphere  is  undergoing  rigid  body 
translation  as  a  result  of  acoustic  plane  wave  excitation.  The  equivalent  electro-mechanical  impedance  circuit  is  also  shown. 


ms+mi  Cs  Rs  mt 


FIG.  4.  Schematic  of  a  velocity-based  spherical  sensor  with  external  suspension  system.  Sphere  is  undergoing  rigid  body  translation 
as  a  result  of  acoustic  plane  wave  excitation.  The  equivalent  electro-mechanical  impedance  circuit  is  also  shown. 
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The  open-circuit  receiving  response  of  the  spherical  accelerometer  is  determined  from  Fig.  (3)  to  be 
^oc=0Z(m)/(/m)2Co.  Likewise,  from  Fig.  (4),  the  receiving  response  of  the  spherical  geophone  is  determined 
to  be  eoc=BlZ(co).  For  both  of  these  expressions,  the  term  Z(co)  is  defined  as: 


Z(fi>)  = 


3  Po 


2ps+Po 


,  ^  CO,  CO, 

1-2 

C 0  CO 


2  Y 


A 


,  „  CO,  CO, 

1-2  jC,— — y 

0)  CO2 


171,  +  171; 


-  COj  COf 

2j£(— +-^ 
CO  CO2 


-1 


(4) 


where  tq  and  £  are  the  resonance  frequency  and  damping  factor  associated  with  the  inertial  transducer 
inside  the  sphere.  These  parameters  are  further  defined  as  C0i={mtCtym  and  £ \-Rt/2mto\ .  To  assess  the 
significance  of  these  results,  the  receiving  response  of  both  transducers  is  plotted  versus  frequency  f=  cof2n. 
It  is  assumed  that  each  sphere  is  lightly  damped  and  exhibits  a  low  resonance  frequency  such  that 
f=CDs/27T=5  Hz  and  £,=0.01.  It  is  also  assumed  that  the  resonance  frequencies  of  the  accelerometer  and 
geophone  are  /r=tq/2;r=25  kHz  and  ft=  30  Hz  respectively,  and  each  transducer  has  a  damping  factor  of 
£= 0.1.  Both  spheres  are  assumed  to  be  neutrally  buoyant  and  the  mass  ratio  mt/(ms+m t)  is  assumed  to  be 
0.01.  The  results  are  plotted  in  Fig.  5  and  normalized  so  that  the  magnitude  of  the  response  is  unity  in  the 
pass-band.  In  this  way,  assumptions  do  not  have  to  be  made  regarding  the  scalar  values  for  0,  C0,  and  Bl. 


FIG.  5.  Response  of  spherical  acceleration  and  velocity  sensors  when  subjected  to  an  acoustic  plane  wave. 

Fig.  5  shows  that  for  each  transducer,  the  magnitude  of  the  response  is  unity  in  the  pass-band  (e.g.,  the 
region  that  is  far  removed  from  the  resonances  in  the  system).  This  figure  also  shows  that  the  phase 
response  is  ‘flat’  within  the  pass-band,  although  a  180°  phase  shift  is  evident  in  the  response  of  the 
accelerometer  above  fs.  Nevertheless,  these  results  provide  the  transducer  designer  with  some  guidelines  on 
how  to  design  the  external  suspension  system  and  select  the  appropriate  transducer  to  meet  the  bandwidth 
requirements  of  an  intensity  probe. 

B.  Experiment 

Positive,  negative,  and  near-neutrally  buoyant  spherical  sensors  were  fabricated  and  tested  so  that  the 
theory  presented  in  the  previous  section  could  be  evaluated.  More  specifically,  commercially  available 
birch  wood,  bronze,  and  urethane  spheres,  having  a  diameter  of  5.08  cm,  were  machined  to  incorporate  an 
Oceana  Sensor  Technologies  Model  AP1BCN  accelerometer.14  These  accelerometers  have  a  nominal 
sensitivity  of  10  mV/g  and  a  mounted  resonance  frequency  of  about  25  kHz.  They  contain  a  shear-type 
sensing  element  and  require  constant  current  DC  power  for  an  internal  FET  based  preamplifier.  They  are 
ideally  suited  for  this  application  because  they  are  designed  to  be  imbedded  into  structures  for  vibration 
monitoring. 

Each  sensor  contains  a  compliant  spring  that  consists  of  a  small  inextensible  tether  that  connects  the 
body  of  the  sphere  to  a  cantilever  beam.  Fig.  6  shows  the  basic  design  of  the  sensor  in  its  physical  and 
lumped  parameter  representations.  Based  on  the  figure,  the  system  can  be  modeled  as  a  mass-loaded 
cantilever  beam  undergoing  rigid  body  translation  in  the  vertical  direction.  Now,  if  the  internal  sensor 
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dynamics  of  the  accelerometer  are  neglected,  then  the  lumped  parameter  resonance  frequency  of  the  system 
is  predicted  by15 

fs 


where  L,  E ,  and  I  are  the  length,  Young’s  modulus,  and  moment  of  inertia  of  the  cantilever  beam.  Note  that 
the  basis  for  neglecting  the  internal  sensor  dynamics  stems  from  the  fact  that  the  resonance  frequency  of  the 
accelerometer  is  well  beyond  the  frequency  range  of  interest. 


3  El 


27C^(ms  +mi)Cs  2;r  y  ( ms  +  mi)L 3 


(5) 


FIG.  6.  Basic  design  of  compliantly  suspended  spherical  acceleration  sensor  in  its  physical  and  lumped  parameter  representations. 

The  logarithmic  decrement  method16  was  used  to  experimentally  determine  the  resonance  frequency 
and  damping  factor  of  each  sphere.  This  was  done  by  submerging  the  sphere  inside  a  vertically  oriented 
column  of  water  and  allowing  it  to  vibrate  freely  once  it  was  displaced  from  its  equilibrium  position.  The 
transient  time-series  data  from  the  accelerometer  embedded  within  the  sphere  was  measured  using  an  HP 
35670A  dynamic  signal  analyzer.  During  the  experiment  the  in-situ  dimensions  of  the  cantilever  beam 
were  measured  with  a  micrometer.  Using  these  parameters  in  conjunction  with  the  physical  properties  of 
the  beam  and  the  sphere  enabled  a  theoretical  prediction  to  be  made  using  Eq.  (5).  The  results  of  the 
experimental  and  theoretical  analyses  are  presented  in  Table  I.  It  can  be  inferred  from  the  table  that  the 
resonance  frequencies  determined  by  experiment  are  in  good  agreement  with  theory  and  that  the  spheres 
are  lightly  damped. 


TABLE  I.  Physical  properties  for  spherical  acceleration  sensors. 


Sphere 

Buoyancy 

ms  (g) 

ps  (g/cm3) 

3p0/(2ps+p0) 

Suspension 

Cs  (mm/N) 

£* 

fs  (Hz)a 

fs  (Hz)b 

Birch  wood 

Positive 

40.0 

0.58 

1.39 

Steel  rod 

4.74 

0.02 

8.8 

8.5 

Bronze 

Negative 

499.0 

7.27 

0.19 

Brass  strip 

1.41 

0.01 

5.3 

5.8 

Urethane 

Near  neutral 

77.5 

1.13 

0.92 

Steel  rod 

10.5 

0.06 

4.0 

4.7 

a  jj - 

Experimental  results,  theoretical  results. 


The  second  part  of  the  experiment  concerns  the  performance  of  the  spheres  in  a  standing  wave  field. 
This  was  done  consistent  with  the  approach  taken  by  Gabrielson2  in  which  the  velocity  of  an  inertial  sensor 
was  compared  to  that  inferred  from  a  reference  hydrophone  positioned  at  the  same  depth  within  a 
waveguide  having  an  acoustic  source  at  one  end  and  a  pressure-release  surface  at  the  other  end.  The 
waveguide  used  for  this  testing  was  made  of  acrylic  and  had  an  inside  diameter  of  10.16  cm,  a  length  of 
100  cm,  and  wall  thickness  of  3.18  mm.  The  walls  of  the  waveguide  are  compliant,  thus  the  sound  speed  is 
substantially  reduced  from  the  bulk  speed  found  in  open  water.  The  sound  speed  was  determined  via 
experiment  to  be  about  381  m/s.17  This  result  compares  favorably  with  theory  for  ‘slow’  waveguides.18"22 

Fig.  7  presents  a  schematic  of  the  experimental  set-up  and  shows  how  the  spheres  and  the  reference 
hydrophone  were  positioned  in  the  waveguide.  A  Reson  TC-4013  hydrophone,  having  a  measured 
sensitivity  of  28.2  ^uV/Pa,  was  used  as  the  reference  standard.  For  each  test  case  the  measurement  depth 
was  7.62  cm  from  the  pressure-release  surface.  The  sound  projector  at  the  bottom  of  the  waveguide  was  a 
Wilcoxon  Research  F4  shaker  that  was  modified  to  incorporate  a  7.62  cm  diameter  aluminum  piston. 
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FIG.  7.  Experimental  set-up  for  tests  performed  in  a  water-filled  waveguide. 

Under  the  assumption  of  a  loss  less  waveguide  having  a  pressure  release  boundary  condition,  it  can  be 
shown  that  the  ratio  of  the  sphere’s  velocity  relative  to  the  acoustic  particle  velocity  is 

Vs  As  pQc  tan  kd 


where  P0  is  the  acoustic  pressure  at  the  measurement  depth,  c  is  the  sound  speed  in  the  waveguide,  k  is  the 
acoustic  wave  number,  and  d  is  the  measurement  depth.  Eq.  (9)  was  implemented  by  measuring  the 
transfer  function  between  the  output  voltages  generated  by  the  accelerometer  in  the  sphere  and  the 
reference  hydrophone.  This  was  done  using  a  HP  3 5 670 A  analyzer  operating  in  frequency  response  mode. 
The  waveguide  was  ensonified  with  band-limited  random  noise  over  the  0-1600  Hz  frequency  range  and 
the  transfer  function  was  measured  using  a  1  Hz  analysis  resolution,  128  rms  averages,  and  a  Hanning 
window.  The  transfer  function  was  subsequently  adjusted  by  the  sensitivities  of  the  transducers. 

The  results  of  the  experiment  are  presented  in  Fig.  8,  which  compares  Eqs.  (2)  and  (6)  over  the  40- 
1000  Hz  frequency  range.  This  figure  shows  that  the  experimental  data  tracks  well  with  theory,  in  that 
varying  the  density  of  the  sphere  causes  a  predictable  change  in  the  corresponding  velocity  amplitude. 
Statistical  analysis  of  the  test  data  revealed  that,  on  average,  the  response  of  the  spheres  deviated  from 
theory  by  up  to  5%. 


FIG.  8.  Result  of  plotting  experimental  data  obtained  on  spheres  in  waveguide  versus  theory  predicted  by  Eq.  2.  The  experimental 
data  is  plotted  with  colored  lines  and  the  theoretical  predictions  are  plotted  as  black  dashed  lines. 

It  is  noted  that  the  analysis  of  the  data  shown  in  Fig.  8  was  limited  to  40-1000  Hz  because  (1)  the  F4 
shaker  could  not  provide  sufficient  levels  of  excitation  below  about  40  Hz  and  (2)  in  Eq.  (6),  the  slope  of 
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tan  kd  changes  rapidly  for  frequencies  greater  than  1000  Hz  and  consequently  puts  excruciating  demands 
on  the  accuracy  of  the  measured  sound  speed,  depth,  and  knowledge  of  any  losses  in  the  waveguide. 

II.  PRESSURE-ACCELERATION  INTENSITY  PROBE 

A.  Practical  intensity  measurements  in  one  dimension 

For  one-dimensional  sound  propagation  having  e*™  dependence,  Fahy6  defines  the  complex  acoustic 
intensity  spectrum  at  a  single  point  in  space  as  C( co)=I( co)+jJ( co) ,  where  I(co)  and  J(cd)  are  the  active  and 
reactive  intensity  spectra: 


/(©)  =  Re{Gp„  (©)}=  2Re{p(ffl)[/(ffl)* }  (7) 

J(co)  =  Imfj  pu  (©)}=  }.  (8) 

In  Eqs.  (7)  and  (8),  Gpu(cri)  is  the  single-sided  time-averaged  rms  cross-spectrum  between  the  acoustic 
pressure  P(cd)  and  the  particle  velocity  U(co).  Intensity  probes  that  contain  discrete  sensors  for  the 
independent  measurement  of  these  parameters  can  be  used  in  conjunction  with  a  standard  two-channel 
dynamic  signal  analyzer  to  compute  the  quantities  I(co)  and  J(co).  For  the  case  of  the  p-u  probe,  Eqs.  (7) 
and  (8)  can  be  implemented  directly,  but  for  the  p-a  probe  the  particle  velocity  is  computed  from  A(co)/jco, 
where  A  (to)  is  the  particle  acceleration. 

B.  Design  and  fabrication 

Two  different  versions  of  the  p-a  probe  were  designed,  fabricated,  and  tested  using  commercial-off- 
the-shelf  parts.  The  designs  conform  to  the  schematic  presented  in  Fig.  9,  which  shows  that  the  sensors  are 
basically  the  same,  but  one  contains  a  slotted  cylindrical  shell  that  surrounds  the  exterior  of  the  probe  and 
each  employ  a  different  type  of  suspension  spring.  The  intent  here  is  to  evaluate  the  performance  of  both 
sensors  so  that  the  effect  of  the  shell  and  the  suspension  springs  can  be  quantified.  An  exploded  view  of 
the  probe  containing  the  slotted  shell  is  shown  in  Fig.  10  and  photographs  of  some  of  the  prototypes  are 
shown  in  Fig.  11. 


Slotted  cylindrical  shell  potted  in  polyurethane  Cylindrical  housing 

made  of  positively 
buoyant  material 

FIG.  9.  Basic  design  concepts  for  the  p-a  probe.  The  probe  on  the  left  is  compliantly  suspended  inside  a  slotted  cylindrical  shell.  The 
probe  on  the  right  is  of  the  same  design  but  contains  a  different  type  of  suspension  spring. 
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The  design  of  the  probe  is  centered  around  embedding  one  bi-directional  accelerometer  and  two  omni¬ 
directional  hydrophones  in  a  cylinder  of  syntactic  foam.  The  accelerometer  is  identical  to  that  discussed 
previously  in  Section  I.  The  hydrophones  are  air-backed  piezoceramic  bender  disks  that  have  a  urethane 
coating  on  the  side  that  is  exposed  to  water.  They  are  electrically  connected  in  parallel  to  effectively 
measure  the  acoustic  pressure  in  the  geometric  center  of  the  probe  where  the  accelerometer  is  located.  The 
electrical  connection  also  enables  the  hydrophones  to  be  relatively  insensitive  to  acceleration.  The 
hydrophones  are  identical  to  those  employed  by  the  AN/SSQ-53D  sonobuoy  and  have  a  nominal  sensitivity 
of  about  -200  dB  re  1  V/^Pa.  They  exhibit  an  in-water  resonance  frequency  of  about  8  kHz.  The  electrical 
signals  from  the  transducers  are  routed  from  the  probe  to  analysis  equipment  via  a  4-conductor  shielded 
cable  having  a  PVC  jacket  (Cooner  Wire:  CW250323). 

The  cylinder  is  machined  from  a  block  of  Syntech  Materials24  AM-24  syntactic  foam,  which  has  a 
nominal  density  of  0.384  g/ cm3.  The  density  is  selected  to  provide  a  buoyancy  force  that  counters  the 
weight  of  the  fully  assembled  sensor.  This  is  done,  to  the  maximum  extent  practical,  so  that  the  probe  is 
neutrally  buoyant  in  water  at  20°C.  As  shown  in  Fig.  10,  the  cylinder  is  comprised  of  three  pieces:  two 
end-caps  and  a  body.  Assembly  steps  are  taken  to  insert  the  hydrophones  into  the  end-caps  and  to  insert 
the  accelerometer  into  the  body.  During  the  final  assembly  step,  the  end-caps  are  inserted  into  the  body. 
Epoxy  resin  is  used  to  maintain  the  integrity  of  the  joint.  The  fully  assembled  probe  has  a  diameter  of  2.54 
cm,  an  aspect  ratio  of  unity,  and  a  mass  of  about  15  g.  These  parameters  translate  into  an  average  density 
of  about  1.17  g/cm3,  which  indicates  that  the  probe  is  slightly  negatively  buoyant. 


Slotted  shell 
(end-cap) 


Suspension  spring 


Slotted  shell 
(body) 


Urethane  potting 
for  hydrophone 


Accelerometer 


Hydrophone 


10-32  heli-coil 
insert 


FIG.  10.  Exploded  view  of  p-a  probe  containing  slotted  cylindrical  shell  and  circular  leaf  spring  suspension. 


The  probe  shown  in  Fig.  10  is  compliantly  suspended  inside  the  slotted  cylindrical  shell  using  two 
circular  leaf  springs  that  rely  on  the  mass-loaded  cantilever  beam  principle.  The  inner  and  outer  bolting 
rings  of  the  springs  are  captured  between  the  mating  surfaces  of  the  syntactic  foam  cylinder  and  the  slotted 
cylindrical  shell,  respectively.  Locater  pins  provide  the  alignment  for  the  inner  bolting  ring  and  a  network 
of  00-90  screws/inserts  facilitate  the  alignment  of  the  outer  bolting  ring.  The  springs  are  designed  to  be 
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compliant  in  the  axial  direction  (e.g.,  the  principle  axis  of  sensitivity)  and  stiff  in  the  radial  direction.  This 
is  done  to  facilitate  an  in- water  probe/suspension  resonance  frequency  of  about  10  Hz.  The  springs  are 
photo-etched  from  0.13  mm  beryllium-copper  shim  stock. 

The  slotted  cylindrical  shell  is  machined  from  a  rod  of  black  nylon-6,  whose  material  properties  exhibit 
a  characteristic  acoustic  impedance  similar  to  that  of  bulk  water.  The  shell  is  used  to  limit  the  axial 
excursions  of  the  probe  and  to  serve  as  a  ‘flow  shield’  when  it  is  encapsulated  in  a  nylon  stocking  or  some 
other  semi-porous  material.  The  object  of  the  flow  shield  is  to  minimize  the  sensor’s  intrinsic  sensitivity  to 
flow  induced  noise  by  displacing  the  turbulent  boundary  layer  (TBL)  away  from  the  body  of  the  probe. 
The  TBL  is  created  when  water  flows  over  the  probe.  A  10-32  heli-coil  insert  is  installed  on  the  flat  face  of 
the  shell  end-cap  to  provide  an  attachment  point  for  an  external  positioning  system.  The  shell  has  a 
diameter  of  36.2  mm,  a  length  of  34.3  mm,  and  a  wall  thickness  of  3.18  mm. 


FIG.  11.  Various  prototype  p-a  probes  and  related  parts. 

As  shown  in  Fig.  11,  a  variant  of  the  design  is  simply  a  ‘bare’  probe  that  employs  a  cantilever  beam  as 
the  external  suspension.  This  embodiment  is  more  or  less  designed  to  be  consistent  with  the  overall 
concepts  presented  in  Section  I.  Moreover,  it  serves  as  an  ‘interim’  design  to  help  quantify  the  effects  of 
slotted  shell,  circular  leaf  springs,  and  acoustic  cavities  that  are  characteristic  of  the  other  probe. 

C.  Experiment 

Several  experiments  were  conducted  to  categorize  the  performance  of  the  p-a  probe.  The  first 
concerns  the  notion  that  a  properly  designed  p-a  probe  should  have  the  ability  to  independently  measure  the 
acoustic  pressure  and  particle  acceleration  without  contamination.  In  other  words,  the  hydrophone  should 
be  able  to  discriminate  against  acceleration  and  the  accelerometer  should  be  able  to  discriminate  against 
pressure.  An  indication  of  whether  this  is  occurring  (or  not)  can  be  deduced  by  scanning  a  standing  wave 
field  for  the  locations  of  the  nodes  and  antinodes.  That  is,  a  perfect  standing  wave  exhibits  pressure  nodes 
and  acceleration  antinodes  at  the  same  location,  and  vice-versa.  Also,  the  acoustic  pressure  and  particle 
acceleration  are  180°  out  of  phase.  Accordingly,  the  p-a  probe  (equipped  with  the  leaf  springs  and  slotted 
shell)  was  positioned  in  the  ‘slow’  waveguide  along  with  a  Reson  TC-4013  hydrophone.  The  waveguide 
was  excited  with  a  250  Hz  tone  and  the  standing  wave  field  was  scanned  over  a  distance  of  500  mm 
starting  from  100  mm  below  the  pressure-release  surface.  Power- spectra  were  acquired  at  5  mm 
increments  using  a  HP  3 5 670 A  spectrum  analyzer.  The  data  were  analyzed  over  the  0-400  Hz  frequency 
range  using  a  1  Hz  analysis  resolution,  16  rms  averages,  and  a  Hanning  window.  Fig.  12  shows  the  result 
of  the  experiment  and  presents  two  separate  plots  which  show:  (1)  the  raw  uncalibrated  voltage  outputs 
from  the  p-a  probe  and  the  reference  hydrophone  as  a  function  of  distance,  and  (2)  the  steady-state  time 
traces  of  the  raw  uncalibrated  voltage  output  from  the  p-a  probe  at  a  distance  76.2  mm  below  the  pressure- 
release  surface.  The  ‘raw  data’  is  presented  so  that  unbiased  conclusions  can  be  made  without  hesitation. 
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It  can  be  inferred  from  Fig.  12  that  the  pressure  measured  by  the  p-a  probe  hydrophone  tracks  very 
well  with  that  of  the  reference  hydrophone.  The  figure  also  shows  that  the  pressure  maxima  and 
acceleration  minima  are  almost  perfectly  aligned.  The  spatial  offset  was  measured  to  be  about  Ax= 20  mm 
which  translates  in  to  a  phase  aberration  of  about  kAx- 4.7°.  This  result  was  verified  by  computing  the 
phase  of  Gpa(co ),  in  that  the  measured  phase  at  250  Hz  was  175.3°,  or  about  4.7°  away  from  perfect  a 
standing  wave.  This  latter  point  can  be  realized  by  viewing  the  time  series  data  presented  in  Fig.  12.  The 
reason  for  the  quasi-perfect  standing  wave  field  is  believed  to  be  associated  with  the  fact  that  the  walls  of 
the  waveguide  are  compliant,  not  rigid,  hence  the  wall  motion  radiates  sound  and  constitutes  a  loss 
mechanism  for  the  system.  Based  on  these  results  it  can  be  concluded  that,  on  a  relative  basis,  the  p-a 
probe  is  working  properly. 


Distance  from  Pressure  Release  Surface  (mm)  Time  (s) 

FIG.  12.  Performance  of  p-a  probe  (equipped  with  leaf  springs  and  slotted  shell)  in  standing  wave  field  at  250  Hz. 

The  second  experiment  concerns  evaluating  the  performance  of  the  accelerometer  in  the  p-a  probe 
using  the  approach  taken  for  the  spherical  sensors  in  Section  I.  The  intent  here  is  to  compare  the  result  of 
Eq.  (2)  versus  Eq.  (6),  however,  this  requires  a  priori  knowledge  of  the  induced  mass  for  the  p-a  probe. 
This  is  an  important  concept  because  the  geometry  of  the  probe  is  not  spherical,  thus  Eq.  (2)  is,  in  general, 
not  applicable.  The  more  general  form  of  Eq.  (1)  is  applicable,  yet  no  closed  form  solution  exists  to 
determine  the  induced  mass  for  the  case  of  a  right  circular  cylinder. 

Accordingly,  the  induced  mass  of  the  p-a  probe  was  measured  using  two  different  techniques.  The 
first  technique  involves  determining  the  resonance  frequency  of  the  compliantly  suspended  probe  using  the 
logarithmic  decrement  method.  This  was  done  consistent  with  the  approach  described  in  Section  I  for  the 
spherical  sensors,  except  the  tests  were  performed  with  the  probe  in  air  and  in  water.  Now,  in  order  to 
simplify  the  analysis,  it  was  assumed  that  the  induced  mass  in  air  was  zero,  thus  Eq.  (5)  can  be  expressed  as 
follows: 


m-  =  m, 


J  ca 

y  few  J 


(9) 


where  mc  is  the  mass  of  the  probe  and  fca  and  fcw  are  the  fundamental  mass-spring  resonances  of  the  probe 
in  air  and  water,  respectively. 

The  second  technique  relies  on  computing  the  in-air  and  in-water  resonance  frequencies  of  the 
hydrophone  employed  by  the  p-a  probe.  This  can  be  realized  when  one  considers  that  the  induced  mass  is 
identical  to  the  radiation  mass  that  the  hydrophone  is  subjected  to  when  it  is  configured  as  a  source  or  a 
receiver.  In  contrast  to  the  first  technique,  this  method  serves  only  as  an  indirect  way  to  determine  the 
induced  mass.  More  specifically,  the  mechanical  mass  of  the  hydrophone  is  not  known,  but  the  ratio  of  the 
resonance  frequencies  can  be  determined.  This  ratio  should  be  the  same  as  that  delineated  in  Eq.  (9),  or, 
fhJfhw=fcJfcw ,  where  fha  and  fhw  are  the  in-air  and  in- water  resonance  frequencies  of  the  hydrophone.  In  this 
way,  the  second  technique  serves  as  a  data  quality  check  for  the  first.  It  is  noted  that  the  resonance 
frequencies  of  the  hydrophone  were  determined  from  an  electrical  impedance  measurement  in  which  the 
transfer  function  between  the  voltage  drop  across  the  transducer  and  the  drive  current  is  measured.25,26 
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Prior  to  performing  the  experiments  on  the  p-a  probes,  the  first  technique  was  validated  using  the 
bronze  sphere  as  the  test  specimen.  Theory  predicts  the  induced  mass  to  be  34.3  g  (e.g.,  mi=2m?> p0/3).  The 
experimental  value  was  determined  to  be  36.9  g,  which  deviates  from  theory  by  7.6%. 

The  results  of  the  induced  mass  experiments  are  presented  in  Table  II.  The  data  obtained  on  the  bare 
probe  indicates  that  the  induced  mass  is  9.1  g  and  that  there  is  very  good  agreement  between  both 
techniques  as  evidenced  by  the  1%  deviation  in  the  term  (fcJfcw)l(fhJfhw )•  The  data  also  indicates  that  Eqs. 
(1)  and  (2)  appear  to  be  fairly  commensurate  despite  the  differences  in  geometry  between  the  cylinder  and 
the  sphere. 

The  results  for  the  case  of  the  probe  with  the  slotted  shell  are  mildly  ambiguous,  in  that  there  was  a 
34%  increase  in  the  induced  mass  relative  to  the  case  of  the  bare  probe,  yet  this  did  not  appear  to  propagate 
into  the  results  predicted  by  Eqs.  (1)  and  (2).  It  is  speculated  that  the  sensitivity  of  these  equations  to 
values  for  the  induced  mass  can  be  rather  slight  provided  the  probe  is  (near)  neutrally  buoyant.  It  is  noted, 
however,  that  the  34%  increase  is  likely  to  be  related  to  the  acoustic  inertance  of  the  fluid  entrained  in  the 
slots. 


Table  H  Results  of  induced  mass  experiments  performed  on  p-a  probes. 


Test  sensor 

mc  (g) 

fcJfcw 

m  (g) 

fJfhw 

( fcJfcW)/(fha/fhW ) 

( m0+mi)/(mc+mi ) 

3pJ(2pc+p0f 

Bare  p-a  probe 

15.0 

1.28 

9.1 

1.29 

0.99 

0.91 

0.90 

p-a  probe  w/slotted  shell 

14.9 

1.39 

13.7 

1.16 

1.20 

0.92 

0.90 

apc  is  the  average  density  of  the  fully  assembled  syntactic  foam  cylinder. 


The  results  of  the  foregoing  analysis  are  now  applied  to  the  experiments  performed  on  the  probes  in 
which  Eqs.  (1)  and  (6)  are  evaluated.  The  bare  probe  was  tested  in  the  ‘slow’  waveguide,  but  the  probe 
with  the  slotted  shell  was  tested  in  a  USRD  G19  calibrator.27  The  G19  is  a  standing  wave  device  made 
from  an  10.16  cm  diameter  aluminum  tube.  The  sound  speed  in  the  G19  is  estimated  to  be  about  1219 
m/s.28  The  G19  was  used  because  it  was  thought  that  the  acoustic  cavities  in  the  slotted  shell  could 
resonate  and  consequently  contaminate  the  measurements.  Increasing  the  sound  speed  of  the  fluid  helps  to 
push  these  resonances  out  of  the  frequency  range  of  interest. 

The  results  of  this  testing  are  presented  in  Fig.  13.  This  figure  shows  that  the  data  obtained  with  the 
bare  probe  is  commensurate  with  Eq.  (1)  except  for  regions  of  the  spectrum  where  spurious  resonances  are 
evident.  It  is  speculated  that  the  source  of  these  resonances  is  associated  with  torsional  modes  and  higher 
order  transverse  bending  modes  of  the  suspension  spring.29  Torsional  modes  are  suspected  because  the 
attachment  points  of  the  cantilever  beam  and  the  signal  cable  are  perpendicular  (see  Fig.  11).  Thus  the 
eccentric  weight  of  the  cable  creates  a  torque  on  the  probe  which  is  ultimately  transferred  to  the  spring. 
Transverse  bending  modes  are  suspected  because  the  system  can  be  only  modeled  with  lumped  parameters 
for  low  frequencies.  In  other  words,  if  the  fundamental  bending  mode  of  a  mass-loaded  cantilever  beam  is 
10  Hz  (as  is  the  case  here),  then  lumped  parameters  should  be  sufficient  to  describe  the  dynamics  of  the 
system  up  to  about  a  decade  above  resonance.  Beyond  this  point  it  is  reasonable  to  assume  that  higher 
order  bending  modes  may  become  evident  in  the  response.  However,  suspension  springs  that  include  an 
appropriate  level  of  damping  are  anticipated  to  resolve  this  issue. 


FIG.  13.  Results  of  testing  p-a  probes  to  assess  whether  their  velocity  amplitude  is  commensurate  with  theory  for  a  standing  wave. 


NUWC-Newport  Workshop  on  Directional  Acoustic  Sensors 
April  17-18,  2001 


J.  McConnell:  p-a  intensity  probe 
Page  12  of  13 


Fig.  13  also  shows  that  the  data  acquired  on  the  probe  with  the  slotted  shell  does  not  agree  with  theory 
and  appears  to  be  controlled  by  a  resonance  response  occurring  at  about  750  Hz.  This  result  was  repeatable 
over  time  and  was  also  observed  when  another  probe  of  the  same  design  was  tested  in  the  G19.  At  the 
present  time,  the  source  of  the  resonance  is  unclear,  although,  wave  effects  in  the  suspension  springs, 
structural  modes  of  the  slotted  shell,  and  acoustic  cavity  modes  of  the  slotted  shell  are  all  considered  to  be 
possible  sources.  Work  is  on  going  to  resolve  these  issues. 

III.  SUMMARY  AND  CONCLUSIONS 

The  dynamics  of  a  compliantly  suspended  sphere  containing  an  accelerometer  were  evaluated  using 
theory  and  experiment.  Results  show  that  the  performance  of  the  sensor  can  be  predicted  using  lumped 
parameter  theory,  in  that  varying  the  density  of  the  sphere  causes  a  predictable  change  in  the  corresponding 
velocity  amplitude.  This  result  is  specific  to  the  pass-band,  which  covers  the  frequencies  between  the 
resonances  in  the  system. 

The  results  obtained  for  the  spherical  sensors  were  extended  to  the  case  of  developing  a  near-neutrally 
buoyant  p-a  intensity  probe.  Two  p-a  probes  were  evaluated,  one  that  was  compliantly  suspended  in  a  free- 
flooding  slotted  cylindrical  shell  and  another  that  employed  a  simple  cantilever  beam  suspension.  In 
essence,  the  only  difference  between  the  two  sensors  was  how  they  were  packaged.  The  former  is  far  more 
complex  than  the  latter. 

Various  tests  were  performed  on  the  p-a  probes  to  assess  their  performance  in  a  standing  wave  field. 
The  basic  design  appears  to  be  viable,  in  that  (1)  the  hydrophone  is  responsive  to  the  acoustic  pressure  and 
discriminates  against  acceleration  and  (2)  the  accelerometer  is  responsive  to  the  particle  acceleration  and 
discriminates  against  pressure. 

The  probe  containing  the  simple  cantilever  beam  suspension  (e.g.,  the  ‘bare’  probe)  exhibited  the  best 
performance  of  the  two  probe  designs  that  were  evaluated.  This  occurred  because  an  in-band  resonance 
corrupted  the  data  obtained  on  the  probe  that  had  the  slotted  cylindrical  shell. 

It  is  clear  from  this  study  that,  novel  and  ‘acoustically  friendly’  packaging  concepts  are  critical  to  the 
successful  development  of  p-a  probes.  This  is  particularly  important  when  the  probes  are  required  to  have 
high  fidelity  over  a  large  bandwidth. 
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Abstract.  A  computational  study  is  described  where  a  2-D  elliptic  cylinder  is  insonified  by  a  plane, 
monochromatic  acoustic  wave.  The  elliptic  cross  section  of  the  cylinder  has  a  fineness  ratio  of 
5:1,  the  incidence  angle  of  the  plane  wave  is  60°  relative  to  the  major  axis  of  the  ellipse,  and  ka  = 
20,  where  a  is  the  major  axis  of  the  elliptic  cross  section  and  k  is  the  acoustic  wavenumber.  The 
calculations  are  performed  using  the  finite  element  method  of  solution  for  partial  differential 
equations.  The  MATLAB®  Partial  Differential  Equations  Toolbox  was  used  to  formulate  and 
solve  the  Helmholtz  equation  with  reflection-free  conditions  imposed  on  the  computational  outer 
boundary,  and  rigid  conditions  imposed  on  the  surface  of  the  scatterer.  Of  particular  practical 
interest  in  this  study  is  the  spatial  distribution  of  the  total  active  acoustic  intensity,  i.e.,  the  sum 
of  the  incident  and  scattered  intensity  components.  Active  intensity  amplitude,  and  the  phase 
between  pressure  and  particle  velocity,  are  computed  and  compared  to  pressure  amplitude  only. 

The  results  show  that  there  is  significant  phase  distortion  in  the  forward  scattered  direction  that 
could  be  useful  in  localizing  objects  in  active  bistatic  operations  if  p-u  type  acoustic  intensity 
probes  were  employed. 

PACS  numbers:  43.20.Fn,  43.30.Yj 

I.  INTRODUCTION 

Acoustic  scattering  from  objects  is  an  important  area  of  underwater  acoustics.  In 
practical  applications,  sound  is  created  by  a  source  operating  under  either  steady-state  or 
transient  conditions;  it  can  be  radiating  discrete-frequency  or  broadband  sound.  The  sound 
scattered  by  the  object  is  received  by  sensors  located  at  various  positions  around  the  object. 

When  the  receiver  is  co-located  with  the  source,  the  scattering  is  known  as  monostatic.  When 
the  sound  is  received  elsewhere,  the  situation  is  called  bistatic.  Pressure  hydrophones  are  the 
standard  sensor  used  in  these  applications.  Acoustic  intensity,  on  the  other  hand  is  the  product  of 
acoustic  pressure  and  acoustic  particle  velocity.  It  is  a  vector  quantity  that  describes  the  acoustic 
energy  per  unit  area  per  unit  time  that  is  radiated  from  an  acoustic  source,  or  is  scattered  by  an 
object  insonified  by  that  source.  It  is  not  a  traditional  measure  in  the  underwater  acoustics  arena, 
primarily  because  there  were  no  reliable  means,  until  now,  to  measure  it.  Because  state-of-the-art 
underwater  acoustic  intensity  vector  sensors  can  measure  the  real-time  phase  shifts  between 
pressure  and  particle  velocity,  a  whole  new  area  of  acoustic  data  interpretation  is  now  available 
through  use  of  these  novel  sensors. 

The  use  of  intensity  measurement  methodology  in  acoustic  scattering  systems  is  not  well 
understood  at  this  time.  Because  it  is  very  difficult  to  distinguish  between  the  direct  field  and  the 
scattered  field  in  bistatic  operations  where  the  target  is  in  the  direct  line  between  receiver  and 
source,  there  is  a  need  for  research  to  determine  whether  intensity  sensors  offer  any  solution  to 


2 


this  problem.  The  primary  emphasis  of  the  research  described  in  this  paper  is  to  determine, 
theoretically,  the  advantages,  or  disadvantages,  of  using  underwater  acoustic  intensity  sensors  in 
underwater  scattering  scenarios.  The  issue  is  whether  the  acoustic  intensity,  determined  by  a 
direct  measurement  of  the  particle  velocity,  the  pressure,  and  the  phase  between  them,  provides 
additional  new  information  on  the  scattering  characteristics  of  an  object  over  that  obtained  from 
conventional  pressure  or  acoustic  velocity  hydrophones  used  alone. 

II.  APPROACH 

As  sketched  in  Fig.  1,  a  rigid,  2-D  elliptic  cylinder  is  the  scattering  object.  The  scattered 
pressure,  particle  velocity,  and  acoustic  intensity  are  determined  from  numerical  solutions  of  the 
Helmholtz  equation.  The  MATLAB®  Partial  Differential  Equations  Toolbox,  that  utilizes  the 
finite  element  method  (FEM)  of  computation,  is  used  for  this  exercise.  We  concentrate  on  the 
total  acoustic  intensity  distributed  throughout  the  field.  This  is  the  complex  sum  of  the  incident 
intensity  from  an  arbitrarily  located  distant  source  and  the  intensity  scattered  from  the  surface.  It 
would  be  representative  of  the  intensity  signal  measured  when  the  sound  source  is  broadcasting  a 
single-frequency  wave  under  continuous,  steady-state  conditions. 


FIG.  1.  Scattering  of  a  plane  wave  by  an  elliptic  cylinder. 

We  assume  that  a  pressure-acoustic  velocity  (p-u  type)  intensity  sensor1  measures  the 
field.  These  probes  measure  intensity  directly  and  make  no  estimate  of  pressure  or  particle 
velocity  using  finite-difference  solutions  of  either  the  continuity  equation  or  the  linearized 
Euler’s  equation,  as  is  necessary  in  velocity-velocity2,3  ( u-u )  and  pressure-pressure4  (p-p)  probes, 
respectively.  The  phase  between  p  (acoustic  pressure)  and  u  (acoustic  particle  velocity),  or  p  and 
a  (acoustic  particle  acceleration)  is  always  available.  These  two  quantities  are  measured  at  the 
same  point  in  space;  thus,  we  can  determine  the  phase  field  at  all  locations  around  the  scatterer. 

The  MATLAB®  triangular  mesh  grid  generator  was  used  to  create  the  scatterer  and 
thousands  of  field  points.  We  require  at  least  10  mesh  points  per  acoustic  wavelength.  A  5:1 
fineness  ratio  (a/b  in  Fig.  1)  was  selected  for  the  elliptic  cylinder.  For  acoustic  radiation  and 
scattering  problems,  the  wave  equation,  which  is  hyperbolic,  is  solved.  The  incident  wave  is  a 
plane  wave  propagating  at  angle  6  to  the  major  axis  of  the  ellipse.  The  frequency  parameter  is 
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ka,  where  k  is  the  acoustic  wavenumber.  The  particle  velocity  is  calculated  using  Euler’s 
equation,  which  is  then  conjugated  and  multiplied  by  one-half  the  pressure  to  get  intensity: 

I  =  i  pu  (1) 

The  real  part  of  this  expression  is  the  active  intensity,  and  the  imaginary  part  is  the  reactive 
intensity.  In  the  numerical  calculations,  we  emphasize  the  magnitude  of  the  active  intensity  and 
the  phase  between  pressure  and  particle  velocity. 

The  FEM  modeling  and  coding  were  verified  by  comparing  the  numerical  calculations 
with  the  theoretically  exact  calculations5  of  an  infinitely-long,  rigid  cylinder  insonified  by  a  plane 
wave.  The  agreement  between  the  two  approaches  was  nearly  perfect  for  ka  =  3  and  using  15 
terms  in  the  exact  modal  series  solution. 

III.  RESULTS 

Calculations  of  the  active  intensity  and  pressure  due  to  the  scattering  of  a  plane  wave  by  a 
rigid,  2-D  elliptic  cylinder  are  considered  for  several  plane  wave  incidence  angles  and 
frequencies.  Presented  here  are  the  results  for  ka  =  20  and  9  =  60°.  Figure  2  shows  the 
computed  total  pressure  field.  Periodic  color  changes  represent  pressure  magnitude  changes  as 
the  compression  and  rarefaction  parts  of  the  acoustic  wave  pass  various  points  in  space.  The 
computational  boundaries  are  large  enough  so  that  the  patterns  shown  near  them  represent  far- 
field  patterns.  The  incident  pressure  field  is  described  by  the  pattern  in  the  vicinity  of  the  arrow, 
which  indicates  the  direction  from  the  source.  It  is  noted  that  the  pattern  in  the  forward-scattered 
farfield  is  very  similar  to  that  of  the  incident  field.  That  similarity  makes  the  detection  of  the 
elliptic  cylinder  in  the  forward- scattered  direction  difficult,  if  a  pressure  hydrophone  is  used  for 
the  sensor  and  if  it  is  located  in  the  farfield  of  the  scatterer. 

Turning  our  attention  now  to  the  computed  intensity,  Fig.  3  shows  the  magnitude  of  the 
active  acoustic  intensity  field.  The  vectors  show  precisely  how  the  acoustic  energy  interacts  with 
the  surface  and  diffracts  around  it  to  form  a  distorted  intensity  field  nearly  everywhere  around  the 
object.  Of  particular  interest  is  how  different  the  forward-scattered  intensity  field  is  from  the 
incident  field.  This  means  that  an  intensity  measurement  performed  in  the  forward-scattered 
direction  may  reveal  features  not  identified  by  acoustic  pressure  measurements  alone.  Although 
forward-scattered  highlights  appear  in  the  intensity  field  plots  of  Fig.  3,  they  are  actually 
“lowlights”  because  the  magnitudes  are  smaller  than  that  of  the  incident  intensity;  i.e.,  the 
forward-scattered  intensity  reveals  a  strong  shadow  region. 

The  beam  pattern  expresses  farfield  magnitude  in  a  polar  format.  The  active  intensity 
beam  pattern  shows  little  more  information  than  the  beam  pattern  for  the  total  acoustic  pressure, 
e.g.,  Fig.  4.  Thus,  simply  measuring  the  magnitude  of  the  active  acoustic  intensity  in  the  farfield 
may  not  provide  any  advantage  over  measuring  pressure  (or  particle  velocity)  alone.  The 
important  information  is  in  the  phase  between  the  pressure  and  particle  velocity. 
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FIG.  2.  Total  acoustic  pressure  scattered  by  an  elliptic  cylinder. 
The  source  is  situated  in  the  lower  left  as  indicated  by  the  arrow 
at  0  =  60°. 


Arrow:  Direction  of  Total  Active  Intensity  (  ka=20,  a=1,  b=0.2,  angle=60  ) 
Background  Color:  Magnitude  in  dB//lncomming  Intensity 


FIG.  3.  Total  active  acoustic  intensity  scattered  by  an  elliptic 
cylinder.  The  source  is  situated  in  the  lower  left  as  indicated  by  the 
origin  of  the  vectors;  6  =  60°. 
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Beam  Pattern  of  Total  Pressure  and  Radial  Active  Intensity  in  dB 


FIG.  4.  Beam  patterns  for  the  magnitude  of  the  active  intensity  in  the  forward 

scattered  directions  compared  to  that  of  the  total  pressure;  9  =  60°. 

Figure  5  shows  the  phase  (radians)  between  the  x-component  of  the  acoustic  particle  velocity  and 
the  pressure  at  all  field  points  around  the  scatterer.  The  phase  of  the  incident  wave  is  simply  kr 
which  grows  linearly  with  distance,  r.  The  plot  shows  this  linear  growth  as  a  sawtooth  wave 
because  the  range  of  the  phase  is  bounded  between  (-k,k).  A  remarkable  deviation  from  the 
sawtooth  pattern  is  evident  in  the  forward- scattered  direction.  This  phase  distortion  should  be 
easy  to  measure  with  a  p-u  or  a  p-a  type  underwater  acoustic  intensity  probe.  Calculations  of  the 
phase  between  the  y-component  of  particle  velocity  and  pressure  show  features  and  trends 
analogous  to  those  of  Fig.  5. 

IV.  CONCLUSIONS 

Presented  here  are  numerical  calculations  of  the  acoustic  scattering  of  a  plane  wave  by  a 
rigid,  2-D  elliptic  cylinder.  The  particle  velocity  and  pressure  are  computed  and  combined  to 
form  the  active  intensity.  The  interpretation  of  these  results  in  terms  of  using  intensity  sensors  in 
scattering  systems  is  as  follows.  Under  steady-state,  pure  tone  source  excitation  in  a 
reverberation-free  environment,  the  total  (incident  +  scattered)  active  intensity  field  contains 
valuable  information  in  the  phase  between  particle  velocity  and  acoustic  pressure.  In  the 
forward-scattered  direction,  this  phase  fluctuates  rapidly  with  location  suggesting  that  an  object 
could  be  easily  identified  in  this  direction  using  a  p-u  type  intensity  probe. 
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Phase  between  Pressure  and  Velocity  in  x-direction  ( ka=20,a=1  ,b=0.2,angle=60) 


FIG.  5.  Phase  between  the  x-component  of  particle  velocity  and  the  acoustic 
pressure  at  collocated  points  in  space  during  the  scattering  of  a  plane  wave  at 
60°  to  an  elliptic  cylinder. 
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Abstract.  This  brief,  provided  at  the  Directional  Acoustic  Sensors 
Workshop  on  17,18  April  2001,  is  intended  to  give  the  attendees  a 
description  and  an  historical  review  of  the  Directional  Frequency  Analysis 
and  Recording  (DIFAR)  Acoustic  Subsystem,  which  is  either  currently  or 
previously  used  by  several  Navy  air  antisubmarine  warfare  (ASW) 
platforms. 


I.  INTRODUCTION 

The  focus  of  the  brief  is  on  the  directional  acoustic  sensors  employed  and  how  their 
data  are  handled:  from  its  hydrophone’s  inputs,  through  the  AN/SSQ-53  DIFAR 
sonobuoy,  the  radio  frequency  (RF)  data  transmission  link  to  the  monitoring  aircraft,  and 
the  aircraft’s  DIFAR  Acoustic  Signal  Processing  and  Recording  Subsystems.  Due  to  the 
classified  nature  regarding  many  of  the  details  related  to  the  DIFAR  Acoustic  Signal 
Processing  Subsystems,  and  the  requirement  that  the  material  presented  at  this  workshop 
be  unclassified,  the  description  of  the  DIFAR  subsystem  provided  in  this  brief  given  at 
the  workshop,  and  herein,  is  both  cursory  and  nonspecific  in  nature. 

II.  SUMMARY 

DIFAR  acoustic  sensors  have  been  in  operational  use  by  the  Navy  in  its  airborne 
antisubmarine  warfare  (ASW)  platforms  since  1969.  They  were  first  introduced  to,  and 
used  by,  the  Navy’s  P3C  Orion  Aircraft.  DIFAR  capability  has  also  been  integrated  into 
other  air  ASW  platforms,  such  as  the  S2G,  S3  A,  and  S3 13  fixed-wing  aircraft,  and  by  the 
SH-60B  (LAMPS),  SH-60F  (CV  Helo),  and  the  SH-60R  (ALFS)  rotary  wing  aircraft.  The 
Navy’s  DIFAR-capable  air  ASW  platforms,  which  will  be  in  operational  use  for  the  near 
future,  are  the  P3C  and  the  SH-60R.  The  SH-60R  helicopter  is  the  replacement  for  both 
the  SH-60B  and  SH-60F  platforms.  The  SH-60R  incorporates  the  capabilities  of  both  the 
SH-60B  and  the  SH-60F  platforms,  and  also  incorporates  upgrades  to  the  acoustic 
subsystems  of  both  the  SH-60B  and  the  SH-60F. 

Other  DIFAR-like  acoustic  sensors  currently  in  use  by  the  Navy’s  air  ASW  platforms 
include  the  passive  Vertical  Line  Array  DIFAR  (VLAD)  (AN/SSQ-77),  and  the  active/ 
passive  Directional  Command  Activated  Sonobuoy  System  (DICASS)  (AN/SSQ-62) 
sonobuoy.  The  VLAD  sensor  employs  a  vertical  line  array  of  omnidirectional  hydrophones 
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to  form  a  steered  beam  in  the  vertical  direction  for  use  as  the  omni  channel  component  of  a 
composite  DIFAR  acoustic  signal.  The  DICASS  sensor  is  an  active/passive  DIFAR-like 
sonobuoy,  which  is  used  by  air  ASW  platforms  for  the  precise  target  localization  and  attack 
phases  of  the  ASW  problem. 

The  types  of  directional  (i.e.,  pressure  gradient  (PG))  hydrophones  used  in  DIFAR/ 
DICASS  sonobuoys  are  listed  and  discussed.  They  include 

Type  1:  Ceramic  vane. 

Type  2:  Moving  coil. 

Type  3:  Closely  spaced  pairs  of  omnidirectional  hydrophones  (OHs)  (which 
measure  the  integrated  pressure  over  their  surface), 

Type  4:  Segmented  ceramic  cylinder,  and 
Type  5:  Ceramic  cylinder  (multimode). 

Types  1  and  2  have  been  predominately  used  for  passive  sensor  (i.e.,  DIFAR/VLAD) 
application;  use  of  types  4  and  5  has  predominated  for  the  DICASS  application.  The 
rationale  for  PGH-type  selection  vs.  sensor  application  and  the  effects  of  flow  noise  upon 
sensor  performance  are  also  discussed  in  the  presentation. 

The  presence  of  fluid  flow  in  the  immediate  vicinity  of  a  PGH  sensor  causes  local 
particle  motion,  thus  acting  as  a  nearfield  source  of  additional  acoustic  noise.  This  effect 
is  especially  true  at  low  frequencies,  thus  reducing  the  signal-to-noise  ratio  of  acoustic 
signals  arriving  from  distant  (i.e.,  farfield)  sources  at  these  frequencies.  When  this  effect 
is  observed,  the  acoustic  sensor  is  said  to  be  “flow-noise  limited”  at  these  frequencies. 

As  the  local  relative  fluid  flow  velocity  increases,  the  frequency  below  which  the  acoustic 
sensor  becomes  “flow-noise  limited”  also  increases.  The  effect  of  this  flow-induced  noise 
also  reduces  the  effective  dynamic  range  of  the  sensor  with  respect  to  acoustic  signals 
arriving  from  farfield  sources  at  all  frequencies.  As  the  flow  velocity  is  increased  still 
further,  the  sensor’s  effective  dynamic  range  is  also  further  reduced,  until,  finally,  the 
sensor’s  output  becomes  saturated/nonlinear,  and,  in  effect,  unusable.  For  these  reasons, 
flow-noise  reduction  components,  such  as  compliant  hydrophone  suspension  systems  to 
reduce  vertical  acoustic  sensor  motion,  plus  “flow  socks/screens”  to  reduce  the  effects  of 
local  turbulence,  are  incorporated  into  DIFAR  sonobuoy  designs. 

The  methods  employed  in  the  DIFAR  composite  signal  multiplexing  process 
performed  within  the  sonobuoy  and  the  DIFAR  composite  signal  demultiplexing  process 
performed  within  the  monitoring  aircraft’s  acoustic  signal  processor  are  outlined  and 
described  in  moderate  detail  within  the  briefing  provided  at  the  Workshop.  A  brief 
discussion  of  one  of  the  methods  used  to  compute  the  magnetic  bearing  associated  with 
detected  signals  is  also  provided. 
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DIFAR  System  Overview 
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FRONT  VIEW 


Ceramic  Vane  Directional  Hydrophone  Layout 

Note  :  Addition  of  an  Omni  Hydrophone  (e.g.  Ceramic  Disk)  to 

the  Central  Mass  unit  is  an  example  of  a  DIFAR  Hydrophone 
array  used  in  some  versions  of  the  AN/SSQ-53  Sonobuoy 
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Abstract.  Analytical  models  of  a  polyurethane-encapsulated  accelerometer- 
based  velocity  sensor,  along  with  its  mounting  bracket,  are  developed  to 
determine  the  interaction  between  the  various  components  of  the  sensor 
assembly.  These  models  are  used  to  understand  the  function  of  the  polyurethane 
encapsulant  and  the  relative  effects  of  accelerometer  buoyancy,  and  to  optimize 
the  design  of  the  sensor-mounting  bracket.  The  two  primary  functions  of  the 
brackets  are  to  isolate  the  sensor  from  ship  structural  vibration  and  to  provide 
sufficient  compliance  for  the  sensor  to  respond  to  an  acoustic  wave.  Based  on 
these  criteria,  the  designs  of  the  sensor  elements  and  the  mounting  brackets  are 
optimized.  Attempts  are  made  to  correlate  the  model  predictions  with  the  results 
of  recent  sensor  tests. 

I.  INTRODUCTION 

The  work  reported  in  this  paper  deals  with  two  topics  within  the  conformal  hull 
array  backfit  program:  sensor  design  and  operation,  and  a  re-design  of  the  sensor 
mounting  brackets.  It  was  identified  early  in  the  program  that  the  anticipated  use  of 
velocity  sensors  in  place  of  hydrophones  would  impose  new  constraints  on  the  sensor 
mounting  brackets.  Support  bracket  vibration  was  not  a  significant  noise  source  when  the 
sensor  was  a  pressure  hydrophone,  owing  to  the  typically  low  acceleration  sensitivities  of 
hydrophones.  The  situation  changes  dramatically  in  the  case  of  a  velocity  sensor,  and  the 
mount  must  provide  significant  isolation  from  ship  structural  vibration.  The 
determination  of  the  required  amount  of  isolation  is  discussed  below. 

The  sensor  considered  in  this  study  is  being  designed  and  fabricated  by  EDO 
Corporation.  Although  several  implementations  of  the  sensor  are  being  considered,  most 
consist  of  a  CAVES-type  accelerometer  encapsulated  in  syntactic  foam  to  achieve  neutral 
buoyancy.  This  neutrally  buoyant  sensor  is  a  single-axis  sensor,  and  three  must  be 
integrated  into  a  single  package  in  order  to  achieve  a  triaxial  sensor.  This  has  been 
accomplished  by  casting  the  sensors  together  using  polyurethane.  Part  of  the  work 
reported  here  addressed  the  importance  and  impact  of  several  design  features  of  the 
triaxial  sensor  package. 

II.  SENSOR  MODEL 

In  order  to  examine  the  fundamental  interactions  between  the  sensor  components,  a 
simple  spherically  symmetric  sensor  model  was  developed.  A  plane  wave  in  the 
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surrounding  water  is  incident  upon  the  sensor  package.  The  details  of  the  model 
development  are  contained  in  a  paper  that  is  currently  under  review  [1],  This  model  is  an 
extension  to  the  earlier  work  done  by  Leslie  et  al.,  who  considered  a  rigid,  yet  neutrally 
buoyant,  package  [2]. 

The  results  of  this  analysis  showed  that  at  low  frequencies,  the  entire  polyurethane 
assembly  can  move  essentially  as  a  rigid  body  with  the  particle  velocity  of  the  incident 
wave.  In  this  case,  the  sensor  obviously  moves  with  the  particle  velocity  as  well.  As  the 
frequency  increases,  internal  dynamics  involving  deformation  of  the  polyurethane  are 
observed.  For  a  freely  suspended  sensor,  the  first  internal  resonance  corresponds  to  a 
nodal  point  at  the  center  of  the  sphere  -  this  implies  a  zero  velocity  response  for  the 
accelerometer.  Alternatively,  at  a  fixed  frequency,  there  is  a  critical  polyurethane  coating 
thickness  above  which  the  internal  sensor  dynamics  appear.  This  fact  may  be  seen  from 
Fig.  1,  which  shows  the  magnitude  and  phase  of  the  ratio  of  sensor  velocity  to  particle 
velocity  as  a  function  of  frequency  for  the  case  of  a  relatively  thick  coating.  Based  on 
these  results,  it  was  deduced  that  the  excitation  of  internal  dynamics  should  be  avoided  in 
the  sensor  design.  This  conclusion  indicates  the  need  for  relatively  thin  polyurethane 
coatings. 


FIG.  1.  Velocity  ratio  magnitude  and  phase  vs.  frequency  for 
1.0-inch  diameter  sensor  and  1.25-inch  coating. 


Although  in  the  ideal  case  the  foam-encapsulated  accelerometer  is  neutrally 
buoyant,  packaging  constraints  sometimes  make  this  impossible.  The  model  was  used  to 
examine  the  effects  of  negative  buoyancy.  It  was  found  that  doubling  of  the  sensor 
weight  results  in  an  approximate  2-dB  reduction  in  the  velocity  response.  The  variations 
in  the  phase  response  were  generally  quite  small,  within  a  few  degrees. 

If  the  internal  dynamics  of  the  assembly  are  negligible,  and  the  assembly  moves 
essentially  as  a  rigid  body,  then  the  effects  of  a  sensor  mount  on  the  sensor  performance 
may  be  easily  examined.  In  fact,  this  is  a  well-known  problem,  treated  by  Leslie  and 
others.  If  the  suspension  is  imagined  to  be  spring-damper  combination,  then  the  sensor 
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velocity  response  may  be  calculated  as  a  function  of  the  ratio  between  the  signal 
frequency/ and  the  in-vacuo  natural  frequency  of  the  sensor-spring  system,  f_s.  These 
results  are  shown  in  Fig.  2  for  three  different  values  of  the  damping  ratio.  It  is  clear  that 
accurate  measurement  of  the  particle  velocity  requires  a  compliant  sensor  mounting 
system.  The  lowest  signal  frequency  of  interest  should  be  at  least  three  times  greater  than 
the  in-vacuo  mounted  sensor  natural  frequency.  This  fact  imposes  a  constraint  on  the 
sensor  mounting  brackets,  since  the  brackets  must  be  sufficiently  compliant  to  allow  the 
sensor  package  to  move  with  the  particle  velocity.  As  shown  in  the  following  section, 
there  are  other  constraints  driving  the  mount  design  as  well. 


Frequency  Ratio  (f/f_s) 

FIG.  2.  Velocity  ratio  magnitude  vs.  frequency  ratio  for 
three  values  of  suspension  damping. 

III.  SUPPORT  BRACKET  DESIGN 

As  mentioned  earlier,  the  sensor  support  bracket  must  also  provide  isolation  for  the 
sensor  from  ship-borne  vibrations.  In  order  to  address  this  issue,  it  was  first  necessary  to 
determine  how  the  necessary  amount  of  isolation  should  be  defined.  It  was  decided  that 
the  sensor  output  due  to  ambient  background  noise  in  the  water  should  be  no  smaller  than 
the  sensor  output  due  to  transmitted  ship  structural  vibrations.  This  defines  a  situation 
that  may  be  described  as  ‘ambient  noise  limited’,  although  this  is  a  bit  of  a  misnomer 
since  there  are  other  noise  sources  that  are  omitted  from  the  analysis.  In  order  to 
implement  this  analysis,  it  is  necessary  to  define  both  the  background  noise  levels  and  the 
levels  of  ship  structural  vibration. 

Ambient  ocean  noise  levels  may  be  obtained  from  standard  sources  based  on  sea 
state  and  shipping  density.  For  purposes  of  demonstrating  the  analysis  technique, 
conditions  corresponding  to  sea  state  2  and  shipping  density  4  are  assumed.  Quantifying 
ship  vibration  levels  at  the  sensor  mount  locations  are  more  difficult.  The  existing 
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hydrophones  are  attached  to  a  truss  system  that  wraps  around  the  bow  of  the  submarine. 
This  truss  is  of  variable  geometry  depending  on  location,  and  is  structurally  tied  to 
supports  such  as  bulkheads  and  frames.  There  is  limited  experimental  data  available  that 
measures  the  vibration  levels  at  these  hydrophone-mounting  locations.  Recent  tests  were 
performed  to  make  some  representative  measurements  at  several  ship  speeds.  These  can 
be  used  to  implement  the  analysis  described. 

The  analysis  proceeds  as  follows.  The  ambient  noise  in  the  ocean  induces  a  certain 
output  level  from  the  sensor.  This  assumes  that  the  sensor  suspension  provides  sufficient 
compliance,  as  discussed  above.  Then,  for  a  given  level  of  ship  structural  vibration,  one 
can  calculate  the  required  amount  of  isolation  the  sensor  mount  must  provide. 
Alternatively,  one  may  calculate  the  sensor  mount  stiffness  or  compliance  necessary  to 
achieve  this  ambient-noise  limited  condition.  This  calculation  is  then  performed  over  the 
operating  band  of  the  sensor,  nominally  100  -  2000  Hz.  Thus,  for  a  given  sensor  mass, 
one  may  calculate  the  largest  mounting  stiffness  that  can  be  tolerated. 

This  value  of  stiffness  may  be  compared  with  the  value  necessary  to  satisfy  the 
sensor  operation  constraints  described  in  the  preceding  section.  Finally,  these  values  of 
stiffness  may  be  compared  against  each  other  to  develop  an  engineering  design  tool. 
Typical  results  from  this  analysis  are  shown  in  Fig.  3  for  the  case  of  low  ship  speed,  and 
in  Fig.  4  for  high  speed.  In  these  figures,  the  solid  dark  line  indicates  the  boundary 
between  self-noise  limited  conditions  (mounting  stiffness  values  above  the  solid  line)  and 
ambient-noise  limited  conditions  (mounting  stiffness  values  below  the  dark  line). 
Similarly,  the  dashed  line  shows  the  maximum  allowable  mounting  stiffness  that  will 
allow  for  accurate  velocity  sensor  performance  based  on  the  results  shown  in  Fig.  2. 
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FIG.  3.  Sensor  mount  stiffness  requirement  for  low  ship  speed. 
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High  Speed  Case 
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FIG.  4.  Sensor  mount  stiffness  requirement  for  high  ship  speed. 


As  a  point  of  reference,  the  existing  mounting  brackets  that  hold  the  DT-276 
pressure  hydrophones  in  the  conformal  array  have  a  mounting  stiffness  of  approximately 
800  lb/in.  It  is  clear  from  these  results  that,  for  these  two  cases,  the  mounting  stiffness  is 
constrained  by  the  requirement  for  proper  sensor  performance,  i.e.,  a  low  mounted  natural 
frequency.  If  this  constraint  is  met,  then  the  conditions  for  ambient-noise  limited 
operation  will  also  be  met. 

Using  the  low-speed  results  shown  in  Fig.  3,  and  assuming  a  nominal  sensor 
assembly  weight  of  2.5  lb,  a  recommended  mounting  stiffness  of  100  lb/in.  may  be 
deduced.  To  demonstrate  that  this  stiffness  may  be  readily  achieved,  a  parametric 
redesign  of  the  current  hydrophone-mounting  bracket  was  undertaken.  Finite  element 
models  of  the  existing  bracket  were  generated  and  validated  to  assure  their  accuracy.  The 
critical  dimensions  were  then  systematically  varied  to  produce  a  design  that  provided  the 
necessary  compliance.  Pictures  of  the  original  mount  (on  the  left)  and  the  redesigned 
mount  (on  the  right)  are  shown  in  Fig.  5. 


FIG.  5.  Comparisons  of  the  original  and  the  redesigned 
sensor  mounting  brackets. 
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IV.  SUMMARY 


An  analytical  model  of  a  neutrally  buoyant,  polyurethane-encapsulated  velocity 
sensor  was  developed  and  analyzed  to  understand  the  interrelationships  between  the 
system  components.  Based  on  the  results  obtained  from  this  model,  it  was  evident  that 
desirable  velocity  sensor  characteristics  could  be  achieved  if  the  entire  sensor  unit  was 
free  to  move  essentially  as  a  rigid  body.  This  requirement  allowed  constraints  on  the 
maximum  sensor  mounting  stiffness  to  be  generated.  In  addition,  another  constraint  on 
the  sensor  mounting  stiffness  was  generated  based  on  the  desire  to  have  the  mount  supply 
sufficient  vibration  isolation  so  that  the  sensor  was  ambient-noise  limited  rather  than  self¬ 
noise  limited.  These  mounting  stiffness  constraints  were  combined  to  form  an 
engineering  design  result  that  showed  the  mounting  stiffness  necessary  for  any  given 
weight  of  the  sensor  assembly.  Preliminary  estimates  using  feasible  sensor  assemble 
weights  identified  a  sensor  mounting  bracket  stiffness  of  100  lb/in.  A  mounting  bracket 
redesign  showed  how  the  existing  hydrophone  mounting  bracket  (800  lb/in.  stiffness) 
could  be  readily  modified  to  provide  the  necessary  value  of  stiffness. 
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ABSTRACT 

We  summarize  our  recent  analysis  and  developments  on  acoustic  pressure  gradient  sensors  of 
both  the  motion  and  fixed  types.  Acoustical  motion  types  generally  refer  to  sensors  whose 
outputs  are  proportional  to  acoustically  induced  motion  such  as  acceleration  (accelerometers), 
velocity  (velocimeters)  or  displacement  (seismometers).  We  also  report  on  pressure  gradient 
sensors  that  are  fixed  in  space.  One  standard  fixed  type  may  be  realized  by  using  pairs  of 
pressure  sensors  (hydrophones)  separated  in  space  and  operating  in  a  differential  mode  to 
determine  a  direct  pressure  gradient,  Ap/Az.  The  second,  less  common  approach,  is  based  on 
using  single  sensors  fixed  in  space  that  respond  to  the  pressure  differences  produced  by  the 
diffraction  path  length  introduced  by  the  sensor.  We  have  fabricated  a  number  of  these  devices 
including:  a)  the  Symmtree™,  a  compact  acoustical  motion  accelerometer  utilizing  piezoelectric 
plates  in  a  bimorph  cantilever  tree-like  fashion,  b)  the  Aronov-Cube™,  a  cubical  device  using 
flexural  bimorphs,  c)  fixed  pressure  gradient  sensors  of  the  differential  type  and  d)  single 
element  diffraction  type  transducers.  The  two  acoustic  (motion)  accelerometers  are  each 
packaged  in  cylindrical  near  neutrally  buoyant  aluminum  housings  having  attached 
omnidirectional  hydrophones  so  that  cardiod  beam  patterns  may  be  formed.  The  four  separate 
channels  (3-accelerometers,  1 -hydrophone)  have  the  same  co-located  acoustical  centers. 

I.  INTRODUCTION 

Much  attention  has  been  directed  in  recent  years  to  pressure  gradient  directional  sensors  of  the 
first  order  that  have  a  figure-of-eight  (dipole)  directivity  pattern.  For  some  reason,  the 
majority  of  the  development  has  been  focused  on  one  type  of  pressure  gradient  sensor,  namely 
the  acoustical  motion  sensor  (so-called  particle  velocity  sensor),  as  if  they  possessed  significant 
advantages  over  other  pressure  gradient  types.  In  designing  the  pressure-gradient  sensor  for 
particular  applications,  there  exists  common  and  unique  problems  that  can  arise  and  be  solved  in 
somewhat  different  ways  for  different  sensor  types.  Therefore,  we  first  make  a  brief  review  of 
the  pressure-gradient  transducer  types. 

II.  IDEAL  SENSORS 

An  analysis  of  sound  fields  requires  sensors  that  measure  both  the  sound  pressure  p  (a  scalar 
quantity)  and  the  acoustical  field  motion  which  is  most  often  characterized  by  the  acoustic 
“particle”  velocity  u  (a  vector  quantity).  However,  the  acoustical  motion  of  the  fluid  is  the 
result  of  the  action  of  a  gradient  in  the  sound  pressure  field  as  is  well  known  and  described  by 
the  linear  force  equation  w  =  -p0  ^ .  Therefore,  it  is  equivalent  to  measure  the  sound  pressure 

and  sound  pressure  gradient.  These  ideal  sensors  are  often  characterized  by  their  directivity 
patterns  or  receive  sensitivity  as  a  function  of  angle,  which  we  may  characterize  as  a  sinusoidal 
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function  of  order  n,  H(<|))  =  cosn((|)),  where  <\>  can  represent  an  azimuthal  angular  dependence. 

For  this  reason  a  pressure  sensor  (hydrophone)  is  called  a  “zero  (n  =  0)  order”  sensor,  an 
“omnidirectional”  sensor,  or  equivalently  a  “monopole”  sensor.  Likewise,  an  ideal  motion  or 
pressure  gradient  sensor  having  maximum  sensitivity  in  a  designed  direction  and  zero  sensitivity 
in  the  perpendicular  direction,  H(c|))  =  cos(c|)),  is  known  as  a  first  (n  =  1)  order  sensor,  “figure-of- 
eight”  directional  sensors,  or  equivalently  “dipole”  sensors. 

Ideal  Acoustic  Pressure  Sensor  (Zero  order):  An  ideal  acoustic  pressure  sensor,  or  zero  (0)  order 
sensor,  must  be  sufficiently  smaller  than  acoustical  wavelength  such  that  %d/\  « 1  where  d  is 

the  maximum  dimension  of  the  sensor.  If  the  sensor  is  not  small,  diffraction  effects  will  alter  the 
ideal  response.  An  ideal  pressure  sensor  has  an  omnidirectional  beam  pattern  response  and  an 
electrical  output  that  is  proportional  to  the  sound  pressure  only.  The  sensor  must  be  fixed  in  the 
sound  field  or  made  insensitive  to  motion.  An  illustration  of  an  ideal  sensor  is  provided  in 
Figure  la. 

Ideal  Pressure  Gradient  Acoustic  Motion  Sensor  PG-AMS  (First  order):  An  ideal  acoustic 
motion  sensor  can  be  imagined  as  a  small  body,  completely  free  to  move  in  a  sound  field  due  to 
the  action  of  a  force  arising  from  the  gradient  of  the  sound  pressure  field  whereby  the  small  body 
has  an  internal  means  to  detect  this  motion  such  as  having  a  small  accelerometer  or  velocimeter 
inside.  An  ideal  Acoustic  Motion  Sensor  has  a  first  order  (figure-of-eight)  directivity  pattern,  is 
only  sensitive  to  motion  in  a  designed  direction  and  insensitive  to  motion  in  a  perpendicular 
direction,  and  has  an  electrical  output  that  is  proportional  to  the  sound  pressure  gradient  or 
acoustical  acceleration  (velocity,  or  displacement).  The  ideal  device  must  be  insensitive  to  the 
actions  of  sound  pressure  and  the  possible  deformation  of  the  body  arising  from  this  action. 

Ideal  Pressure  Gradient-Fixed  Differential  Sensor  PG-FIX  (First  order):  An  ideal  pressure 
gradient  sensor  of  the  fixed  type  can  be  imagined  as  a  pair  of  small  (nd/'k«l)  ideal  and  identical 

acoustic  pressure  (zero  order)  sensors  separated  by  a  distance  d  also  small  compared  to  the 
acoustic  wavelength  and  electrically  connected  in  opposite  as  illustrated  in  Figure  lb.  The 
electrical  output  is  proportional  to  the  pressure  gradient  only  and  the  directivity  pattern  is  of  first 
order,  H(t(>)  =  cos1  (j). 

Ideal  Pressure  Gradient  -  Single  Element  Diffraction  Type  Sensor  PG-FIX  (First  order):  An 
ideal  single  element  diffraction  type  sensor  consists  of  a  single  element  whose  boundary  is  fixed 
in  the  sound  field  but  has  a  transduction  element  that  responds  due  to  the  diffraction  induced 
pressure  gradient.  The  electrical  output  is  proportional  to  the  pressure  gradient  only  and  the 
directivity  pattern  is  of  first  order. 

We  note  that  for  all  of  these  sensors,  since  small  is  defined  relative  to  the  acoustical  wavelength 
and  real  bodies  introduce  diffraction  effects  that  may  be  wavelength  dependent,  there  is  a  limited 
frequency  range  where  one  can  expect  a  near-ideal  response.  Also,  that  for  ideal  first  (1)  order 
sensors  of  the  motion  or  fixed  types  the  ideal  figure-of-eight  response  can  only  be  truly 
obtained  under  the  assumption  that  the  sensitivity  is  infinitely  large  because  the  sensitivity  “off- 
axis”  approaches  zero.  We  will  note  here  for  future  use  that  the  linear  force  equation  does 
predict  that  the  acoustical  particle  acceleration  of  a  wave,  jcou,  is  also  proportional  to  the 
pressure  gradient  (VP  =  -p  jam) . 
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Figure  la.  Illustration  of  an  ideal  monopole  (0  order)  sensor  and  characteristic  beam  pattern, 
b.  Illustration  of  an  ideal  dipole  (1  order)  sensor  and  characteristic  beam  pattern. 
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Figure  2.  Illustrations  of  effective  unwanted  actions. 
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III.  REAL  SENSORS,  UNWANTED  ACTIONS,  AND  THE  SYSTHESIS  OF  CARDIODS 


Real  acoustical  sensors  are  never  quite  “ideal”.  As  there  can  be  no  “point”  sensors  and 
sensitivity  is  proportional  to  the  size  of  the  element  even  ideal  omnidirectional  sensors  are  hard 
to  realize.  Also,  it  is  not  possible  to  truly  “fix”  a  sensor  in  space  without  adding  structures  of 
significant  size  that  will  alter  the  sound  field.  In  helping  to  understand  these  limits  it  is  useful  to 
establish  the  response  to  ideal  actions  and  to  unwanted  actions.  Some  of  the  effects  of  unwanted 
noise  are  illustrated  in  Figure  2. 

Combined  outputs:  The  interest  in  dipole  sensors  comes  largely  from  the  fact,  that  the 
combination  of  the  output  of  an  ideal  dipole  (pressure  gradient)  sensor  with  the  output  from  an 
ideal  omnidirectional  (pressure)  sensors  and  can  offer  a  compact  means  for  determining  the 
direction  of  an  acoustical  wave.  Such  a  combination  of  equal  weighted  outputs  produces  a  0+1 
or  cardiod  response  (valid  for  one  axis).  This  beam  pattern  is  illustrated  in  Figure  3a,  b  for  the 
case  when  the  outputs  of  ideal  sensors  of  equal  weight  are  combined  and  in  case  c)  when  there  is 
a  mismatch  in  amplitude. 


l+COS  *((()) 


dB  Scale 


Figure  3.  Illustration  of  an  ideal  cardiod  beam  pattern  response  derived  from  an  ideal  monopole 
(0  order)  sensor  and  dipole  (1  order)  sensor  on  a)  linear  scale  and  b)  dB  Scale.  It  is  seen  that  the 
deepness  of  Cardiod  null  is  arrived  at  by  the  matching  of  the  0  mode  and  1  mode  response,  and  it 
has  little  to  do  with  the  deepness  of  the  dipole  null.  In  c)  we  have  an  ideal  cardiod  plotted  (blue) 
and  the  case  of  a  mismatch  in  amplitude  matching  of  the  omni  and  dipole  of  ldB  (red)  and 
2dB(Green).  The  plots  are  all  re-normalized  for  comparisons. 
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IV.  COMPARISON  OF  PRESSURE  GRADIENT  ACOUSTIC  MOTION  SENSORS  (PG- 
AMS)  WITH  THE  PRESSURE  GRADIENT  FIXED  SENSORS  (PG-FIX)  OF  THE 
DIFFRACTION  TYPE 


Introduction:  Pressure  gradient  (PG)  sensors  fixed  in  space  may  be  a  preferred  alternative  to  the 
use  of  acoustic  motion  sensors  for  certain  acoustic  low-frequency  applications.  In  both  types,  the 
sensors  are  pressure  gradient  sensors,  and  the  output  approaches  zero  (or  the  noise  floor)  at  low 
frequencies.  It  is  well  known  that  a  pair  of  sensors  fixed  in  space  may  be  used  to  measure  the 
pressure  gradient  but  what  is  less  known  is  that  a  single  sensor,  fixed  in  space,  may  be  used  as 
a  pressure  gradient  device  by  exploiting  the  diffraction  caused  by  the  sensor  element  and  its 
support.  In  certain  applications  this  may  be  more  desirable  than  using  pressure  gradient  - 
acoustic  motion  sensors  (PG-AMS)  (whether  accelerometers,  velocimeters,  or  seismometers)  as 
these  devices  need  a  “free”  or  soft  suspension.  In  contrast,  a  pressure  gradient  sensor  can  be 
“fixed”  (PG-FIX)  rigidly  to  a  frame  or  structure.  We  will  show  that  motion  and  fixed  sensors  of 
the  same  size  will  have  nearly  the  same  sensitivity.  There  is  no  intrinsic  advantage  in  using 
motion  sensors  over  pressure  gradient  sensors  to  offset  this  reality.  Therefore,  we  suggest  our 
preference  for  the  fixed  type  sensors  for  the  applications  of  Low  Frequency  Bow  Arrays  (LFBA) 
and  Integrated  Bow  Conformal  Arrays  (IBCA)  as  the  sensors  can  be  fixed  to  the  structure.  An 
analysis  follows. 

Consider  two  identical  circular  piezoelectric  (bimorph)  disks.  One  is  configured  in  a  near 
neutrally  buoyant  case  and  used  below  resonance  (accelerometer);  the  second  sensor  will  employ 
the  same  disk  supported  at  its  periphery  but  fixed  in  space  and  used  below  resonance.  These  two 
cases  are  shown  in  the  Figure  1.  Under  action  due  to  the  acceleration  of  the  case,  the  distributed 
uniform  mechanical  force  is 

f  =  ptu  (1) 

The  sound  pressure  action  P  may  be  assumed  to  act  uniformly  over  the  face.  Therefore,  we  can 
write  the  specific  acoustic  accelerometer  sensitivity  yu  in  terms  of  the  pressure  sensitivity  yp  as 

Yu=P?-Yp  (2) 


2a 


Figure  4.  Illustration  of  the  same  supported  disk  used  as  a  pressure  gradient  device  in 
a)  a  PG  Acoustic  Motion  Sensor  (PG-AMS)  mounted  in  a  moveable  case  and 
b)  a  PG  sensor  of  the  Fixed  (PG-FIX)  and  diffraction  type. 
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We  will  assume  the  sphere  vibrates  in  a  plane  wave  sound  field  with  velocity  u  according  to 

u=  P/pc.  (3) 

In  the  first  case  of  the  pressure  gradient  acoustic  motion  sensor  (PG-AMS)  (denoted  A),  we  have 
the  voltage  output  and  sensitivity  of  in  terms  of  pressure 
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CO 


=  y - P  * 

output  l  u  ’ 

1  pc 
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PnC 


(4) 


In  the  second  case  of  the  pressure  gradient  fixed  (PG-FIX)  diffraction  type  (denoted  B), 
we  have  the  voltage  output  and  sensitivity  in  terms  of  pressure  as 

V  =y  -2 ^--ka-P  (5) 


output 
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3k 


y®=y„  —  ka-na2  -0.45  , 

3k 

where  (8/3tc)  ka  is  the  classical  diffraction  coefficient  responsible  for  the  differential  path 
length.  Thus,  comparing  the  pressure  sensitivity  of  the  two  cases  we  have 


(6) 
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where  the  ratio  of  the  density  of  the  piezoelectric  ceramic  to  water  is  approximately  p  /p  =  7.7 . 

Thus,  we  conclude  that  a  PG-FIX  diffraction  type  device  having  relative  dimensions  of 
radius  of  the  outer  support  fixture  (a)  to  the  thickness  of  the  disk  (t)  on  the  order  of  10,  will  have 
a  sensitivity  that  is  equal  to  the  sensitivity  of  the  acoustic  motion  sensor  PG-AMS  type.  This 
may  be  expressed  (for  frequencies  below  resonance)  as  follows: 

Yp=Yp’~;  if- =  10,  thenYp=Yp  .  (8) 

at 

If  a  PG-FIX  sensor  can  produce  comparable  sensitivities  as  PG-AMS,  the  practical  matters  of 
vibration  isolation  and  mounting  may  be  quite  advantageous.  Of  course  if  a/t  is  made  larger  than 
10,  additional  advantage  is  achieved.  Also,  for  the  PG-AMS  the  case  has  to  made  to  withstand 
hydrostatic  pressure,  while  there  is  no  such  requirement  for  the  PG-FIX  single  element  diffraction 
type  sensor. 


V.  REAL  EXAMPLES:  PRACTICAL  DEVICES 

Aronov  has  developed  and  prototyped  a  number  of  practical  multiaxis  pressure  gradient  sensors 
some  of  which  are  highlighted  here  and/or  in  the  accompanying  presentation.  A  detailed 
summary  is  beyond  the  scope  of  this  summary  paper.  Generally,  for  low  frequency  applications, 
and  for  this  context,  we  will  define  as  those  frequencies  below  10  kHz,  it  is  desirable  to  use 
transduction  elements  that  are  flexural  devices  most  commonly  of  the  piezoelectric  bimorph 
approach.  These  are  illustrated  in  the  figure  below. 
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Figure  5.  Illustration  of  flexural  transducer  elements  used  for  Pressure  Gra’dient  Sensors  of  both 

motion  and  fixed  types. 


COMPACT  SYMMTREE™  CANTILEVER  PG-AMS: 


Aronov  has  developed  a  cantilever  based  multiaxis  PG-AMS  that  employs  symmetrically 
mounted  transduction  elements  on  a  common  tree-like  base  hence  the  name  Symmtree™.  Two 
proof-of-concept  devices  have  been  fabricated  using  conventional  transduction  materials  (PZT- 
4)  and  patents  are  pending.  The  devices  also  have  integrated  hydrophones  in  the  form  of  flexural 
disks  and  the  attributes  that  all  four  channels  (3  accelerometers  and  1  hydrophone)  have  the  same 
acoustical  center.  The  device,  packaged  in  a  neutrally  buoyant  housing,  is  compact  having 
approximate  dimensions  of  height  2.9  in  (7.5  cm)  and  diameter  2.375  (5.8cm).  With  PZT-4 
plates  device  BT-Syml  has  a  sensitivity  of  y  =  80  mV/g,  capacitance  1.8  nF  ,  and  resonance 
frequency  of  5.5  kHz.  A  photo  of  a  Symmtree™  accelerometer  is  illustrated  in  the  figure  below. 


Figure  6.  Photos  of  a  BTECH  Symmtree™  PG-AMS.  Patents  pending. 


ARONOV  CUBE  PG-AMS  (ACOUSTICAL  ACCELEROMETER): 

Aronov  has  also  developed  a  flexural  bimorph  multiaxis  PG-AMS  that  is  compact  and  has 
demonstrated  high  sensitivity  (377  mV/g).  The  device  has  integrated  hydrophones  in  the  form  of 
flexural  devices  and  has  the  attributes  that  all  four  channels  have  the  same  acoustical  center.  The 
device,  packaged  in  a  neutrally  buoyant  cylindrical  housing,  has  approximate  dimensions  of 
height  2.9  in  (7.5  mm)  and  diameter  2.375  (5.8  mm).  With  PZT-4  plates  device  BTECH-Cube-01 
has  a  channel  sensitivity  of  462  mV/g  @  1kHz  and  a  resonance  frequency  of  3.5  kHz.  A  sensor 
illustration  is  shown  in  Figure  7.  The  capacitance  is  1.2  nF  per  channel.  Beam  patterns  of  a  tested 
device  are  shown  in  Figure  8.  Modifications  of  this  design  have  produced  a  sensor  with  a 
sensitivity  of  1250  mV/g  at  1  kHz  as  shown  in  Figure  9. 
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Figure  7.  Illustration  of  a  BTECH  ARONOV-CUBE™  PG-AMS.  Patents  pending 
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Figure  9.  Frequency  /  sensitivity  response  of  the  Aronov  Cube  with  different  damping 
treatments.  Data  shown  for  one  plate  only.  Two  plates  per  channel  produce  a  sensitivity  of 

1250  mV/g  @  1  kHz. 

Regarding  Quality  Factor,  Resonance  Frequency  and  Damping: 

Often  the  quality  factor  (Q)  needs  to  be  sufficiently  low  in  a  receiver  in  order  to  prevent  the 
device  from  ringing  and  to  extend  the  useable  frequency  range  particularly  when  that  the 
resonance  of  the  transduction  device  is  near  to  the  frequency  band  of  interest  or  operation. 
Increasing  the  receiver  sensitivity  while  maintaining  a  high  resonance  frequency  is  difficult,  as  it 
can  be  shown  that  the  two  are  inversely  related,  or  more  exactly  the  product  of  the  sensitivity  and 
square  of  the  resonance  frequency  is  often  a  configuration  design  constant.  The  Aronov  Cube 
Accelerometer  device  affords  means  for  the  convenient  introduction  of  damping  materials  and 
the  application  specific  optimisation  of  this  transducer  is  under  development.  Figure  9  illustrates 
the  frequency  response  of  the  transducer  with  the  application  of  damping  materials. 
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YI.  PRESSURE  GRADIENT  FIXED  DIFFERENTIAL  AND  DIFFRACTION  TYPES 

A  number  of  Pressure  Gradient  Fixed  Differential  and  Diffraction  Type  Transducers  are  under 
development  and  will  be  discussed  at  the  meeting  including  a  three-axis  accelerometer  in  a 
planar  package 
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Research  (ONR  321 SS),  internal  funding  (BTECH)  for  prototype  development,  and  industry 
contributions  from  EDO  Ceramics  for  PZT  materials  and  testing. 


9 


Advances  in  Acoustic  Particle 
Velocity  Sensorics 

Paul  A.  Wlodkowski 
Fred  Schloss 
Wilcoxon  Research,  Inc. 

21  Firstfield  Road 
Gaithersburg,  Maryland  20878 


Advances  in  Acoustic  Particle  Velocity  Sensorics 


Paul  A.  Wlodkowski  and  Fred  Schloss 
Wilcoxon  Research,  Inc. 

21  Firstfield  Road 
Gaithersburg,  Maryland  20878 


Abstract.  The  development  of  a  high-sensitivity,  low-noise  acoustic  particle 
velocity  sensor  in  a  small  package  remained  a  formidable  technical  challenge  for 
many  years.  With  the  advent  of  single-crystal  piezoelectric  material,  however, 
that  goal  has  been  achieved.  Wilcoxon  Research  has  successfully  developed  a 
Vector  sensor  that  meets  U.S.  Navy  requirements  through  a  synergy  of 
technology  encompassing  single  piezoelectric  crystals,  a  novel  amplifier  design 
with  low-noise  characteristics,  and  a  unique  method  of  inducing  strain 
amplification  in  the  piezoelectric  crystal. 


I.  INTRODUCTION 

Acoustic  particle  velocity  transducers,  also  known  as  Vector  Sensors,  have  a  history 
rich  in  innovation.  During  World  War  II,  Bell  Laboratories1  developed  a  pressure  gradient 
hydrophone,  comprised  of  electromagnetic  inertia-type  transducers  mounted  in  a  spherical 
aluminum  shell.  Each  orthogonal  pair  of  sensors  has  a  figure  eight  or  cosine  pattern.2 
T.  J.  Schultz  of  the  Acoustics  Laboratory  at  Harvard  University  constructed  an  acoustic 
wattmeter3  for  measurements  in  air.  Concurrently,  an  effort  was  undertaken  at  the  Naval 
Ordnance  Laboratory4  in  which  the  effect  of  buoyancy  on  the  particle  velocity  sensor  was 
discussed.  At  the  Diamond  Ordnance  Fuse  Laboratory,5  an  ingenious  method  was 
established  to  relate  the  sound  intensity  to  temperature  fluctuations  in  air.  G.  L.  Boyer6  first 
incorporated  piezoelectric  materials  into  an  acoustic  intensity  meter  at  the  David  Taylor 
Model  Basin.  F.  W.  Desiderati7  at  the  Naval  Ship  Research  and  Development  Center 
demonstrated  a  bona  fide  pressure  gradient  transducer,  consisting  of  two  pressure-sensitive 
elements  separated  by  a  fixed  distance.  Moreover,  a  technique  for  correcting  near-field 
effects  was  provided  using  electrical  compensation  in  which  an  RC  differentiating  circuit 
was  placed  at  the  output  of  the  pressure  gradient,  or  acceleration  channels,  of  the 
respective  transducers. 

The  dormancy  in  this  field  for  the  next  thirty  years  can  be  attributed  to  limitations  in 
signal-processing  capability.  Once  signal-processing  technology  evolved  to  a  mature 
level,  further  requirements  were  placed  on  the  sensor  design.  The  development  of  a  high- 
sensitivity,  low-noise  acoustic  particle  velocity  sensor  in  a  small  package  presented  a 
formidable  technical  challenge  to  researchers  for  many  years.  By  harnessing  the 
exceptional  electromechanical  properties  of  single  piezoelectric  crystals8,9  with  a  novel 
amplifier  with  low-noise  characteristics,  and  a  unique  method  of  strain  amplification  in 
the  piezoelectric  crystal,  that  elusive  goal  was  achieved. 


1 


II.  DESIGN  CONCEPT  AND  PROTOTYPE 


Having  recognized  that  the  exceptional  properties  of  single  piezoelectric  crystals, 
belonging  to  the  relaxor  group  of  ferroelectric  materials,  would  be  ideally  suited  for 
developing  low-noise  vibration  transducers,  engineers  at  Wilcoxon  Research,  Inc.  have 
applied  a  method  of  strain  amplification10  to  enhance  the  performance.  The  method  of 
strain  amplification  utilizes  a  sensing  structure  comprised  of  a  bending  plate  that  has  a 
castellated  surface,  which  raises  the  piezoelectric  sensing  material  that  is  mounted  there 
above  the  substrate’s  neutral  axis.  Strain  in  the  crystal  is  amplified  via  mechanical 
leverage.  Resultantly,  stress  levels  in  the  piezoelectric  materials  are  augmented, 
generating  more  charge  output  and  lower  noise  than  traditional  designs.  Without  the 
castellated  surface,  the  configuration  would  be  that  of  a  conventional,  flexural  disc  mode 
accelerometer.11 

The  castellated  height  is  adjusted  depending  upon  the  electromechanical  properties 
of  the  sensing  material,  and  it  amplifies  the  stress  in  the  crystal  and  generates  a  higher 
charge  output.  Application  of  this  technique  results  in  the  increase  of  the  charge  output 
for  the  transducer  by  a  factor  of  two  (6  dB),  irrespective  of  the  piezoelectric  material 
chosen.  Most  significantly,  this  is  accomplished  without  adding  inertial  mass  to  the 
sensor  and  by  keeping  the  resonance  frequency  constant.  A  generalized  approach  to  the 
method  of  strain  amplification,  incorporating  single  piezoelectric  crystals,  has  been 
published  recently  by  Wilcoxon  Research,  Inc.12 

Wilcoxon  Research,  Inc.  initiated  development  on  its  single  piezoelectric  crystal- 
based  acoustic  particle  velocity  sensor  in  October  2000.  Five  months  later,  in  March  2001, 
the  first  prototype  was  assembled,  and  measurements  were  performed  at  the  Acoustic  Test 
Facility  (ATF)  at  the  Naval  Undersea  Warfare  Center  (NUWC)  in  Newport,  Rhode  Island. 
A  photograph  of  the  Wilcoxon  Research  Vector  sensor  is  illustrated  below  in  Fig.  1. 


III.  EXPERIMENTAL  SETUP  AND  MEASURED  RESULTS 


Measurements  on  the  Wilcoxon  Research  Prototype  01  Vector  Sensor  were  taken  in 
accordance  with  ANSI  SI. 20- 1988  (Revision  of  SI. 20- 1972).  The  water  temperature 
was  17°C.  Projector  J-9  #125  and  hydrophone  H-52  #101  were  used  to  meet  USRD 
standards.  Receive  voltage  sensitivities  and  phase  on  all  four  channels  were  recorded. 
Beam  patterns  were  generated  in  the  XY  plane  for  the  X  and  Y  accelerometer  channels 
and  the  pressure  transducer;  and  in  the  XZ  plane  for  the  X  and  Z  accelerometer  channels 
and  the  pressure  transducer.  The  distance  between  the  projector  and  Vector  sensor  was  a 
constant  1.62  meters  over  the  swept  frequency  range  of  200  Hz  to  2  kHz. 

The  received  voltage  sensitivities  for  the  four  channels  of  the  Vector  sensor  are 
illustrated  below  in  Fig.  2.  . 
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FIG.  1.  Vector  sensor  (Wilcoxon  Research  Prototype  01). 
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FIG.  2.  Receive  voltage  sensitivities  (Wilcoxon  Research  Prototype  01). 
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As  evident  from  the  graph  in  Fig.  2,  the  receive  voltage  sensitivities  of  the 
accelerometer  channels  differ  by  only  a  fraction  of  1  dB,  and  the  slope  of  the 
accelerometer  sensitivity  curve  follows  theory  according  to  Newton’s  law  for  water 
motion. 

Beam  patterns  in  the  XZ  plane  are  shown  below  in  Fig.  3.  A  typical  cosine  pattern 
was  generated,  and  most  notably  nulls  of  greater  than  -35  dB  were  obtained.  This  result 
translates  into  a  transverse  sensitivity  of  less  than  2%. 
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FIG.  3.  Beam  pattern  in  XZ  plane  (Wilcoxon  Research  Prototype  01). 


The  phase  response  between  the  Y  channel  of  the  accelerometer  and  the  pressure 
transducer  is  depicted  in  Fig.  4. 
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FIG.  4.  Phase  response  of  Wilcoxon  Research  Vector  Sensor  Prototype  01. 


As  the  graph  indicates,  the  phase  response  is  90°±4°,  as  specified  by  the  U.S. 
Navy  in  its  design  objectives.  Given  the  measurement  limitations  of  the  Acoustic  Test 
Facility,  it  is  expected  that  the  phase  response  will  be  significantly  improved  in  a  free 
field  with  a  continuous  wave  source.  Low  frequency  phase  variations  may  be  caused  by 
output  cable  restraint.  Moreover,  a  focused  effort  is  underway  at  Wilcoxon  Research  to 
enhance  the  transducer  performance,  which  will  further  mitigate  phase  between  the 
channels. 

Evaluation  of  the  inherent  noise  levels  of  acoustic  particle  velocity  sensors  is 
critical  for  a  robust  design.13  It  has  also  been  shown  that  accelerometer  noise  is  inversely 
proportional  to  charge  output,  especially  at  low  frequencies.14  The  measured 
accelerometer  noise  of  the  Vector  sensor,  converted  to  the  equivalent  pressure  level,  is 
shown  in  Fig.  5  as  a  function  of  frequency.  Superimposed  on  the  graph  is  the  Knudsen 
Sea  State  Zero  level,  extrapolated  to  50  Hz.  If  necessary,  it  will  be  possible  to  lower  the 
noise  floor  further  through  a  combination  of  electrical  and  mechanical  design  efforts. 
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FIG.  5.  Measured  accelerometer  noise 
(Wilcoxon  Research  Vector  Sensor  Prototype  01). 


IV.  CONCLUSIONS 

An  advance  in  acoustic  particle  velocity  sensorics  has  been  achieved  with  the 
development  of  a  Vector  sensor,  incorporating  single  piezoelectric  crystals,  a  novel 
amplifier  with  low-noise  characteristics,  and  a  unique  method  of  inducing  strain 
amplification  in  the  piezoelectric  crystal.  Directivity  patterns,  phase  response,  receive 
voltage  sensitivities,  and  the  equivalent  pressure  level  of  accelerometer  noise  demonstrate 
that  U.S.  Navy  requirements  have  been  met.  This  is  evidenced  by  a  Vector  sensor 
exhibiting  deep  nulls  exceeding  35  dB  and  phase  response  between  hydrophone  and 
accelerometer  channels  of  90°±4°.  It  is  anticipated  that  further  improvements  to  the 
performance  can  be  expected  from  enhancements  to  transducer  and  amplifier  design. 
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Abstract.  EDO  has  developed  acoustical  directional  sensors  that  address  the 
needs  of  two  U.S.  Navy  hull  arrays,  Conformal  Acoustic  Velocity  Sonar 
(CAVES)  and  Low  Frequency  Hull  Array  (LFHA).  These  sensors  are  based 
upon  a  flexural  disc  accelerometer  design.  The  CAVES  sensor  is  a  uni-axial, 
neutrally  buoyant  accelerometer  design.  The  LFHA  sensor  consists  of  tri-axial 
accelerometers  and  pressure  sensors.  This  “dyadic”  configuration  is  capable  of 
generating  a  directivity  index  of  9.5  dB.  The  MR  A  can  be  steered  in  any 
direction  by  adjusting  the  channel  weights.  The  designs  have  been  tested 
extensively  at  U.S.  Navy  test  facilities.  The  accelerometer  design  passes  the 
MIL-90 1-C  explosive  shock  test. 

I.  DEVELOPMENT  HISTORY 

The  development  of  directional  acoustic  sensor  technology  at  EDO  was  driven  by 
the  U.S.  Navy  need  for  hull  mounted  motion  sensors.  A  chronological  history  of  events 
is  provided  in  Fig.  1.  A  uniaxial  accelerometer  design  was  utilized  in  both  the  Conformal 
Acoustic  Velocity  Sonar  (CAVES)  and  CAVESAVide  Aperture  Array  (WAA)  programs. 
An  outgrowth  of  this  design  was  later  implemented  into  vector  and  dyadic  sensor  designs. 


1995  Motion  Sensor  Workshop 

-  EDO  Paper:  Concept  for  a  Low  Profile  Mold-in-Place  Accelerometer 

1996  CAVES  Spherical  Neutral  Buoyant  Accelerometer 

-  Flexural  Disc  Design 

1997/1998  CAVES  Cylindrical  Neutral  Buoyant 
Accelerometer 

-  Flexural  Disc  Design,  Syntactic  Foam 

1998  Paper  Published 

-  NUWC/EDO  Collaboration 

-  Moffett,  Trivett,  Klippel,  Baird,  A  Piezoelectric,  Flexural-Disc,  Neutral 

Buoyant,  Underwater  Accelerometer,  IEEE  Transaction  on  Ultrasonics, 
Ferroelectrics,  and  Frequency  Control,  Vol.  45,  No.  5,  1998 

1998  CAVES/ ASTO  Patch  Test  Accelerometers 

-  Delivered  120  Accelerometer  Assemblies 

1998-Present  LF  Bow  Array  Vector  Sensor  Development 

-  September  2000:  NUWC  Test/Technology  Demonstration 

-  May  2001 :  NUWC  Dyadic  Sensor  Test  -  Seneca  Lake 

2000-Present  CAVES  WAA  Preliminary  Design 

-  Subcontractor  to  General  Dynamics,  Electric  Boat _ 


FIG.  1.  EDO  Motion  Sensor  History 
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II.  DIRECTIONAL  ACOUSTIC  SENSOR  DESIGN  OBJECTIVES 


The  development  of  a  directional  acoustic  sensor  is  driven  by  applications  where  an 
artificial  baffle  is  required.  In  these  applications,  the  use  of  a  passive  resistive  baffle  is 
either  impracticable  or  negatively  impacts  the  performance  of  other  sonar  systems. 


One  such  application  would  be  the  Low  Frequency  Hull  Array  (LFHA)  as  illustrated 
in  Fig.  2.  The  LFHA  is  a  horseshoe-shaped  array  of  sensor  staves  that  is  conformal  to  the 
bow  of  a  688-class  submarine.  The  LFHA  must  be  acoustically  trans-parent  to  the  other 
sonar  arrays  that  reside  behind  it.  Currently,  the  LFHA  utilizes  the  DT-276,  an  omni¬ 
directional  pressure  sensor.  Processing  against  ownship  self-noise  can  occur  only  at  the 
array  level.  As  indicated  in  Fig.  2,  many  staves  are  located  within  the  ballast  tanks.  As 
a  result,  the  self-noise  field  varies  with  position.  The  use  of  directional  sensors  in  this 
application  would  permit  the  processing  against  local  self-noise  at  each  sensor  location. 


Stave 


FIG.  2.  Low  Frequency  Hull  Array  Configuration 

(Courtesy  NUWC  Division  Newport) 

The  objective  of  this  development  is  to  design  a  sensor  that  will  be  capable  of 
generating  a  baffle  that  can  be  steered  in  a  desired  direction.  The  sensor  should  also  be 
capable  of  generating  cardioid  directivity  patterns  at  a  minimum.  The  form  factor  of  the 
sensor  shall  also  be  constrained  to  the  existing  size  of  the  DT-276.  Additional  require¬ 
ments  are  detailed  within  a  NUWC  Division  Newport  motion  sensor  specification. 

III.  UNDERLYING  THEORY 

The  dyadic  or  tensor  sensor  outputs  three  quantities:  pressure,  acceleration 
(proportional  to  gradient  of  pressure),  and  the  gradient  of  acceleration.  A  scalar  sensor  or 
hydrophone  utilizes  only  the  pressure  quantity.  A  vector  sensor  utilizes  the  pressure  and 
the  acceleration,  whereas  the  dyadic  utilizes  all  three  quantities.  The  directivity  index  of 
each  of  these  sensor  configurations  is  provided  in  Fig.  3. 


Sensor  Type 

Tensor  Rank 

Directivity  Index  (dB) 

Scalar 

0 

0 

Vector 

1 

4.8-6. 0 

Dyadic 

2 

7. 0-9. 5 

FIG.  3.  Sensor  Directivity 
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The  vector  sensor  can  generate  a  standard  cardioid  directivity  pattern  with  one 
null,  or  by  selecting  the  weights,  generate  higher  directivity  index  with  two  nulls.  The 
dyadic  sensor  due  to  its  spatial  sensor  offset  can  generate  a  higher  directivity  index.  In 
order  to  realize  the  performance  improvement  at  low  frequencies,  dyadic  sensors  must  be 
very  closely  matched  in  sensitivity. 

IV.  DYADIC  SENSOR  DESIGN 

The  dyadic  sensor  contains  both  accelerometers  and  pressure  sensors.  The  accelero¬ 
meters  are  arranged  to  provide  orthogonal  outputs.  The  accelerometers  are  both  co-located 
and  distributed  relative  to  the  acoustic  center  of  the  pressure  sensor.  This  configuration 
(see  Fig.  4)  provides  the  capability  to  generate  scalar,  vector  and  dyadic  outputs. 


mzK-i 


FIG.  4.  Dyadic  Sensor  Design 


V.  ACCELEROMETER  DESIGN 

The  accelerometer  design,  as  illustrated  in  Fig.  5,  is  composed  of  an  inertial  mass 
substrate,  two  ceramic  discs,  and  a  pressure-tolerant  housing.  The  flexural  disc  design  is 
electrically  terminated  to  provide  voltage  doubling. 


FIG.  5.  Accelerometer  Design 
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VI.  ACCELEROMETER  PERFORMANCE 


Accelerometer  performance  parameters  that  have  been  measured  include  resonance 
frequency,  capacitance,  electrical  dissipation,  and  sensitivity  to  motion  and  pressure.  The 
nominal  characteristics  of  the  sensor  without  cable  are 

Capacitance:  3400  pF 

Electrical  Dissipation:  0.014 

Resonance  Frequency:  11.2  kHz 

Acceleration  Sensitivity:  -17.25  dB 

The  accelerometer’s  insensitivity  to  pressure  over  the  operating  range  of  hydrostatic 
pressure  has  been  demonstrated  using  NUWC’s  System  J.  The  motion  sensitivity  has 
been  demonstrated  using  shaker  head  and  in-water  measurements.  Shaker  head 
performance  at  NUWC  with  30  feet  of  cable  is  illustrated  in  Fig.  6. 


FIG.  6.  Shaker  Head  Acceleration  Sensitivity 

In  support  of  the  CAVES  Program,  EDO  manufactured  120  accelerometer 
assemblies.  These  sensors  were  manufactured  without  binning  of  components.  The 
manufacturing  process  capability  associated  with  that  product  is  provided  in  Fig.  7. 
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Capacitance  (nF) 

Standard  Deviation  =  0.063 


Resonance  Frequency  (kHz) 

Standard  Deviation  =  0.072 


Sensitivity  (dB//V/g) 

Standard  Deviation  =  0.19 


FIG.  7.  Manufacturing  Performance 


Through  the  use  of  component  selection,  using  characteristics  such  as  substrate  mass 
and  ceramic,  the  stated  performance  standard  deviations  can  be  reduced  by  a  factor  of  3:1. 


VII.  DYADIC  SENSOR  PERFORMANCE 
Directivity 

The  directivity  performance  of  the  dyadic  sensor  at  three  frequencies  (200,  1000, 
and  5000  Hz)  is  provided  in  Figs.  8  through  10.  The  measurements  at  200  and  1000  Hz 
were  made  at  the  NUWC  ATF  at  a  distance  of  1.55  m.  The  measurements  at  5000  Hz 
were  gathered  at  the  EDO  ATF  at  6  m.  The  sensor  in  all  cases  was  attached  to  the 
mounting  fixture  illustrated  in  Fig.  4. 
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FIG.  8.  Dyadic  Sensor  Directivity  at  200  Hz 
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FIG.  9.  Dyadic  Sensor  Directivity  at  1000  Hz 
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FIG.  10:  Dyadic  Sensor  Directivity  at  5000  Hz 

Acoustic  scattering  from  surfaces  in  close  proximity  to  the  sensor  will  impact  the 
directivity  performance.  The  LFHA  mounting  fixture  will  impact  the  directivity 
performance  at  higher  frequencies.  An  example  of  this  impact  upon  directivity  is 
provided  in  Fig.  11.  The  performance  without  the  LFHA  mount  is  ideal  dipole. 


With  Mount  Without  Mount 


FIG.  11.  Dyadic  Sensor  Directivity  at  5000  Hz,  X-Axis, 
Mounted  With  and  Without  LFHA  Mount 
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Receive  Sensitivity 


The  nearfield  receive  sensitivity  of  the  X,  Y,  and  P  channels  is  provided  in  Fig.  12. 


100  1000  10000 

Frequency  (Hz) 

FIG.  12.  Dyadic  Sensor  Receive  Sensitivity,  Horizontal  Plane 


The  sensitivity  performance  below  1000  Hz  does  not  fall  off  at  6  dB  per  octave 
as  predicted  by  theory.  The  higher  measured  sensitivity  is  the  result  of  the  nearfield 
position  of  the  sound  source.  The  separation  between  sensor  and  sound  source  was 
1.55  m.  The  nearfield  particle  velocity  for  spherical  wave  is  given  by 


u  = - cos(ft)t  —  kr — (p)  • 

per  kr 

When  corrected  for  the  nearfield  effect  at  200  Hz,  the  difference  between  the  measured 
and  predicted  is  0.4  dB. 


Electronic  Noise 

The  predicted  electronic  noise  performance  of  dyadic  sensor  is  provided  in  Fig.  13. 
The  noise  model  used  has  been  validated  by  measurements.  As  illustrated  in  the  inset 
graph,  the  measured  values  (blue)  and  model  predictions  (red)  are  in  close  agreement. 
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The  design  provides  compliant  performance  relative  to  the  maximum  acceptable  level  of 
noise  specified  by  the  NUWC  requirement.  At  the  lower  frequencies  the  Johnson  noise 
of  the  sensor  dominates  the  performance. 


50  250  450  650  850  1050  1250  1450  1650  1850 

Frequency (Hz) 

FIG.  13.  Dyadic  Sensor  Electronic  Noise  versus  Frequency 


VIII.  ENVIRONMENTAL  PERFORMANCE 

The  accelerometer  utilized  in  the  dyadic  sensor  design  has  undergone  extensive 
environmental  testing.  This  testing  includes  hydrostatic  pressure,  non-operating  pressure, 
and  underwater  explosive  shock.  Figure  14  defines  the  conditions  under  which  the 
accelerometer  has  been  demonstrated  to  operate  and/or  survive. 


Operating  Pressure  Range 

•  0-1000  psi 

Non-Operating  Temperature 
•-100°  to  250°  F 

Explosive  Shock,  MIL-901  C,  Grade  A 

•  Pressure  Release  Mount 


FIG.  14.  Accelerometer  Environmental  Performance 
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IX.  CONCLUSIONS 


EDO  has  demonstrated  directional  motion  sensors  based  on  a  common  building 
block  that  can  be  produced  in  a  manufacturing  environment,  and  whose  performance 
addresses  both  CAVES  and  LFHA  motion  sensor  requirements.  The  EDO  dyadic  sensor 
supports  vector  sensor  processing.  The  demonstrated  performance  is  summarized  in 
Fig.  15. 


•  Acceptable  Dipole  and  Omni  Directivity 

•  Acceptable  Sensitivity  Matching 

•  Extended  Bandwidth  50-6000  Hz  (50-2000  Hz  Navy  Baseline) 

•  Acceptable  Electronic  Noise  Performance 

•  Accelerometer  Insensitivity  to  Pressure 

•  Acceptable  Hydrostatic  Pressure  Performance 

•  Broad  Operating/Non  Operating/Process  Temperatures 

•  Accelerometer  Passes  MIL- 901  C  Explosive  Shock  Test,  Grade  A 

FIG.  15.  Dyadic  Sensor  Demonstrated  Performance 

The  extended  bandwidth  of  the  EDO  dyadic  sensor  provides  the  potential  to  support 
target  ranging  using  triangulation  techniques  with  two  or  more  dyadic  sensors  in  an  array. 
This  feature  in  a  forward-looking  array  may  be  of  significant  interest  to  provide 
submarines  a  means  of  maintaining  tactical  advantage  in  close  encounter  situations. 

EDO  has,  in  addition,  demonstrated  a  capable  manufacturing  process.  The  process 
is  capable  of  producing  hardware  to  the  very  tight  tolerances  required  by  the  NUWC 
LFHA  Motion  Sensor  Specification. 

Lastly,  the  EDO  dyadic  sensor  provides  the  capability  to  support  a  dyadic  signal 
process  that  would,  in  theory,  provide  a  directivity  index  as  great  as  9.5  dB,  or  as  much  as 
3.7  dB  greater  than  a  vector  sensor. 
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Abstract.  The  use  of  directional  hydrophones  in  towed  systems  is  summarized.  The 
acoustic  and  flow  noise  performance  of  particle  velocity,  particle  acceleration,  and 
pressure  difference  or  gradient  on  hydrophones  are  compared.  Primary  noise 
mechanisms  and  fundamental  limits  for  directional  hydrophones  in  a  towed  array 
system  are  summarized.  Approaches  for  noise  mitigation  in  towed  systems  are 
presented. 


I.  INTRODUCTION 

Directional  hydrophones  are  those  that  intrinsically  have  a  non-axisymmetric 
directivity  pattern.  Two  common  examples  are  dipole  and  cardioid  sensors,  both  of 
which  are  available  as  commercial  audio  equipment.  Dipole  sensors  respond  to  one 
component  of  a  vector  field  (e.g.  the  particle  velocity)  and  cardioid  sensors  have  a 
maximum  response  at  a  single  angle.  Examples  of  beampattems  from  omnidirectional, 
dipole,  and  two  types  of  cardioid  sensors  are  shown  in  Fig.  1 . 

Two  salient  benefits  of  using  directional  hydrophones  are  the  ability  to  determine 
the  arrival  angle  of  an  acoustic  source  and  a  directivity  index  greater  than  0  dB.  A  linear 
array  of  directional  hydrophones  can  be  used  in  the  same  manner  as  conventional 
omnidirectional  hydrophones  to  provide  additional  gain  against  ambient  noise. 
Volumetric  (e.g.  multi-line  towed  arrays)  can  effectively  create  the  performance  of 
directional  hydrophones  through  beamforming. 

The  Office  of  Naval  Research  (Code  321SS)  funded  the  SuperDirective  Array 
(SDA)  and  Multi-Mode  Directional  Towed  Array  (MMDTA)  projects  at  the  former 
Naval  Underwater  Systems  Center  and  the  Naval  Undersea  Warfare  Center  Division 
Newport.  The  Surface  Ship  Torpedo  Defense  (SSTD)  program  also  funded  the  SDA 
effort.  The  primary  objective  of  the  projects  was  to  design  and  demonstrate  a  tactical 
single-line  towed  array  that  used  directional  hydrophones.  The  capabilities  of  the  array 
included  left/right  bearing  discrimination  and  enhanced  array  gain.  From  this  research 
came  an  understanding  of  the  important  noise  mechanisms  for  directional  hydrophones 
and  ways  of  minimizing  the  non-acoustic  noise  to  maximize  performance. 
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FIG.  1.  Sensor  beampattern  examples. 


Other  countries  have  pursued  development  of  directional  towed  arrays  for  reception 
of  active  sonar  signals.  This  implementation  is  particularly  appealing  since  the  left/right 
position  of  a  target  can  be  identified  using  a  simple  omnidirectional  or  axis-symmetric 
source  and  a  single  line  array  of  directional  hydrophones. 

II.  DIRECTIONAL  HYDROPHONE  DESIGN  IN  TOWED  SYSTEMS 

The  directional  hydrophones  developed  and  tested  under  the  SDA  and  MMDTA 
projects  consisted  of  an  omnidirectional  hydrophone  combined  with  two  dipole 
hydrophones.  One  such  design  is  depicted  in  Fig.  2.  The  dipole  hydrophones  were 
oriented  to  measure  the  acoustic  particle  velocity  along  two  axes  (x  and  y)  perpendicular 
to  the  array  axis  (z). 

In  this  example,  the  dipole  hydrophones  consisted  of  geophones  encapsulated  in 
syntactic  foam  to  achieve  approximate  neutral  buoyancy.  The  geophones  had  a 
nominally  flat  response  of  23.6  V/(m/s)  above  a  highly  damped  (C,  =  0.33)  28-Hz 
resonance.  The  reader  is  referred  to  Abraham1  for  further  details  of  this  sensor  and  to 
Gabrielson  et  al.2  for  design  of  a  similar  sensor. 
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Cardioid  hydrophone  groups  can  be  formed  in  the  same  manner  as  for  conventional 
omnidirectional  hydrophones.  An  “extended  sensor”  or  “high  wavenumber  filter”  is 
simply  a  number  of  hydrophones  distributed  within  half  of  the  design  acoustic 
wavelength  (XJ2  =  c/2fd).  Each  group  is  wired  together  in  either  series  or  parallel  (or  a 
combination  of  both)  to  minimize  electrical  noise  and  maximize  sensitivity.  The  omni 
and  dipole  outputs  are  digitized  by  three  separate  analog  to  digital  converters  that  are  part 
of  the  array  telemetry  system.  Additionally,  at  least  one  roll  sensor  is  needed  to  provide  a 
reference  to  gravity.  Some  systems  integrate  the  roll  sensor  with  the  dipole  sensor 
outputs  and  only  digitize  the  dipole  response  steered  to  the  horizontal.  This  configuration 
preserves  the  left/right  bearing  discrimination  capability  but  loses  the  full  elevation  angle 
steering  capability. 

Several  physical  mounting  systems  have  been  used  for  directional  hydrophones  in 
towed  arrays.  These  designs  include  suspension  in  open-cell  foam,  rubber  lobe-mounts, 
and  suspension  along  internal  axial  strength  members.  These  are  illustrated  in  Fig.  3. 

The  open-cell  foam  provides  excellent  mechanical  isolation  due  to  its  low  stiffness  but 
the  sensors  may  rotate  with  respect  to  the  array  boot  leading  to  angular  orientation 
uncertainty.  The  lobe  mounts  contact  the  boot  directly,  and  thus  sensor  orientation  is  not 
an  issue  but  a  more  direct  coupling  of  the  wall  motion  will  exist  compared  to  the  foam 
mount.  The  axial  strength  members  provide  alignment  of  all  the  hydrophones  in  the 
array  module  at  the  cost  of  coupling  axial  motion  into  the  response.  The  use  of  axial 
strength  members  with  acceleration  canceling  omnidirectional  hydrophones  is  appealing 
because  the  dipole  sensors  as  oriented  in  Fig.  2  will  have  a  null  in  the  axial  direction. 

The  off-axis  sensitivity  of  accelerometers  and  geophones  is  typically  4%  (-28  dB)  or  less. 
None  of  the  mounting  systems  are  optimal,  and  sensor  orientation  and  ultimately  acoustic 
performance  must  be  weighed  against  self-noise  performance. 


Syntactic  foam 
encapsulation 


omnidirectional  sensor  dipole  sensor 

(e.g.  Lead  Titanate  or  Lead  (e.g.  Geospace  GS-14) 
Zirconate  Titanate) 


FIG.  2.  Example  of  directional  hydrophone  components. 
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FIG.  3.  Directional  hydrophone  mounting  schemes. 


III.  DIRECTIONAL  HYDROPHONE  ARRAY  ACOUSTIC  OPERATION 

By  combining  the  calibrated  pressure  outputs  from  the  omnidirectional  and  dipole  hydro¬ 
phones,  a  cardioid  response  at  any  elevation  angle,  \| /,  can  be  obtained.  This  is  expressed 
mathematically  as, 


D  ,  ,  Pomni  +  Px  COS(V)  +  P y  sin(\J/) 

PcardioidW  = - 2 - - •  (1) 

The  directivity  index  (DI)  for  a  single  dipole  hydrophone  is  4.8  dB.  If  a  cardioid 
response  is  formed  as  in  Eq.  (1),  then  it  will  also  have  a  DI  of  4.8.  Weighting  the  dipole 
response  by  a  factor  of  three  greater  than  the  omnidirectional  response  forms  a 
beampattern  with  a  DI  of  6.  However,  this  is  achieved  at  the  cost  of  losing  the  uni¬ 
directional  response  null. 

Apart  from  forming  the  cardioid  (or  dipole)  response,  array  processing  for  the 
directional  hydrophones  parallels  that  of  omnidirectional  hydrophones.  Conventional  and 
adaptive3  beamforming  methods  have  been  used  to  examine  the  noise  field  to  which 
omnidirectional  and  dipole  hydrophones  are  subject  inside  a  fluid-filled  towed  array 
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module.  The  MMDTA  project  emphasized  different  modes  in  which  a  directional 
hydrophone  array  could  be  operated:  conventional,  dipole-only,  or  cardioid.  For  target 
searches  near  broadside  (e.g.  +/-  45°),  the  dipole-only  mode  foregoes  the  left/right 
bearing  discrimination  but  gains  excellent  end-  and  forward-fire  noise  rejection. 
Additionally,  the  vertical  distribution  of  ambient  noise  from  wave  breaking4  generally  has 
maxima  at  angles  pointing  towards  the  surface  and  bottom,  so  the  horizontal  dipole 
response  can  achieve  a  DI  in  excess  of  4.8  on  the  group  level. 

The  two-axis  dipole  hydrophone  array  discussed  thus  far  has  the  obvious 
disadvantage  that  the  acoustic  particle  velocity  along  the  array  axis  is  not  measured.  One 
could  construct  a  three-axis  dipole  sensor  for  this  purpose,  and  rejection  of  own-ship 
noise  is  one  resulting  advantage  for  beams  steered  towards  end-fire. 

IV.  DIRECTIONAL  HYDROPHONE  SELF-NOISE  PERFORMANCE 

Self-noise  is  defined  as  any  signal  in  the  sensor  response  that  is  uncorrelated  with 
the  desired  acoustic  signal.  For  both  omnidirectional  and  dipole  hydrophones,  electronic 
noise,  cable  strum,  array  strum,  turbulent  boundary  layer  noise,  and  hose  wall  vibration 
all  contribute  to  self-noise.  The  electronic  noise  performance  of  dipole  hydrophone 
parallels  that  of  conventional  hydrophones  with  the  exception  of  moving  coil  geophones, 
which  have  a  low  output  impedance  unlike  most  piezoelectrics.  Proper  buffer  electronics 
must  be  selected  to  ensure  sensor  output  compatibility  with  the  array  telemetry  system. 
Additionally,  if  the  dipole  sensor  consists  of  a  motion  sensor  with  an  in-band  resonance, 
particular  care  must  be  taken  to  assure  that  the  resonance  does  not  electrically  saturate  the 
array  analog  to  digital  converter  during  towing  operations. 

The  transverse  component  of  cable  strum  that  is  not  damped  out  by  a  Vibration 
Isolation  Module  (VIM)  will  be  measured  by  dipole  hydrophones  as  noise.  Additionally 
if  the  array  acoustic  module  itself  is  not  neutrally  buoyant,  then  it  may  strum  at  low  tow 
speeds  and  further  contribute  to  dipole  hydrophone  self-noise. 

The  turbulent  boundary  layer  that  develops  over  an  array  contributes  significantly  to 
the  self-noise  levels  of  both  omnidirectional  and  dipole  hydrophones.  The  fluctuating 
wall  pressure  and  wall  shear  stress  excite  the  array  hose.  Evanescent  and  propagating 
waves  are  excited  in  the  array  that  contribute  to  sensor  self-noise.  Peloquin5  gives  the 
response  of  shell  with  inner  and  outer  fluid  loading.  Examples  of  the  ratio  of  outer 
pressure  (e.g.,  from  the  turbulent  boundary  layer)  to  inner  pressure  (what  the  sensor 
responds  to)  for  a  7-cm  outer-diameter  flexible  cylinder  are  shown  in  Fig.  4.  The  n  =  0 
response  is  the  fluid  or  breathing  wave.  The  n  =  1  response  is  a  bending  or  shear  wave. 
The  n  =  2  lobar  response  cuts  on  at  approximately  46  Hz  for  this  cylinder.  All  three  of 
these  waves  are  dispersive.  The  dipole  hydrophones  will  response  primarily  to  the 
bending  wave  because  they  are  essentially  uniaxial  motion  sensors  (or  equivalently 
pressure  gradient  sensors).  Because  the  turbulent  boundary  layer  wall  pressure  is 
homogeneous  across  the  circumference  of  the  array  (e.g.,  there  is  no  preferred  orientation 
or  statistical  variation),  the  “X”  and  “Y”  dipole  hydrophones  should  have  the  same  self¬ 
noise  response. 
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(c)  n  =  2 

FIG.  4.  Examples  of  ratio  of  inside  pressure  to  outside 
pressure  for  infinitely-long  fluid-filled  cylinder. 


The  circumferential  distribution  of  pressure  or  shear  stress  on  the  outside  of  an 
array  can  be  represented  as  an  infinite  sum  of  “modes.”  The  mode  n  pressure  field  can  be 
defined  as, 

Pn(¥)  =  ^iP(¥)e“in¥dt|/,  (2) 


where  P(\)/)  is  the  pressure  at  angle  \j/.  A  useful  semi-empirical  model  of  the  turbulent 
boundary  layer  for  homogenous  flow  over  a  flat  plate  is  given  by  Chase6.  This  model  is 
given  in  the  wavevector-frequency  domain.  For  the  cylindrical  array,  the  spanwise 
wavenumber  at  which  the  mode  n  response  is  centered  is  kz  =  n/a  where  a  is  the  cylinder 
outer  radius.  To  first  order,  the  flat  plate  turbulent  boundary  layer  model  can  be  assumed 
to  hold  for  flow  over  a  cylinder  except,  perhaps  for  very  large  ratios  of  the  boundary 
layer  thickness  to  the  cylinder  radius,  5/a.  The  reader  is  referred  to  Ackroyd7  for 
estimates  of  the  boundary  layer  properties  for  turbulent  flow  over  long  cylinders. 

The  mode  n  wall  pressure  axial  wavenumber-frequency  spectrum  can  be  estimated  from 
the  wavevector-frequency  spectrum  as 

00  ? 

<Fn(kx,co)=  |d>(kx,kz,(o)|Hn(kz)|  akz  ,  (3) 

—  00 

where  Hn  is  the  response  function  for  a  mode  n  sensor  and  is  given  by 


Hn(kz)  = 


sin[(kz  -n/a)rca] 
(kz  -n/a)rca 


(4) 


Figure  5  shows  an  example  of  the  estimated  mode  1  wall  pressure  field  present  on 
the  outside  of  a  7  cm  outer  diameter  cylinder  towed  at  speed  of  20  knots  with  a  140  cm 
boundary  layer  thickness.  The  majority  of  the  energy  from  the  boundary  layer  occurs 
along  the  convective  ridge  (kx  =  co/Uc).  This  relatively  high  wavenumber  energy  is  easily 
filtered  by  the  array  hose  response  and  group  designs.  The  residual  levels  at  low 
wavenumbers  (e.g.,  lkxkco/c)  remains  as  the  limiting  self-noise  (the  “low  wavenumber 
limit”). 


Using  the  mode  n  component  of  the  turbulent  boundary  layer  wall  pressure 
excitation  and  the  ratio  of  the  inner  to  outer  pressure  field  (in  the  wavenumber-frequency 
domain),  the  frequency  autospectrum  of  the  mode  n  response  can  be  estimated.  For  a 
dipole  hydrophone,  the  mode  1  pressure  is  related  to  the  farfield  pressure  (Pg)  as 

|pff|  =  ^|pl|.  <5) 

where  k  =  co/c  and  Pi  is  the  mode  1  pressure  at  the  dipole  hydrophone  radius  a.  Bokde  et 
al.8  estimated  the  mode  n  wall  pressure  from  measured  data  around  the  circumference  of 
a  cylinder  in  axial  flow.  They  noted  the  mode  0  and  mode  1  pressure  spectra  had  similar 
levels.  The  equivalent  farfield  pressure  for  a  dipole  hydrophone,  however,  will 
inherently  be  noisier  than  for  an  omnidirectional  hydrophone  due  to  the  scaling  of  (5). 
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By  nature,  the  dipole  hydrophones  respond  to  the  spatial  gradient  of  the  pressure  field,  so 
it  is  not  unexpected  that  the  dipole  hydrophones  will  have  higher  self-noise  levels  than 
omnidirectional  hydrophones. 
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FIG.  5.  Example  of  predicted  mode  1  wavenumber-frequency 
spectrum  for  external  flow  over  a  cylinder. 


An  example  of  the  predicted  mode  1  self-noise  is  shown  in  Fig.  6  for  tow  speeds  of 
5,  10,  15,  and  20  knots.  This  estimate  is  only  for  the  contributions  from  the  turbulent 
boundary  layer  wall  pressure  to  a  single  hydrophone  (length  5.1  cm  and  diameter  1.4  cm) 
inside  a  7-cm  outer-diameter  flexible,  fluid-filled  cylinder.  In  general,  the  response  is 
dominated  by  excitation  of  the  bending  wave.  At  higher  frequencies,  the  sensor  finite 
length  begins  to  filter  the  high  wavenumber  energy  of  the  bending  wave.  A  longer 
hydrophone  or  a  group  of  hydrophones  has  been  shown  to  further  spatially  filter  this 
wave  energy.  At  low  speeds,  the  dipole  hydrophone  response  may  also  be  limited  by 
cable  strum,  ambient  acoustic  noise,  electronic  noise,  or  a  combination  of  these. 

Experimental  measurements  of  the  dipole  self-noise  in  a  towed  array  also  showed 
the  possibility  of  “leakage”  of  energy  of  higher  order  modes  into  the  mode  1  response. 
Near  the  cut-on  frequency  of  a  lobar  mode,  the  impedance  of  the  cylindrical  structure  is 
very  low,  and  the  high  vibration  levels  can  bleed  into  other  modes  due  to  imperfect 
circular  geometry.  For  example,  elevated  self-noise  levels  around  the  mode  5  cut-on 
frequency  have  been  measured.  These  noise  augmentations  were  not  speed  dependent 
but  were  found  to  decrease  when  the  array  fill  fluid  was  pressurized  to  minimize 
geometric  distortions. 
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Modeling  and  measurements  of  the  performance  of  dipole  hydrophones  in  towed 
arrays  suggest  that  the  self-noise  field  in  a  fluid-filled  cylinder  excited  by  a  turbulent 
boundary  layer  precludes  meaningful  use  at  low  frequencies.  But,  the  exact  frequencies 
of  this  usefulness  depend  on  many  factors  including  primarily  the  tow  speed,  array 
diameter,  and  ambient  noise.  Signal-free  reference  noise  cancellation  has  been 
implemented  using  accelerometers  designed  with  low  acoustic  sensitivity  to  achieve  up  to 
20  dB  of  noise  reduction  at  low  frequency.  To  date,  effective  use  of  dipole  hydrophone 
technology  in  towed  systems  has  been  with  mid-frequency  arrays  at  low  tow  speeds.  The 
choice  of  the  type  of  dipole  sensor  (e.g.  accelerometer,  geophone,  displacement  sensor) 
depends  essentially  on  electronic  performance  and  does  not  have  an  effect  on  the  flow- 
induced  self-noise  limits.  Using  solid-filled  towed  arrays  with  dipole  sensors  is  not 
expected  to  improve  self-noise  because  the  bending  or  shear  wave  would  still 
significantly  contribute  to  the  dipole  hydrophone  response. 


FIG.  6.  Example  of  dipole  hydrophone  self-noise  prediction. 


V.  CONCLUSIONS 

The  design  and  use  of  directional  hydrophones  for  fluid-filled  towed  arrays  was 
discussed.  The  primary  noise  mechanism  for  dipole  hydrophones  was  shown  to  be  a 
bending  or  shear  wave  in  a  fluid-filled  cylinder.  The  turbulent  boundary  layer  flow 
around  the  array  excites  many  modes  of  response  include  the  n  =  0  breathing  wave  and 
the  n  =  1  bending  wave.  The  dipole  hydrophone  responds  primarily  to  the  mode  1 
pressure  field  inside  the  array  (related  to  the  fluid  velocity).  The  self-noise  is  highly 
speed  and  frequency  dependent.  The  self-noise  can  be  minimized  by  using  an  extended 
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hydrophone  to  filter  out  high  wavenumber  energy.  The  dipole  hydrophone  performance 
at  low  frequency  is  fundamentally  limited  by  the  low  wavenumber  energy  present  in  the 
towed  array  system. 
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Abstract- Directional  transducers  will  improve  the 
performance  of  undersea  acoustic  modems.  Current  designs 
provide  only  vertical  directionality  with  omni-directionality  in 
the  horizontal  plane.  Horizontal  directionality  can  yield  a  higher 
source  level  and  higher  received  signal  to  noise  through  the 
increase  in  the  DI  or  equivalently  the  narrower  beam  pattern  of 
the  device.  Moreover,  a  horizontal  steered  directive  beam  allows 
selective  and  more  secure  communications  between  the  nodes  of 
the  system.  This  paper  reports  the  development  of  a  low 
frequency  prototype  device  that  achieves  desired  directivity  by 
the  controlled  combining  of  the  three  most  fundamental 
extensional  modes  of  vibration  of  a  piezoelectric  cylinder. 

I.  Introduction 

The  multimode  directional  telesonar  transducer  provides 
horizontal  directionality  through  extensional  modal  addition. 
Results  on  a  high  frequency,  3.5  inch  diameter  telesonar  were 
presented  earlier1.  The  present  transducer  is  a  vertical  array 
which  is  composed  of  a  larger  quadrupole  33  mode  staved 
cylinder,  a  striped  tangentially  poled  dipole  cylinder  and  a 
pair  of  31  mode  omni  cylinders.  By  combining  these  three 
modes,  omni-directional,  dipole  and  quadrupole,  horizontal 
directionality  is  excited  through  each  rings  extensional 
mode2.  With  this  method  identical  beam  patterns  may  be 
obtained  over  a  range  of  frequencies.  In  addition  to 
horizontal  directionality  the  beam  can  be  electrically  steered. 
It  has  been  shown  that  highly  directive  beams  may  be 
obtained  from  spherical  transducers3.  We  extend  this  to  the 
case  of  a  cylindrical  radiator.  Analytical  and  finite  element 
modeling  were  performed  to  verify  expectations  before  the 
design  was  implemented  to  achieve  cardioid  and  super 
cardioid  horizontal  directional  patterns.  Measured  results  are 
presented  and  are  shown  to  agree  reasonably  well  with 
calculated  results. 

II.  TRANSDUCER  OPERATION 

A.  Ring  Modes 

We  consider  first  the  transducers  modes  of  operation  and 
resonant  frequencies.  Figure  1  shows  the  first  three 
extensional  modes  of  operation.  With  f0  the  resonant 
frequency  of  the  "breathing"  omni-mode,  the  resonant 
frequencies  of  the  higher  order  extensional  modes  are  given 
by  fn  =  f0(l  +  n2)1/2 , -where  n  is  the  mode  number  and 


Fig.  1  First  three  extensional  ring  modes 


Fig.  2  Image  Acoustics,  Inc.  Multimode  Telesonar 
Directional  Transducer  (3  ring  design) 

f0  =  c/7lD  with  c  the  sound  speed  in  the  ring  of  mean  diameter 
D.  The  mode  numbers  n  =  0,  1  and  2  give  the  frequencies  of 
the  omni,  dipole  and  quadrupole  modes  respectively;  thus,  fi 
=  f0V2  and  f2  =  f0V5. 

Figure  2  shows  the  prototype  multimode  transducer.  A 
resonant  frequency  of  10.5  kHz  was  chosen  for  each  ring.  A 
diameter  of  4.5  inches  was  chosen  for  the  omni  mode,  a 
diameter  of  5.45  inches  for  the  dipole  ring  and  a  diameter  of 
9  inches  for  the  quad  ring.  Since  this  transducer  diameter 
was  limited  to  6.5”  the  quad  ring  was  too  large.  To  decrease 
the  quad  diameter,  twelve  glass-reinforced  G-10  wedges 
along  with  twelve  pair  of  33  mode  piezoelectric  bars  were  use 
to  construct  the  quad  ring.  The  G-10  is  stable  with  a  lower 
sound  speed  than  PZT  and  was  chosen  to  yield  a  reasonably 
small  sized  6.2  inch  diameter  quad  mode  cylinder  with  a 
resonant  frequency  in  the  vicinity  of  10.5  kHz.  The  smaller 
diameter  10.5  kHz  striped  tangentially  polled  dipole  ring 
yields  high  coupling  and  a  12  pair  electrode  segmentation 
corresponding  to  the  quad  ring  allowing  steering  every 
360/12  =  30°.  The  segmentation  allows  the  excitation  of  the 
quad  ring  into  the  quad  mode  of  vibration  by  reversing  the 
phase  in  each  successive  quadrant.  The  dipole  mode  is 
achieved  by  reversing  the  phase  on  one  half  of  the  cylinder. 
The  smaller  diameter  omni  mode  cylinders  are  resonant  at 
nearly  the  same  frequency  and  needs  no  segmentation  for 
their  uniform  radial  mode  motion.  Omni  directional  receiving 
response  is  obtained  from  the  omni  cylinders. 


(3) 


B.  Multimode  Approach 


In  the  two-mode  method  the  omni  mode  has  a  directivity  of 
1  while  the  dipole  mode  has  a  figure-eight  directivity 
function  cos  0  where  0  is  the  cylindrical  coordinate  azimuth 
angle  in  the  horizontal  plane.  If  we  simply  add  the  two,  and 
normalize,  we  get 


P(0)  = 


1  +  cos  6 

2 


(1) 


which  is  a  cardioid  beam  pattern  function  with  -3dB  value  at 
±  65.5°  and  a  deep  back  null.  This  may  be  better  visualized  in 
Fig.  2.  Here  the  cardioid  is  formed  by  adding  the  omni  and 
dipole  modes.  Both  have  unity  amplitude.  Since  the  “+”  half 
of  the  dipole  is  in  phase  with  the  corresponding  omni  half, 
they  both  add,  doubling  the  amplitude.  The  half  of  the 
dipole  is  180°  out  of  phase  with  the  corresponding  omni  half 
causing  cancellation  yielding  the  cardioid  pattern. 

In  the  three-mode  method  the  quadrupole  pattern  function 
may  be  written  as  cos  20.  If  we  combine  this  with  an  omni  of 
strength  1  and  a  dipole  of  strength  3,  we  get  the  "super 
directive"  beam  pattern  function 


F(Q)  =  I  +  2cos0 +  cos20  (2) 

To  better  visualize  the  three-mode  method,  Fig.  4  shows 
the  three-mode  pattern  synthesis.  First  the  strength  1  omni 
mode  pattern  of  (a)  is  added  to  the  strength  1  quadrupole 
pattern  (b)  to  get  the  bidirectional  pattern  (c)  of  strength  2 
with  lobes  of  the  same  phase.  If  this  is  now  combined  with  a 
strength  2  dipole  pattern  (d),  we  get  the  desired  resulting 
"super  directive"  pattern  (e)  of  strength  4  which  conforms  to 
Eq.  (2).  The  result  may  be  normalized  by  dividing  by  4  to 
yield  a  unity  value  at  0  =  0.  We  also  can  use  trigonometric 
identities  to  get  the  normalized  result  shown  in  Eq.  (3). 


GXD 


Dipole 


Cardioid 


l  +  2cos0  +  cos(20)  cos(0)[l  +  cos(0)] 


As  seen,  the  beam  pattern  function  for  this  “unique” 
distribution  is  the  product  of  a  dipole  and  a  cardioid  function. 

Adjusting  the  amplitude  of  the  quadrupole  to  0.414,  while 
leaving  an  amplitude  of  1  for  the  omni  and  dipole  modes, 
results  in  a  “minimalist”  super  cardioid.  We  call  this  a 
minimalist  super  cardioid,  since  it  achieves  a  desirable  90°  - 
3dB  beam  pattern  and  15  dB  front-to-back  ratio  with  minimal 
dependence  on  the  higher  modes.  Equation  (3)  can  be 
rewritten  in  a  general  form  as 


P(0)  = 


1  +  Acos0  +  2?cos(20) 
1 +  A  +  B 


(4) 


For  the  case  of  the  cardioid  A=1  and  B=0.  For  the  unique 
super  cardioid  A=2  and  B=1  and  for  the  minimalist  super 
cardioid  A=1  and  B=0.414,  our  current  case  of  interest. 

III.  MEASURED  RESULTS 

The  measured  Transmitting  Voltage  Response,  TVR,  for 
the  omni,  dipole  and  quad  sections  driven  separately  are 
shown  in  Fig.  5.  The  measurements  on  the  dipole  and  quad 
units  were  made,  for  convenience,  with  half  of  each 
transducer  driven  in  series  with  the  second  half.  The  TVR 
results  would  be  6  dB  higher  if  the  sections  were  all  driven  in 
parallel.  A  comparison  with  theoretical4,5  and  finite  element 
models  indicates  that  the  measured  omni  response  is  very 
close  to  prediction.  The  dipole  resonance  is  approximately 
0.5  kHz  higher  than  expected  with  reduced  output  below 
resonance.  The  quad  resonance  is  approximately  1.0  kHz 
higher  than  expected.  (The  first  experimental  quad  mode  ring 
yielded  a  resonance  closer  to  the  calculated  value.)  The 
receiving  free  field  voltage  sensitivities  for  each  ring  are 
shown  in  Fig.  6.  The  response  of  the  omni  ring  falls  off  above 
resonance  and  is  reasonably  flat  below  resonance  as 
expected.  These  resonant  frequencies  are  higher  than  the 
TVR  values  because  of  the  open  circuit  conditions.  All  the 
sensitivities  are  above  -190  dB  in  the  range  from  8  kHz  to  13 
kHz. 


Fig.  3.  Two-mode  beam  synthesis 
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Fig.  4.  Three-mode  pattern  synthesis 


Frequency  (kHz) 

Fig.  5  Multimode  Telesonar  Directional  Transducer  (3 
ring  design)  transmitting  voltage  response 
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Fig.  6  Multimode  Telesonar  Directional  Transducer  (3 
ring  design)  receiving  response 

The  horizontal  measured  beam  patterns,  at  10.5  kHz,  for 
the  dipole  and  quad  modes  are  shown  in  Figs.  7  and  8.  As 
seen,  the  dipole  and  quad  nulls  vary  from  -15  dB  to  an 
exceptional  -30  dB  values.  The  lobe  levels  are  essentially  the 
same  for  the  dipole  case  and  nearly  the  same  for  the  quad. 
The  horizontal  beam  pattern  results  for  summing  the  omni 
with  amplitude  1,  dipole  with  amplitude  1  and  quad  with 
amplitude  0.414  are  shown  in  Fig.  9.  As  seen  These  results 
are  very  close  to  the  ideal  levels  of  -3  dB  at  45°,  -12  dB  at 
90°  and  -15  dB  at  180°  given  by  Eq.  (4)  with  A=1  and 
B=0.414. 


Fig.  8  Multimode  Directional  Telesonar  Transducer 
Quad  cylinder  horizontal  beam  pattern  at  10.5  kHz 
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Fig.  9  Multimode  Directional  Telesonar  Minimalist  Super 
Fig.  7  Multimode  Directional  Telesonar  Transducer  Cardioid  horizontal  beam  patterns  (10  dB/div) 

Dipole  cylinder  horizontal  beam  pattern  at  10.5  kHz 


IV.  SUMMARY  AND  CONCLUSIONS 


We  have  presented  results  that  have  demonstrated  the 
feasibility  of  achieving  a  constant  beamwidth  horizontal 
directional  response  using  a  combination  of  the  omni,  dipole 
and  quadrupole  extensional  modes  of  vibration  of  co-resonant 
cylindrical  transducers.  The  measured  results  substantiated 
predicted  results.  Over  the  frequency  range  of  8  to  13  kHz 
the  beamwidths  were  90°  with  -10  dB  at  90°  and  a  front-to- 
back  ratio  of  15  dB. 
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Directional  Hydrophones  for  Towed  Arrays 
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Abstract  The  poster  describes  the  development  of  different  hydrophones 
capable  of  measuring  the  vector  component  of  the  acoustic  wave,  thus 
providing  directional  information  orthogonal  to  the  axis  of  the  towed 
array.  The  preferred  solution  is  the  DERA  Low  Profile  Accelerometer 
(LPA)  hydrophone,  which  comprises  a  biaxial  accelerometer  to  form  two 
orthogonal  dipoles  and  two  omnidirectional  hydrophones  wired  in  parallel 
to  give  a  single  high  capacitance  output.  When  matched  in  phase  and 
amplitude  these  sensors  can  provide  a  steerable  unidirectional  (cardioid) 
output.  Sea  trials  have  shown  that  these  hydrophones  can  provide 
instantaneous  resolution  of  left/  right  ambiguities  and  offer  increased  gain 
with  respect  to  isotropic  noise.  This  work  was  funded  by  the  Ministry  of 
Defence  Corporate  Research  Programme  contract  number  TGI  2.2/0 12/u. 


Executive  Summary 

A  number  of  hydrophones  have  been  developed  by  DERA  under  the  MoD 
sponsored  Sonar  Transducer  Technology  programme.  These  hydrophones  have  an 
inherent  directivity  despite  their  small  dimensions  and  typically  consist  of  2  orthogonal 
sensors  to  measure  the  vector  component  of  the  acoustic  wave,  together  with  a 
conventional  pressure  hydrophone.  When  suitably  matched  in  phase  and  amplitude  a 
unidirectional  (cardioid)  beampattern  can  be  obtained  by  combining  the  dipole  output  of 
the  vector  sensor  with  the  monopole  output  of  the  pressure  hydrophone.  A  key  challenge 
is  that  these  sensors  must  be  compact  (less  than  30mm  diameter),  while  achieving  the 
required  directivity  and  sensitivity,  against  a  background  of  high  vibrational  noise,  within 
the  towed  array.  They  must  withstand  hydrostatic  pressures  of  100  MPa  (equivalent  to  an 
ocean  depth  of  1000m)  while  sensing  dynamic  pressures  as  small  as  100  jaPa. 

Geophone,  accelerometer  and  pressure  gradient  designs  have  been  developed  and 
have  been  extensively  tested  both  in  the  laboratory  and  at  sea. 

Traditional  pressure  gradient  designs  offer  a  poor  low  frequency  sensitivity,  are 
susceptible  to  noise  and  are  unsuitable  for  small  diameter  arrays.  The  geophone-based 
hydrophones  performed  well  especially  at  lower  frequencies  and  were  the  lowest  cost  at 
approximately  $150  per  unit  using  Geospace  GS1L9  geophones.  The  accelerometer- 
based  hydrophones  offered  the  best  bandwidth/  sensitivity  product  and  were  the  preferred 
solution  especially  for  use  as  receivers  in  an  active  sonar  system.  Commercial  devices 


were  too  large  and  expensive  and  a  Low  Profile  Accelerometer  (LPA)  hydrophone  was 
developed  by  DERA. 

A  Vector  Towed  Array  was  developed  to  demonstrate  the  ability  of  an  array  of 
directional  hydrophones  to  resolve  the  bearing  ambiguity  and  to  establish  their  self-noise 
response.  The  array  consisted  of  32  LPA  hydrophones  each  containing  a  pressure 
hydrophone  and  a  pair  of  orthoganal  accelerometer  outputs. 

This  array  was  evaluated  during  a  trial  at  the  Naval  Surface  Warfare  Centre  facility 
in  Idaho,  USA.  The  self-noise  response  of  vector  hydrophones  is  critical  to  their 
application.  A  detailed  analysis  of  their  self-noise  performance  is  beyond  the  scope  of 
this  paper,  however  analysis  of  the  trials  data  show  that  the  array  is  clearly  able  to  resolve 
both  port  and  starboard  contacts  instantaneously  and  under  realistic  operational 
conditions.  The  beamformed  array  resolves  the  target  ambiguity  with  discrimination 
between  port  and  starboard  targets  of  up  to  20  dB 
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Beamformed  response  to  a  starboard  target  of  16  LPA  hydrophones  in 

omni  and  cardioid  modes 


The  DERA  LPA  hydrophone  is  the  smallest,  high-performance  directional  sensor 
tested  under  this  programme  and  its  ability  to  instantaneously  resolve  left  right 
ambiguities  has  been  successfully  demonstrated  at  sea.  The  ability  to  resolve  bearing 
ambiguities  could  provide  the  Navy  with  an  obvious  operational  advantage  particularly  in 
shallow  water  where  manoeuvring  constraints  may  be  greater. 
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